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ABSTMCT 

The  second  generation  iesign  of  ttis  Naval  Postgriduata 
School's  ocean  floor  geomagnetic  data  collectiorj  systeiu  is 
described  with  emphasis  on  modif  ici  tions  for  i[nprcvsrns',.t 
over  its  predecessor.  These  improveasat s  include  15  channel 
parallel  analog-to-serial  palse  codB  oiodulation  conversion, 
fiber  optic  data  link,  and  raiiD  telemetry  to  shor'^ 
recording  eguipment.  The  system  has  flexibility  in  operating 
depth,  an  increase  in  data  icquisition  time,  syn::hronization 
with  land  site  data,  and  a  lata  fDrniat  that  is  readily 
converted   to   digital   for   co3iput=r  assisred   analysis. 
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This  thesis  resedrch  is  part  of  =  continuing  effort  by 
the  Naval  Postqrac.aate  School  to  obtain  a  long  term  data 
base  for  the  improved  unis rsti ndiag  and  interpretation  of 
the  ELF  electrotp.a  gn^t  ic  environment  ii  the  vicinity  of  tne 
sea  fleer.  The  overall  project  emphasizes  the  importance  of 
acquiring  measarements  of  geomagnetic  fluctuations  on  the 
sea  floor  over  a  period  of  sevsral  years  ht  various 
locations  and  depths  while  simultaneously  obtaining  data  at 
a  land  based  sit -3.  Primary  obj  actives  include  the 
inrerpretation  cf  sigre.ls  covering  foar  decades  of  frequency 
froi  O.OIHZ  to  100KZ  through  th2  use  of  total  field 
magnetometers,  induction  coils  and  (JLF/ELF  receivers.  This 
phase  cf  data  collection  covers  ttis  O.OIHZ  to  20HZ  range 
utilizing  induction  coils  :i  s  seasors. 

The  magnetic  noise  in  t a e  sea  is  Df  interest  both  from  a 
geophysical  viewpoint  :-3  w=ll  as  for  laval  applications.  In 
addition  to  undersea  field  ma gneton 3t ry ,  sone  areas  of 
current  geophysical  intarrst  are  th9  leasurement  of  marine 
geomagnetic  anoraalirjs  neaz  centers  of  sea  floor  spreading, 
characterization  of  magnatic  fields  induced  by  ocean  waves, 
geophysical  exploration  and  the  utilization  of  low  frequency 
waveguides  present  uni'er  the  sea  floor.  Applications  of 
interest  to  the  Navy  ara  in  tha  areas  of  mine-warfare, 
submarine  detection,  sibmarine  cd nmunication  using  a 
suparconduct ing  quantum  intarfecence  ievice  (SQUID)  as  the 
sensing  element  for  an  extremely  low  frequency  (ELF) 
recaiving  antenna  system,  and  possibly  as  a  means  of  remote 
verification  of  underground  nuclear  datonations. 
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To  achieve  an  uE:inr3t  andL  ng  of  the  natirs  cf  -the 
fluctuating  under  sea  geomagnetic  s;pectrum,  it  is  necessary 
to  develop  a  data  base  witi  goals  for  high  quality  and  great 
quantity.  The  collection  s/otsm  utilizBi  by  the  predecessors 
[  1  ]  of  this  research  team  -Acknowledged  a  shortfall  in  both 
these  goals.  Their  collection  systsm  consisted  of  two 
channels,  each  with  a  coll  antenna  ssnsor,  preamplifier, 
variable  gain  amplifier,  voltage  controlled  oscillator, 
reference  oscillator,  and  xn  anc.log  cassette  tape  recorder 
configured  as  depicted  in  Figure  1.1,  The  components  were 
encapsulated  in  Benthos  glass  s-pheres  with  the  deployed 
configuration  as  shown  in  figure  1.2.  The  data  collected 
with  this  system  was  limited  by  virtue  of  the  equipment 
restrictions  and  deployoaent  techniques.  The  analog  tape 
recorder  posed  the  greatest  limitetioi.  The  data  acquisition 
time  was  limited  by  the  cassette  tape  length  to  a  maximum  of 
45  minutes  per  channel  per  system  deployment.  The  accuracy 
of  the  data  was  greatly  iepenceat  apDQ  the  stability  of  the 
oscillators.  rtdditionally  the  utility  for  comparison  of  the 
data  with  data  taken  ai  a  nearby  land  station  is 
questionable  since  the  synchronization  of  the  data  could 
only  be  crudely  approxisa t;3i  .. 

The  particular  objectives  of  this  thesis  research  was  to 
improve  the  ocean  floor  geoiiagnetic  collection  system  design 
in  an  effort  to  improve  data  quality  and  quantity  while 
attaining  time  synchronization  of  land-oceaa  data.  A 
secondary  objective  was  to  develop  procedures  for  computer 
assisted  data  analysis.  Lastly  the  associated  thesis  report 
is  to  act  as  a  systems  nanual  in  an  effort  to  provide 
continuity  for  future  systea  users. 

In  the  design  of  an  ocean  floor  lagnetic  data  collection 
system,  it  is  necessary  zd  consider  the  constraints  placed 
upon  the  system  as  a  result  of  its  function  and  the 
environment  to  which  it  is  subjected. 
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Figure   1.1         Data Ac^aisition   Systam   Schematic 
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Figure    1.2        Diagraa  of   OeplDyed   System   Configuration 
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It  is  the  function  of  th^  syotsu  to  sense  gecmagnetic 
fluctuations  near  the  ccsan  floor  in  the  frequency  range  of 
intirest,  process  and  tcansmi"  those  signals  to  the  surface 
for  telemetry  to  shore  basad  receiving  equipment  for  dual 
(ocaan-iand)  r^icording,  dec::diag,  and  analysis  while 
maintaining  the  information  content  of  the  signals.  One 
requirement  of  the  integrated  systen  is  to  ensure  signal 
purity.  That  is,  zhe  system  3hoal3  measure  the  magnetic 
signals  while  minimizing,  to  tha  sxtant  possible,  system 
induced  variations  in  the  signal  reiag  measured.  It  is  well 
known  that  a  magnet,  ic  fieii  cm  b=  set  up  by  an  electric 
field  (iimpere's  Law).  This  -nutual  interference  between 
system  components  must  oa  avoids^  where  possible  and 
minimized  otherwise.  To  a:;hiev3  aiaxii  um  signal  purity,  the 
sensors  must  be  isolated  magnetically  as  much  as  possible 
fron  the  rest  of  the  syst=?i.  BecaLse  of  ths  low  signal 
strsngth  of  the  collectei  iit'i,  seas  =  iplif icatioa  must  take 
place  at  the  sensor  location  to  all:'W  for  transmission  by 
nonferrous        metallic        caoLe        tr  a        remotely        located 

instrumentation  subsystem,  7arther  maza^tic  isolation  may  be 
achieved  by  utilization  of  a  fiber  of'tic  data  link  from  the 
instrumentation  sphere  to  tie  :ilem2t:y  bouy.  Additionally, 
near  surface  interference  is  ^  aatter  of  concern.  Low 
frequency  electromagnetic  cadiation  ca^i  penetrate  the  ocean 
surface  for  up  to  several  meters  with  the  potential  of 
channelling  down  the  link  aid  providing  perturbations  to  the 
signal  being  measured  if  a  sstallic  conductor  cable  were 
used  for  data  tr ansmissioi.  This  "antenna  effect"  can  be 
eliminated   by   employment   of    a   fiber   optic    link. 

Due  to  space  and  power  limitatiois,  transmission  line 
signal  attenuation  may  be  ai  important  consideration  in  deep 
sea  locations.  Low  loss  fibsr  optics  ^an  reduce  Dr  eliminate 
this    concern. 
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The  advantage  of  f  iberr-opti.T  links  have  been 
demonstrated  world-wide  by  expariaeatiil  optical  wave  guide 
systems.  Among  the  importiat  fsaturrss  are  low  transmission 
loss,         wide        transmission      bandwidth,  insensi tivity      to 

electromagnetic      and      raii3         frequrnry      int'rrf  erence,        and 
suitability    for   digital  ^omm uni-ati~n3 ,    and    oulsa    modulation 
methods    (for      amplitude   moiiilation   1 5::haique:5,         fiber    optic 
cable      losses   are      indepenisnt      of   tr?. 'i^rjiss  ion        frequency) 
[2]. 

Optical  fibers  are  attractive  sii^':^  a  single  strand  is 
sufficient  to  carry  the  ii formation  signal  and  cannot  be 
short-circuited  by  exposure  to  si-a  vater.  Mso  recent 
treads,  through  advanced  technical  expe;:ti32  in  fiber 
manufacturing  have  narrowed  the  e::oaorp.LC  g.ip  between  long 
haul  optical  fibers  and  their  for-iriost  competitor,  copper 
wirs    and   cables   [  3]. 

Any  practical  desiga  ultimately  r?st5  upon  the 
requirements  of  the  user,  ais  resourrrs  and  capabilities  as 
well  as  the  technology  of  the  day.  The  major  goals  have 
been  stared  as  that  of  signal  purity  "ini  data  quantity.  It 
has  been  shown,  although  not  rigorously,  that  a  fiber  optic 
data  link  is  a  viable  alternative  to  the  metallic 
transmission  line  in  the  a^aievsment  or  signal  purity.  The 
improvement  of  data  acqui^itioQ  tias  i~  gained  through  RF 
transmission  of  data  to  shore  r5c:2LviD.g  equipaent.  Data 
acquisition  time  is  therebf  limi-ei  only  by  power  supply 
capacity    and   the    schedula    of    the    deploying    ship. 

As  indicated  earlier,  fiber  optic  ^^ystems  are  suitable 
for  digital  and  pulse  modulation  methods  of  data 
transmission.  In  addition,  inalog  transmission  may  be 
accomplished.  There  ace  several  pulse  modulation  methods 
available.  Among  them  are  pulse  position  modulation  (PPM), 
pulse    width      modulation    (PiJi),         pulss    amplitude      modulation 
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(PAM) ,  and  pulse  code  moiilation  (PCS) ,  ?CM  is  ths  best 
known,  and  today  the  most  important  palsa  modulation  system. 
Thrse  separate  operations  are  involved  in  providj.ng  a  PCM 
signal.  The  first  is  to  interrogate  the  inessage  signal  at 
regularly  spaced  intervals  (sampling).  The  sscc^nd  is  to 
approximate  the  measured  amplituda  value  to  th??  nearest 
periitted  voltage  referen::3  level  (gaantizing^  .  The  third 
operation  is  to  represeit  the  aporoxiniatec  (quantized) 
amplitude    values    as    a    series    of    cod=d    pjlses. 

In  PCM,  several  pulsss  per  sampla  are  used  t.o  signify 
the  amplitude  value,  instaad  of  one  pulse  per  sanple  as  in 
the  cases  of  PAM,  PPM,  or  PWM.  Consequently^  PCM  systems 
require  an  increased  bandwidth  (witiin  the  capability  of 
fibsr  optic  systems)  for  their  transmission.  HDws:ver,  any 
small  deformations  in  ths  hsight  or  rfidth  of  tna  nulses  are 
irrslevant  since  it  is  only  necessary  to  know  wiiether  the 
pulse  is  present  or  absent  in  order  to  r  =  tri'iVe-  th^  original 
message.  Moreover,  in  a  PCM  traasmission,  noise  is 
nonaccumulative  because  noisy  PCM  signals  caa  ^rasily  be 
cleaned  up,  when  it  beco^nss  nscessacy,  by  "he  -process  of 
regeneration.  Thus  the  qualizy  of  a  PCM  tr'insai.ssion  is 
dependent  on  the  sampling,  quantizing,  arid  cc-.ling  ])rocesses, 
and  not  the  length  nor  the  noise  of  zh-  trans;Tiis3:.on  media. 
By  contrast,  the  FAM,  PWM,  and  PPM  s^s^ems  are  continuously 
affected  by  noise  and  cannDc  be  cleaned  up,  o'^  regenerated. 
The  noise  is  accumulative  and  the  Longer  th'3  distance  of 
transmission,    the   greater    will    oe   ths    noise   [ '+ ]• 

Additionally,  the  desired  resolition  of  the  measured 
signals  may  be  achieved  in  a  PCM  systam  by  selection  of  the 
sampled  word   length  and  ths    sampling    rate. 

Both  the  digital  and  analDg  transmission  aethods  are 
inferior  to  PCM.  The  loss  of  a  bit  of  digital  data  obscures 
the      value      of      the      samplsi      data.  Analog      data      may      be 
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irrstrievably  altered    by   riDise.      PC^.^       '■-hen,      is    perliaps   rhe 
ideal   transmission  method    for   our   2.ppl  icdtion. 

This  report  delineates  an  ocean  floor  geoa5gn'?*.ic  daza 
collection  system  that  employs  coil  antenna  sensors,  a  pulse 
cods  modulation  system,  fibsr  optic  link,  RF  dati  :.i.:\k,  and 
recording  equipment.  It  also  details  prDcedur^s  tnu  decoding 
and  digitizing  the  data  as  well  thoss  for  compacc-r  assisted 
data    processing   for   eventual    analysis. 
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II.       SrSTEM    DESISN 

A.       EQUIPMENT    CON  FIGU RATI OIJ 

The     general    arrangemsat      cf      the      SensratiDa    II      3yst?:;n 
components   in   their    deployBl    mode    is    shown    in   Figure    2.1, 

"*•      Data   Collection  SaiLMiSi 

A  factional  block  diagram  of  the  data  agaisition 
system  illustrated  in  Figurs  2.2  revails  the  follDwing  ni?jor 
components: 

1.  coil   antenna   sensors    (2) 

2.  preamplifiers    (2| 

3.  preconditioners    (2) 

U.  pulse    code    moiulition    systaa     (1) 

5.  optical  transmitt er  (1 ) 

6.  optical  receiver     (1) 

7.  radio    frequency    transmitter     (1) 

8.  radio    frequency   receiver    (1} 

S.    instrumentation   tape   recorisr    (1) 

The  system  componerxts  are  discussed  in  th=  following 
subsections: 

a-       Sensor  Subsystea 

The  sensor  subsystem  is  currently  assigned  with 
two  sensors  sach  consisting  of  a  c;Dil  antenna  sensor,  a 
preamplifier  and  a  battery  power  supply  housed  ii  a  Benthos 
glass  sphere  with  a  0.404  aster  inner  diameter.  The  sensor 
and  the  instrumentation  suQsystems  ar=  interconnected  by  two 
30  aeter  long,  coaxial  cables  teraiiated  at  each  end  with 
Brantner  connectors  and  Bentos  sphere  penetratDrs  through 
the    glass    walls    to  the   enclosed   equipment. 
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The 


sensor   spheres   ars   mountei    in   Bsnrhos    hsi'lhat 


+  ,-> 


to      tne      equlpiI!9I^ 


enclosures.        The      hardhats      are      bolted 
mounting   assembly   with  nylon    bolts. 

(1)       Equipment      Mo uq tin 3    i£§.§S2l2      The    equip!r.'?^nt 
moulting    assembly   depictad      in   Figurs    2.3    is      a    lon-aeta ilic 


3/4"   Fiber-Glass 
Mounting  Plates 


Free  Flooding 
ABS   Structure 


Figure  2.3        Sensor    Subsystem   Squipmsnt    Hounting    assembly 


stand  that  permits  ease  of  ieploymsnt  and  orovidss 
structural  support  and  protection  fDr  the  senscr  spheref;* 
The  sensor  stand  design  calls  for  a  light,  frse  flooding 
stricture  with  reraovabls  strap  on  weights  to  obtain  the 
required  negative  buoyancy,  but  still  allow  for  portability 
and  ease  of  handling  in  system  depl^yient.  A  aon-metallic 
material    is    necessary    in    ths    stand   construction    to    avoid   tha 


enacts 


or 


induced     currants 
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electro-chemical 
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reactions  at  the  interfa^s  of  a  natal  surface  and  the 
eloc  urically  conductive  seai/ater.  Ths  stand  is  cons-ructei 
of  light  weight,  high  strength, 

Aeryionitrils-But aoliene-St yrene  (ABSI  piping  which  is  well 
suited  for  this  application.  The  weights  are  fashioned  from 
lengths  of  ABS  piping  filled  with  leai,  sand  and  transformer 
oil.  The  pipes  are  capped  off  at  each  end  to  prD/ide  a  water 
tight,  non-metallic,  aai  compact  means  of  adding  the 
necessary   negative  buoyancy    to   the    stind. 

(2)  Q2II.  Antenna  SaisD.r  Each  sensor  is  a 
continuously  wound  coil  antenna  manufactured  from  5460  turns 
of  ^S  gauge  copper  magnet  wire  by  Elna  Engineering  of  Palo 
Alto,  California.  The  coils  weigh  approximately  100  pounds 
eacn  with  dimensions  as  depicted  in  Figure  2.U.  The  average 
enclosed  area  of  each  coil  is  3.0321  square  meters  which  is 
determined  by  taking  the  average  of  the  areas  obtained  by 
coj)patinq  the  area  for  both  the  inner  radius  and  the  outer 
radius.  The  dimensions  of  the  sensor  ire  constrained  by  the 
geci.etry  of  the  glass  spheres.  The  ooil  resistance  is  120 
oh^Aj    with   a    self  -  inductance    of   approximately    9.31    henries. 

(3)  Pream^li f 1 er    The      preamplifier    selected    for 

0 

use  in  the  sensor  sphere  Ls  the  moiel  13-10A  low  noise  ELF 
amplifier  manufactured  by  Dr.  Alan  Phillips  of  SRI 
International.  The  final  stage  of  the  amplifier  contains  an 
active  low-pass  filter  with  a  cutoff  frequency  at  20  Hz. 
Such  filtering  is  necessary  to  reject  the  effects  of  60  Hz 
signals  and  harmonics  due  to  the  presence  of  power  lines  for 
the  similar  system  used  for  land  data  acquisition.  Stray 
frequencies  above  20  Hz  do  not  penetrate  significantly  to 
depths  at  which  the  ocean  system  is  operated  (50-50  meters). 
The  overall  preamplifier  gain  is  approximately  55  d5. 

A  schematic  diagram  of  the  preamplifier  circuit  is  contained 
in  Figure  A.I.     The  gain  and  filtering   characteristics  of 
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Figure  2. U    InductiDn  Coil  Seasor  Dimensions 

the  system  are  illustrated  in  Appendix  C.  and  ace  relevant 
to  the  discussion  conceraiig  system  transfer  func-ions.  The 
preamplifier  is  powered  by  two  .75  AH,  12  volt  mercury 
batteries  mounted  in  a  section  of  1  iich  PVC.  The  bart-eries 
are  arranged  in  such  a  fasiion  as  to  give  the  required  plus 
or  ainus  12  vol-s  and  grouii.  The  preamplifier  is  housed  in 
an  aluminum  case  to  shield  in  fields  developed  by  the 
preamplifier  circuitry.  lie  preamplifier  and  bat-ery  pacl< 
are  attached  to  either  side  of  the  CDil  antenna  sensor  with 
electrical  tape.  The  inter  connections  are  made  by  coaxial 
connection  for  both  the  coil  input  to  the  preaiplif ier  and 
the  preamplifier  output  to  the  Bentios  sphere  penetrators 
while  the  power  supply  input  is  made  with  a  keyed 
three-prong  Jones  connector. 

Power  drain  on  various   cDnponents  of  the  system 
is   a    serious   consideration    when   trying    zo      maximize 
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on-3tation  data  collection  tim=.  It  is  also  desirable  to 
have  the  system  sealed  and  made  ready  for  sea  prior  to 
transportation  to  the  ship.  Thas,  an  externally  actuated 
magnetic  activation  circuit  is  employed.  This  circuit  shown 
in  Figure  A. 2  employs  a  normally  open  reed  switch.  The  reed 
switch  affixed  to  the  top  inside  of  the  glass  sphere  is 
closed  when  a  magnet  is  pla::ed  in  its  vicinity  from  outside 
the  sphere.  The  closed  circuit  biases  transistor  switches 
closed  cutting  off  power  to  the  preamplifier.  This 
arrangement  is  advantageous  as  it  allows  packaging  in 
advance  of  the  of  the  deployment  dats.  It  should  be  noted 
however,  that  this  circuit  does  deplste  a  nominal  amount  of 
power;  so  final  packaging  is  usually  completed  one  day  in 
advance  of  deployment. 

b.   Instrumentation  Subsystsi 

The  function  of  the  instr i mentation  subsystem  is 
to  receive  the  two  analog  outputs  from  the  preamplifiers, 
coniition  them  to  the  requirements  of  the  follow  on 
circuitry,  multiplex  the  channelized  signals,  convert  them 
from  analog  to  digital  format  and  pulse  code  modulate  them 
for  input  to  the  fiber  optic  data  iiiK. 

C)  l21ii2Il§-lL  ll2iiQ.Ii^3  ^il.i.l^ll  T^^  equipment 
mounting  assembly  is  shown  in  Figure  2.5.  Constructed  of 
aluminum  it  supports  two  Benthos  glass  spheres,  housed  in 
Benthos  hardhat  enclosures.  One  sphers  houses  the  instrumen- 
tation subsystem  while  the  other  has  a  battery  pack  to  power 
the  optical  transmitter.  En  additioi,  the  optical  trans- 
mitter assembly  is  attached  to  the  vertical  member  of  the 
mounting  assembly.  The  design  incorporates  quick  release 
clanps  for  rapid  engagemeit  or  disengagement  from  the  main 
support  line  (250  feet  of  1  inch  polypropylene  line) . 
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Figure   2.5 


InstrumentatiDi    Subsystaa    Equipment    llcunting 
Assembly 


(2)  Signal  CoQ^i  itioner  Vhs  signal  ronditionsrs 
recaive  the  analog  signals  from  zhs  ooil  preaaplif iers  in 
parallel  fashion,  amplifias  thsm  on  ths  order  of  30  d3  and 
limits  signals  with  peak  amplitudes  greater  thari  7.5  volts. 
The   circuit    is   presented    in    Figure    \.3. 

(3)  Pulse  Code  lodulation  (PCM)  System  The  pulse 
cods  modulation  system  chosen  for  use  in  the  ocean  floor 
geonagnetic  data  collection  system  is  one  designed  and 
manufactered  by  Dr.  Robert  Lowe,  Lowecom  Ino.  ,  also  of 
Scripp's  Institute  of  Oceaao graphy.  La  Jolla,  California. 
The  device,  schematically  illustrated  in  Figure  k,k  and  A. 5 
features  15  channel  analog  input  capability  and  offers  the 
option   of   selectable    ratss    of   2"samplB3    per    secoad,      where    n 
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may  take  on  integer  valuis  of  3  to  7.  By  appinpria'^ely 
jumpering  the  analog  input  pins,  ta?  sampling  ra-e  may  be 
increased  five  fold.  The  s/stea  is  ^arrently  drployed  wir:h 
the  coil-x  output  connected  t:^  PCM  inpa^  channels 
1,4,7,10,13  and  coil  y  Dutput  coi^iected  to  Pd  input 
channels  2,5,8, 11, 1U  with  the  reiilnSer  of  ths  channels 
grounded.  With  a  selectei  sampling  rate  of  32  i-amoles  per 
second,  this  configuration  results  in  153  samples  pf.-:  second 
of  Bach  coil's  signal.  While  chan.'irl  input  ciPhnitY  and 
sampling  rate  capability  ars  much  greater  rhan  cur-C-nt  needs 
dictate,  they  provide  growth  potaitial.  Th.^  aciitional 
channels  could  be  utilize!  for  a  thrB=  coil  syst-sa,  input 
froa  a  pressure  sensor  fDr  swell  wave  da-a,  temperature 
sensor,  or  real  time  clock  data.  rh=  BCyi  systex.  design 
includes  the  option  for  analog  or  iigital  data  output  from 
channel  14  with  a  minor  change  in  the  exterr.al  pin 
interconnections. 

The  ?C!1  systeoi  incccporatsr  a  cry?it-ii  oscillator 
and  associated  CMOS  integrated  circuitry  tc  lev-ilop  the 
clocking  pulses,  a  16  chanial  C13S  analog  multiplexer,  a  16 
channel,  12  bit  CMOS  a nalog-ro-1 igita i  converter  and 
associa-ced  circuitry  to  provide  tJi=  pulse  zz-lz.c^q.  The 
crystal  clock  oscillator  ooarating  at  a  frequency  cf  24.575 
KHz  produces  a  sguare  wave  output  witi  a  loss  bi-  rate  of  1 
bit  in  one  million.  This  freguenoy  combined  with  rhe 
inherent  delays  of  -he  issocLatei  clock-logic  circuitry 
provides      the    128      samples   per      seconi       sampling    rate.  By 

ex-carnal  pin  interconnections,  the  othsr  sampling  rates  may 
be  selected.  The  clock  pulses  gate  tha  analog  aalt iplexer , 
digital-to-analog  convertsr  and  follow  on  associated 
circuitry  that  form  the  pulse  code  words.  The  basic  output 
is  a  bi-phased  pulse  coied  signal  that  is  adjusted  to  TTL 
standards    by    adjustments    zo    the   offset    and   gain    iirimmers,    P2 
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and  P3  respectively.  The  3ata  is  organized  in  frames.  Z\ch. 
frame  is  headed  by  a  syn^  code  wDrd  which  is  fcilowe'l 
seqaentially  by  the  pulss  coded  sanpLes  from  PCI  charinels  1 
through  15.  The  sync  coda  word  is  a  pjlse  coded  iigital  word 
with  a  decimal  value  betwaen  0  ani  4096,  This  word  ±3 
preselected  and  hardwired  oi  the  circiit  board  by  connecting 
a  logic  high  or  lew,  as  appropriate,  to  pias  a  through  i  and 
k  tlirough  m.  rhs  sync  coie  »ford  is  assaatial  to  ths  decoding 
process.  Additionally,  it  providBs  protection  again-=;^ 
confusing  the  land  acquirad  data  wibh  that  taksn  from.  t'i5 
ocean  floor  as  each  sysiaa  has  it=  own  sync  cod=  wcrl. 
Currently  sync  code  words  35  38  and  2S20  are  associaiei  wi-h 
the  land  and  sea  sites  rsspactively.  Sync  code  word  3155  he.d 
been  used  in  test  and  saa  floor  data  collectioi  before  tie 
September  runs.  The  3  155  board  also  had  plus  or  minus  5  volt 
range  on  the  A/D  converter  whils  all  other  boards  whsT'^  plas 
or  ainus  10  volts. 

Gould  Gelyte  cechargabla  batteries  (plus  and 
minus  12  Volt  and  grounii  are  utilized  for  all  pow=r 
requirements  within  the  inst  rumentatiDn  sphere.  The  cur.v.ent 
consumption  of  the  preconi  itioner  ail  the  PC^l  system  is 
approximately  28  ma.  at  the  32  sample  per  second  rat~. 

The  externally  actuated  magnetic  ac-iv5.ricn 
circuit  described  in  section  II.A.I.a.  (3)  is  utilized  to 
activate  the  instrumentation  sphere  equipment  upon 
deployment. 

The  equipment  i nterconaeot ions  are  all  coaxisl 
except  that  the  power  inpat  is  laie  with  a  keyed  thr^e 
pronged  plug  and  socket  coanector.  -are  must  be  taken  to 
ensure  the  x  and  y  coil  output  leads  are  connected  to  their 
respective  PCM  inputs. 

Appendix  D  discusses  the  details  of  data  flow 
in  the  PCM  system  utilized  in  the  generation  II  geomagnetic 
data  collection  system  design. 
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c.   Optical  Data  Link 

The  optical  lata  liak  sliainates  ths  concarn  of 
antsnna  effect,  discussed  briefly  in  section  I  while  intro- 
ducing a  variable  length  section  in  the  system  zhat  gives 
flexibility  for  operations  at  various  depths.  Fiber  op'ic 
system  design  con sideratiois  as  well  as  an  overview  on  the 
theory  and  state  of  fiber  optic  techiology  are  discussed  in 
Appendix  E. 


m  INPUT 
TRilNSMinER  PHASING 

•  0 

•  in- 


Pnrtr  Adjust 


DIGITAL  Vr 

COMMON 


TTL  OUTPUT 
ANALOG  OUTPUT 
ANALOG  COMMON 
DIGITAL  COMMON 


THRESHOLD  ADJUST 


Figure  2.6 


Simplified  Black  Diagraa  of  the  Optical  Data 

Lin< 


The  fiber  optic  data  link  consists  of  an  optical 
transmitter  and  receiver  aach  house!  in  an  environmental 
enclosure  and  connected  by  a  fiber  optic  cable  (see  Figure 
2.5|.  During  deployment  the  cable  is  stored  on  and  dispensed 
from  a  plas-ic  reel  fashioned  froa  i  20  inch  bicycle  rim 
(sea    figure   2.7    for    detailsi  . 
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(1)  0  ptical  rrin smitter/Rscsiver  Set  rh^  optical 
transmitter  and  receiver  ara  the  Burr-Bicwn  37131  and  3  713R 
respectively.  The    371 3 T      and   3713R      when      ccnasccsd    by      a 

suitable  fiber  optic  cable  form  a  25D  FCbaud  NRZ  fiber  cptic 
data  link  capable  of  operation  to  1 .  7  kilometer.  The  371 3T 
fibsr  optic  transmitter  is  an  elect:  ical-to-cpt?lcal  -rans- 
ducer  designed  for  digital  transmission  over  single  fiber 
channels.  Transmitter  circuitry  converts  TTL  levBi  inputs,  to 
op-^ical   pulses   at   data   ratss    frooi    D3    to    2    Kbaud   NRZ. 

The  3713P  fiber  optic  rscsiver  is  an  opxic3l-to- 
electrical  transducer  designed  for  r5-9ption  of  iigital  c.a-a 
over  single  fiber  channels.  The  racsivar  circuitry  converts 
optical  pulses  to  TTL  laval  outputs  with  a  receiver  s-^nsi- 
tivity   of    15    nW    and    data   rates    to    250    Kbaud    NRZ. 

An  integratad  optical  connector  on  both  the 
3713T  and  3713R  allows  easy  intarfacing  between  aodijl^^s  and 
optical  fiber  without  problems  of  soarce/fiber/derir^or 
alignment.  The  me*al  pacicagas  of  tha  3713T  and  3713R  provide 
immunity  to  alactromagnatic  radiation  and  djirac*  prir.ted 
circuit  board  mounting  with  no  additional  heat  ^.-ink 
regaired. 

A  simplified  block  diagram  of  the  37  132  trans- 
mitter  is  shown  in  Figura  A. 5.  Tie  input  stage  useo  a 
Schimitt  Trigger  Exclusive  DR  Gate  5,  for  noise  ianianity  and 
its  logic  is  configured  so  the  phasing  of  the  -ransmit-.er:  is 
pin  programmable.  When  tha  transmittal  phasing  tarmin?.!  is 
connected  to  theta-0,  the  light  outpit  is  in  phase  with  tha 
digital  input  signal.  Tha  LED  is  on  whan  the  TIL  input  is 
high.  Connecting  the  transmitter  phasing  terminal  to 
theta-180  causes  the  reversa  no  happa^:  the  LED  is  on  for  a 
digital  low.  This  option  is  salectabla  on  the  ocean  floor 
geoaagnetic  data  collection  system  tirough  the  installation 
of   a    single    pole,    double    throw   switch.    The   system    is 
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operated  in  the  in-phass  moie.  The  advantage  cf  the 
theta-180  mode  is  it  makes  it  possible  to  detect  a.  breait  in 
the  fiber  cable  when  the  lata  link  is  idle.  P.Dwev-r,  an 
idle  PCM  signal  does  not  provide  2.  roDtinuous  TIL  low  and 
thereby  obviates  the  advantage.  i^mplifier  il  1  and  -^hs 
current  switch      drive  a      light   emitting      diode    (LSD).  The 

optical  power  output  of  the  371 3r  transmitter  may  be 
adjusted  by  controlling  tiie  resistance  between  ~he  power 
adjast  pin  and  ground.  This  contrDls  the  peak  or  '-'on" 
current  in  the  LED.  When  the  resis-3r  is  miniiam  the  LSD 
output  is  maximum  and  as  the  resistance  is  increased,  the 
LED  output  asyrapt orically  approaches  a  minimum  value.  In  the 
ocean  system  the  resistor  is  adjusted  to  provide  a  240 
millivolt    peak-to-peak  receiver   output    signal. 

A  simplified  block  diagran  of  the  3713R  rec=;iver 
is  shown  in  Figure  A. 7.  Input  light  is  converted  tc  a 
current  by  the  PIN  photodioie  CR^  rfhich  is  connected  in  zh'-^ 
photovoltaic  mode  for  maximum  sensitivity.  A  low  bias 
current  FZT  input  current- to-/ clrage  converter  transf  crras 
the      diode   current      in-o    a      voltage    (S^a)  which    is      further 

amplified  by  A2  and  presented  to  comparator  A3  as  Vb  whera 
it  is  compared  to  the  threshold  voltage  Vt.  For  iiiaxi-iiurR 
noise  immunity  iz  is  desirable  to  have  the  threshold  voltage 
set  to  a  value  corresponding  to  a  level  half  way  between  the 
higa  and  low  value  regardless  of  nae  actual  light  level  at 
the  input.  In  the  3713R,  this  is  accomplishei  by  a  peak 
detector  automatic  thresholi  circuit.  A  pulse  of  light  input 
causes  a  voltage  pulse  at  7b  which  is  stored  ia  the  auto- 
matic threshold  circuit,  iLvided  ia  half,  and  supplied  to 
the  comparator  as  the  thresholi  input  Vt.  Thus,  Vt  is  a 
voltage  corresponding  to  the  midpoint  of  the  light  and  no 
light  conditions  cf  the  dioie.  Since  the  automatic  threshold 
circuit   uses      a    capacitance      hold    technique,         the    threshold 
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voltage  Vt  is  subject  to  decay  whan  the  light  is  ren^cved 
frcM  CR1.  A  no-light  condition  of  approximately  0.5  sscoud 
duration  (a  1  baud  data  rata)  can  bs  jsad  with  no  significnt 
effact  on  noise  immunity.  The  PCM  system  presents  data  to 
the  optical  data  link  at  about  a  5K  baud  rate  aid  is  there- 
fors  within  the  conditions  for  maximized  noise  iamuniiiy.  The 
analog  output  terminal  of  the  3713R  is  the  output  of  ths 
linsar  amplifier  A1.  The  voltage  at  tiis  terminal  is  i-ropor- 
tional  ^o  the  input  power  to  the  receiver.  As  such,  i~  makes 
an  excellent  diagostic  point  for  testing  the  fiber  optic 
cab'.e.  Monitoring  the  analog  output  terminal  gives  a  rela- 
tive measure  of  cable  loss  at  the  transmitted  wavelength  and 
a  direct  measurement  of  receiver  sigi al-to-noise  ratio  when 
the  transmitter  is  off.  In  the  ocean  floor  geomagnetic  data 
collection  system  the  modilation  iaput  is  taicen  from  the 
optical  receiver  analog  output  sin-=  the  output  amplitude 
can  be  adjusted  by  varyiig  the  optical  transmitter  power 
adjust  resistor.  An  output  level  of  233  to  240  millivolts  is 
utilized  to  ensure  RF  traismitter  aodalation  without  over- 
driving  its    circuits. 

(2)  Ziki.1.  QBlLi  ^^ble  The  fiber  optic  cable  is 
described  in  detail  in  appendix  E.  The  fiber  is  a  63.5 
micrometer  radius,  low  n^ise,  giraied  index,  multimcde, 
single  strand  optical  fiber  cladded  and  sheathed  wi-h 
strength        members  as        illustrated         in  Figure        E. 2. 

Approximately  300  feet  of  type  5053  oable  is  stored  on  the 
storage  reel  shown  in  Figire  2. 7. a.  The  fiber  optic  link 
gives  flexibility  in  systen  operating  depths  as  all  other 
system  parameters  are  fixed.  The  cable  is  affixed  to  the 
polypropylene  line  with  electrical  tape  at  aboat  ten  foot 
intervals  with  some  slack  to  allow  for  stretcaing  of  the 
main  support  line.  The  eleotrical  tape  is  fastened  in  such  a 
way  as  to  provide  a  tad  that  aids  la  quick  renoval  during 
reco  very. 
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(3)  Coa^linjs  The  ^oupLiig  system  iiltilizsd  is 
the  AMP  Optimate  Single  Position  Fiber  Dptic  Cable  Connector 
System.    Since  AMP   connectors    were   not   available   wi-h 
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Figure   2.8        Fiber  Optic  Cable   CDnnector   Assembly 

parameters  compatible  with  the  fiber  Dotic  cable  intended 
for  the  system's  optical  iita  link,  the  fiber  optic  cable 
was  given  an  additional  ja-ket  of  lait  shrink  tubing.  The 
heat  shrink  tubing  builds  up  the  oatside  diaister  of  the 
cable  to  near  that  of  the  inside  di=aeter  of  tie  connector 
as  shown  in  Figure2.8.  Iha  optic  cable  is  allowed  to  extend 
aboit  1/16  inch  from  the  f=ce  of  the  -onnector.  The  cable 
and  the  heat  shrink  tubing  is  bonded  to  the  connector  shaft 
by   a    few    drops   of   cyanoazr/l ate      glue.       The    retaining    cap   is 
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secared  in  a  position  about  the  fsrrule  by  tha  crimp  rlrig 
mounted  1/16  inch  behind  ti2  cap.  The  polishing  bushing  is 
placed  ever  the  connector  and  securely  mated  to  the 
retaining  cap.  The  polishing  i=  accDuplished  by  applying  a 
figare  eight  moticn  to  the  fiber  against  various  grit  abra- 
sive paper  while  it  is  subasrged  under  a  half  inch  of  »/ater. 
A  9  X  9  inch  aluminum  baking  pan  is  very  suitable  for  the 
procedure.  Polishing  is  accomplishad  using  5  graders  of 
silicon  carbide  finishing  paper;  320,4  00,600,3  aicron,  0,3 
micron  in  that  order.  rh=  polishing  is  completed  when  the 
surface  of  the  polishing  bjshing  imaeiiately  surrounding  the 
ferrule  is  polished  to  i  glossy  finish.  The  polishing 
bushing   should   not  be   reiis  =  i. 

The  AMP  connectors  reaove  concern  of  optical 
power  losses  due  to  lateral  misalignmsnt ,  angular  misalign- 
ment,   end   separation    and   end    preparatiDn    quality. 

(4)  Envirqnmsi  tal  Encl2.=  iL§  The  optic  trans- 
mitter and  receiver  are  =ach  hcusei  in  an  snvircnmental 
enclosure  that  maintains  water-tight  integrity  in  the  ocean 
environment.  The  enclosure  illustrated  in  Figure  2.9  is 
manufactured  from  6  inch  I.D.  ? VC  ?i?e.  One  end  of  the  pipe 
and  its  associated  end  cap  are  machiiel  to  a  trae  round  as 
near  to  initial  dimensions  as  possible.  Typical  finished 
dimsnsions   are    : 

PVC    PIPE  1.3.        6.090     t    0.003    inches 

O.D.         6.600      +    0.003    inches 

PVC    END    CAP      I.D.         6.615      ±    0.015    inches 

An  "0"-ring  groove  is  machined  to  0.05  3  inches  isep  by  C.160 
inches  wide  located  0.75  inches  fron  the  machined  end  of  the 
pipe  to  conform  with  a  6.25  X  1/8  inci  "0"-ring.  The  pipe  is 
reinforced  with  rings  male  Df  0.75  inch  sheet  ?VZ  material. 
These  rings  are  slotted  3.25  inch  below  ring  center  to  serve 
as  a    support   for    the    optical    transmittar    or    receiver   circuit 
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boards.  The  machined  eni  cap  is  fitted  with  two  Brantner 
connectors  and  a  vent  plug.  Ths  opposite  end  cap  is  fitted 
with  the  fiber  q)tic  throagh-hull  psnetrator.  All  these 
fittings  are  equipped  witi  an  "0"-i:i2g  seal.  The  end  cap 
supporting  the  fiber  optic  through-hull  penetratDc  is  bonded 
to  the  unmachined  end  of  ths  pipe  witii  PVC  cement.  The  fiber 
optic  cable  is  passed  thrDigh  the  peiatrator  and  initially 
fastened  with  cy anoacrylats  glue.  A  final  treatment  of 
DEV30N  5-minute  epoxy  serves  to  secure  the  cable/penetrator 
interface. 

The  fiber  outside  the  penetra-'-or  is  protected 
froi  bending  stresses  by  application  Df  three  layers  of  heat 
shrink  tubing  that  extend  Lato  the  peaetrator.  Each  layer  of 
tubing  is  longer  in  length  than  the  preceding  Dne  so  that 
the  cable  flexibility  increases  as  a  function  of  distance 
frooi  the  penetrator.  Initial  design  tests  iniicated  the 
weaicest  point  in  the  optical  linx  was  the  point  at  which  the 
cable  entered  the  penetrator.  Tests  subsequent  to  the  appli- 
cation of  the  tubing  indicated  this  point  to  be  one  of  the 
strongest. 

(5)  Battery  P2ii£  SaoElies  The  optic  transmitter 
is  powered  by  a  6  volt,  8.5  AH  rechargeable  sealed  lead  acid 
battery  configured  with  a  100  ohm  secies  liraiter  resistor. 
This  battery  is  housed  in  a  3enth03  glass  sphere  mounted 
with  the  insx.rumentation  sphere  in  the  equipment  mounting 
assembly  discussed  in  11 . A. 1 .b  .  ( 1) . 

The  optic  receiver  pDwer  supply  is  a  Gould 
Gelyte  (+  12,  -  12  volt  ani  ground)  battery  pack  mounted  in 
the  can  buoy.  The  battery  pack  is  ratei  for  18  AH. 

d.   Can  Buoy 

The  can  buoy,  illustrated,  in  Figure  2.7, 
provides  71  pounds  of  positive  buoyancy  (this   is  inclusive 
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Figure  2.9   Optic  rraasmitter/Raceiver  Enclosure 
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of  battery  weight)  to  the  s/stem.  It  contains  a  Soald  Gelyte 
rechargable  (+  12,  -  12  volt  and  groiind)  battery  pack  that 
powers  the  optic  receiver.  The  power  is  delivered  through  a 
three  conductor  cable  equipped  with  Bcantner  connectors.  The 
can  buoy  is  fitted  with  ailitional  Brantner  CDnnectors  for 
data  transmission  in  and  oit  to  the  spar  buoy. 

The  axle  of  the  fiber  Dptic  storage  reel  is 
mounted  in  the  center  of  the  can  buoy  cover  containing  the 
electrical  connectors.  T.ia  opposite  cover  is  fitted  with 
quick  release  clamps  which  are  used  to  secure  the  buoy  to 
the  main  support  line. 

e.   Telemetry  Spar  3uoy 

(1)  Descrig^LDn  The  spar  buoy  designed  and 
constructed  by  Dr.  M.E. Thomas,  Applisd  Physics  Laboratory, 
Johns  Hopkins  University  and  LT .  P.H.  RutherfDrd,  Naval 
Postgraduate  School  [5],  is  depictai  in  Figure  2.10.  The 
spar  buoy  is  partitioned  into  several • sections  as  follows: 

(a)  Transmitter  Power  Supply/Lead  Ballast. 

Inclidiag  the  dampening  plate,  this  is 
the  bottom  two  feet  of  the  MPS  spar  buoy.  It  is  constructed 
of  5  inch  PVC  tubing.  This  section  coitains  three  6  volt/8.5 
AH  lead  acid  batteries,  which  proviie  power  to  the  trans- 
mitter, plus  ninety  pounis  of  lead  weights.  Water  tight 
integrity  is  maintained  by  a  3/U  in~h  thick  PV3  plug  with 
two  "0"-rings  to  aid  in  pressure  sealing  (Figure  2.11),  The 
power  leads  from  the  batteries  are  soldered  to  a  pin 
assembly  permanently  housed  in  the  plug.  The  transmitter 
receives  power  through  a  cable  whici  attaches  to  the  pin 
assembly  in  the  plug  by  way  of  a  Brantner  connector.  The 
slip  coupling  adjoining  the  bottom  and  free  flooding 
sections  is  permanently  adhered  to  tie  bottom  section  with 
PVC  cement. 
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(b)  Free-FlDoding   Saction 

The  purpose  of  ths  free  floDding  section 
is  primarily  to  enhance  the  separation  of  ths  center  of 
buoyancy  (cb)  and  the  canter  of  gravity  (eg).  Phis  section 
is  constructed  of  a  5  foot  length  of  6  inch  PV3  pipe.  A 
series  of  holes  are  drilled  -cowaris  the  bottom  of  the 
section  to  act  as  flood  ports.  Holas  at  the  top  are  vents. 
The  power  cable  from  tha  batteries  to  the  transmitter  tra- 
verses this  section. 

( c)  Middle  Section 

This  ssction  is  composed  of  three 
sub-units,  the  lowest  portion  hoasinj  the  transaittar.  The 
traasmitrer  utilized  is  taken  from  a  P3C  Orion  soncbuoy 
system.  It  has  a  power  output  of  oaa  watt  and  operates  in 
channel  29  at  a  frequency  of  171.65  HHz.  Its  watertight 
intagrity  is  maintained  by  a  PVC  plug  identical  to  rhat 
hetween  the  lower  and  fraa  flooding  sections.  The  power 
cable  leading  from  the  bottom  plug  is  attached  to  this  plug 
with  a  Brantner  connector.  Pracautionary  measures  that  are 
to  be  taken  during  this  portion  of  the  assembly  will  be 
discussed  in  the  next,  section  antitlad  transportation  and 
assambly.  Just  above  the  transmitter  housing  is  the  buoyancy 
collar  shown  in  Figure  2.12.  The  collar  providas  the  spar 
tuoy  with  an  additional  aigaty  pounds  of  positive  buoyancy, 
again  enhancing  further  saparation  of  the  cb  and  eg.  The 
remaining  portion  of  tha  niddla  section  is  composed  of  the 
middle  mast  which  houses  the  transmitting  cable  leading  to 
the  antenna.  This  middle  mast  is  rainforced  by  1/2  inch  and 
1  inch  PVC  tubing.  Tha  1/2  inch  sections  are  each 
individually  sealed.  Tha  1  inch  section  serves  as  a  channel 
for  the  transmitting  cable.  This  saction  is  saaled  using  a 
PVC  ring  and  adhesive.  Tha  antenna  cable  protrudes  through 
the  center  and  the  hole  is  nade  watar  tight  with  application 
of  silicon  based  rubber  sealant. 
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(d)  Upper  tlast/&nteana 

The  remiiniag  portion  of  the  telemetry 
spar  bucy  couiis-f  s  of  an  eight  foot  -ection  of  1/2  inch  PVC 
tubing.  It  houses  t.ie  remaining  lenjrh  of  the  transmission 
cable  and  is  capp'9d  by  the  transnit ting  antenna  which 
utilizes  a  reflecting  ground  plane  wiich  is  at  150  degrees 
from  the  vertical. 

(-■)  X£llis20rta  tici/Asseabl^  The  spar  buoy  is 
transferred  to  and  from  the  R/V  Acanii  via  a  piclcup  xruck  in 
most  ins-cances.  This  method  of  transportation  requires  that 
the  spar  buoy  be  easily  broken  down  into  lengths  rhat  are 
compatible  with  the  length  of  the  piocup  bed.  The  following 
subsections  diiicass  the  procedures  fDr  the  assembly  of  each 
affected  section  and  the  precautions  required. 

(a)  Power/Ballast  and  Free  Flooding  Sections 
Place  the  lead  ballast  in  the  bottom 
section  prior  to  inserting  the  batteries.  Ensure  that  the 
"0"-ring£  are  sufficiently  coated  wirh  a  thin  film  of  high 
vacjum  grease,  ynsert  plug  so  that  it  snuggly  fits  below  the 
marked  level.  CiDTION:  prior  to  coniecting  the  power  cable 
jack  to  the  plug,  ensure  that  the  indentation  in  the 
Brantner  pluq  is  linsd  up  with  its  counterpart  on  the  jack. 
Failure  to  do  so  could  iamige  the  transmitter  by  reversing 
the  polarity.  iJext  slip  the  free  flooding  section  into  the 
joint.  The  holes  that  allow  flooding  are  to  be  at  the 
bottom.  Ensure  that  the  narks  on  both  slip  joint  and  free 
flood  section  are  aligned.  Bolt  the  free  flood  section  in 
place, 

(b>  Free  Flooding  and  Middle  Section 

Prior  to  joining  these  two  sections, 
insert  the  transmitter  assemoly.  First,  connect  the 
transmission  cable  then  slide  the  traismitter  in  place  again 
ensuring  a  snug  fit  and  proper  placement  of  the  plug  in  line 


U3 


Hoisting 
Collar 


\. 


^ 


^ 


r 


\ 


1/4"  Al^imir.um 
Plates 


<Z5 


^ 


C3 
c3 


/ 


Polyethylene 
Collar 


/ 


\ 


Mounting 
Collar 


Scale:    1/4 "=1" 


.J 


Figure   2.  12        Telemetry    Spar   Baoy   Buoyancy   Collar 
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with  the  interior  markings.  CAUriON:  It  is  important  to  have 
a  Sit  ohm  dummy  load  on  the  transmitter  at  this  point  of 
asssmtly  since  it  is  beiag  powers!  without  the  antenna 
attached. 

( c)  Middle  Section  to  Opper  Mast/Antenna 

This  pDution  of  the  spar  buDy  is  to  be 
bolted  in  place  while  the  remaining  DDrtion  is  hanging  over 
the  side  of  the  ship  hoisted  up  by  the  Acania's  on  board 
crane.  Ensure  that  the  traasmission  cable  is  conaectsd  prior 
to  ipper  mast  attachment. 

(3)  Degloxmeat/Recover^/Dis assembly  All  sections 
of  th9  spar  buoy  shall  be  assembled  :)nboard  the  licania  with 
excsp'ion  of  the  upper  mast  priDr  to  iaployment.  Attach  ^he 
crane  line  approximately  five  f=9t  up  from  the  tathar  point 
on  the  spar  buoy.  Once  lo^ared  over  the  side,  the  buoy  can 
be  tempcrarily  secured  to  the  ship's  railing  and  lowered  an 
additional  amount  to  enabLa  attachmeit  of  the  upper  mast. 
The  buoy  is  next  lowered  into  the  water  and  tended  along 
sid5  until  the  rest  of  the  3yst.em  CDmponents  are  deployed. 

Recovery  and  disassembly  is  accDaplished  by 
following  the  reverse  of  tl=  above  procedure. 

f.   Receiving  Station  Sguipmeat 

The  receiving  station  agiiipment  consists  of  an 
AN/5.RR-52A  radio  receiving  set,  an  H?  396UA  instrumentation 
tape  recorder  and  two  dual  trare  oscilloscopes.  AN/PRC  77 
radio  sets  serve  as  a  voic=  communications  link.  The  equip- 
ment arrangement  is  as  shown  ia  Figure  2.13.  The  annennas 
used  in  ccnjuction  with  tha  telsmenry  link  are  double,  tuned 
YAGI  type. 
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Figure  2.13   Receiving  Station  Egaipaent  Configuration 
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III.   DAIA  PSOCESSINS 

Currently  the  data  is  processed  tD  develop  various  data 
plots  generated  on  the  NavaL  PDstgcaiiiate  School's  IBM  3033 
computer.  Prior  to  the  actual  gsiaration  of  the  plots 
however,  the  PCM  encoded  data  whica  was  store!  on  analog 
tape  utilizing  the  equipment  described  in  section  II.A.I.f. 
(receiving  station  eguipmeat)  must  be  decoded,  digitalized, 
and  transcribed  onto  digital  tape. 

A.   HARDWARE 

1 •   PCM  to  Digital  Conversion  Bjuigment 

Several  electronic  components  are  utilized  to  decode 
the  PCM  data  to  discrete  analog  fora  and  then  to  digital 
form  which  is  ultimatly  stored  on  9  tcacic  digital  tape.  This 
data  processing  system  is  siown  schsmat ically  in  Figure  3.1, 
illustrating  data  flow  and  control  links.  Central  control 
of  this  process  is  accomplished  with  a  Hewlitt  Packard 
9845A  computer  utilizing  in  operatic  interactive  program 
titled  "PCM  ?HOG".  After  execution  of  the  program,  the 
computer  requests  entry  of  specific  finction  control  parame- 
ters into  the  computer  and  other  equipment.  These  inputs  are 
used  to  control  synchronization  of  egiipment  start,  digital 
tape  drive  speed,  decode  rate,  decoda  time,  and  synchroniza- 
tion code  word  entry  into  trie  decodar.  The  PCM  encoded  data 
is  fed  into  the  system  from  the  HP  3964A  4  ohannel  tape 
recorder  previously  used  to  store  the  data.  Utilization  of  a 
different  recorder  may  introduce  noise  due  to  a  difference 
in  head  alignment  and  drive  speeds.  Such  noise  decreases  the 
sign al-to-noise   ratio   apoi   which   decode   synchronization 
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depends.  Should  decode  syar hronizatiDn  not  occur,  a  lower 
analog  tape  drive  speed  aai  simple  rate  must  be  selected, 
resulting  in  longer  processing  time.  The  decoding  of  the 
PCM  encoded  data  is  accomplished  with  a  Marine  Profiles, 
Inc.  Model  319  PCM  decoder.  Two  Monsanto  AM-64 1 9/USM-368 
oscilloscopes  are  used  to  iisplay  the  PCM  encoded  data  and 
the  discrete  analog  data.  The  discrete  analog  data  is 
processed  by  the  computer  tD  develop  the  digital  data  which 
is  recorded  on  6250  bpi/123  0  ft.  9  track  digital  storage 
tape.  A  Kennedy  Model  99D3  digital  recorder  and  computer 
interface   are  employed   for   that   purpose. 


D-discrete    analog 
C -  control 
DL-  dig  ital 


c  omputer 


T 
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-c-^ 


)'L     \       ( 


1 i    °\\ 

(interface    J     Ydecoder 
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Ke  nnedy 
J  ta  pe 
drive 


Figure   3.1        Decoding   System   Data    Flow   and   Control 
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a.      Decoding    Proceiare 

The  following  steps  illustrate  the  procedure  for 
decoding  the  PCM  encoded  raw  data  into  digitalizsd  data  and 
then  placing  it  onto  a  9  track  digital  tape  whioh  can  than 
be  delivered  to   the   IBM  3033    coraputar    for   data    processing. 

1.  Disconnect  the  coaxial  cable  labled  number  '^ 
from  the  outpun  of  the  HP  39  68A  rapa  recorder  located  in  ths 
equipment   rack   adjacent  to   the   H ?   9345A    computer. 

2.  Connect  cable  7  iato  the  output  channel 
desired  for  decoding  of  the  HP  3954A  U  channel  tape 
recorder. 

3.  Energize  the  HP  396'4A  tape  recorlar,  the  two 
Monsanto  AM-6419/USM-368  osoillo sopes,  the  Kennedy  interface 
and  AANDEEAA  tape  transfer  interfaoB,  the  moi  =  l  319  PC:i 
decoder,  and  the  HP  gS^+SA  conputsr  (the  ON/OFF  swi-ch  is 
located  on  the  right  hand  side).  Also,  on  the  9Z'A  model  3  19 
decoder  place  the  fan  toggle  switch  in  the  UP  ?osi-ion  and 
the    AC/DC/OFF   toggle    switch    -o    the    A3    position. 

4.  Place  into  the  right  hiad  side  taps  reader  of 
the    HP    98a5A   computer    ths    program  nainsd    "PCM    PR33". 

5.  Mount  ths  analog  tape  onto  the  HP  39dUA  tape 
recorder. 

6.  On  the  HP  93U5A  conipjter,  type  the  command 
GET   "MT"    and   press  EXECUTE    to   initiata    processing. 

7.  On  the  PCM  decoler  plaoe  the  following  func- 
tions   to    the   listed    positions: 

SOU  5CE    -    1 

SAMPLE  HATE  -  64  (for  3  3/U  recorder  speed) 

-  128  (for  7  1/2  recorder  speed) 
INVERT/NORMAL  -  NORMAL 
OUTPUT/SAMPLE  RATE  -  3 
RECORDS/FILE  -  INFINII? 
SYNC  CODE  -  930 
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8.  Press  RUN  oi  the  aompater,  ignore  "enter  y  to 
skip  tape  init"  (displayed  on  rhe  computer's  CRTi  ,  and  press 
CONriNUE. 

9.  The  compatar  now  indicates  "loai  tape"  into 
the  Kennedy  unit  and  "put  dq  line".  Td  do  this  energize  the 
Kennedy  ur-it,  load  the  zape  according  tD  the  diagram  located 
on  the  inside  of  xhe  unit's  door,  pc=33  the  LOAD  button  and 
the   0N-LIN3    button   located    on   the    froQt    on    the    Kennedy    unit. 

10.  The  CDmpirer  now  indicates  "enter  sync 
code".  Type  into  the  computer  3658  dc  3155  as  appropriate 
for  the  sea  svstsm  or  2623  for  the  land  site  system  and 
press    CONTINUE, 

11.  Enter  transfer  tima  in  minutes  and  seconds 
into  the  computer.  Example  20  minutes  and  20  seconds  would 
be  typed   in    as    20,20.    After    this    is    d^ns    press    CONTINUE. 

12.  Push  the  STOP  switci  on  the  PCM  unit  and 
push    CONTINUE   on    the    computer,. 

13.  Push  the  START  switch  on  the  PCM  unit  simul- 
tansously  with  the  PLAY  button  on  th=  Hewlett  Packard  tape 
recorder. 

14.  If  data  transfer  must  be  ended  early,  push 
the  KO  button  on  the  computer  to  halt  data  transfer.  If 
this  option  is  elected  taen  I  Bust  be  entered  on  the 
computer  following  the  KO  command  to  write  end  of  tape  to 
the   tape. 

15.  "End  of  run"  will  be  indicated  on  the 
computer  CRT.  At  this  time  deenergizB  all  the  equipment  and 
disconnect  cable  7  from  the  4  chan.i5l  tape  ricorder  and 
reconnect  it  to  the  8  chaaiel  tape  recorder  located  in  the 
equipment    rack. 

16.  Finally,  remove  th=  yellow  ring  from  the 
digital  tape  to  ensure  that  the  tape  is  not  inadvertently 
over    written. 
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2.       Maill  Q23LEi^§I.  Mrdirare 

The  data  is  read  from  the  digital  tape  into  the  IBM 
3033  computer  via  software  ievelopai  by  Lcdr.  J.  Fisher  of 
the  Naval  Postgraduate  SchDol.  The  data  is  unpacked  (placed 
into  a  format  where  it  can  easily  be  accessed  and  manipu- 
lated) and  converted  into  a  spectrum  of  frequencies  using  a 
Fast  Fourier  Transform  routine.  This  data  is  then  plotted 
usiag    a   peripheral  Versatec    plotter   f^r    spectral    display. 

B.       SOFTWARE 

The    programs    are      written    in    the    FORTRAN      IV    programming 
language   and  are    briefly   discussed   below. 

1  •      I§.§1   Prog  ram 

The  test  program  si  own  in  Appendix  F  was  developed 
to  test  the  computer  manipulation  against  known  data  to 
ensure  the  resulting  plots  wouli  m=rt  expectations.  With  the 
exc2ption  of  the  internally  generatsi  lata,  tha  program  is 
essentially  the  same  as  tha  Main  Program  shown  ii  Appendix  H 
and  explained  in  paragraph  3  below. 

2.   Mass  Storage  Program 

The  digital  data  tapes  are  mLabled  with  data  organ- 
ized in  16  word  frames,  tha  first  word  of  which  is  the  sync 
code  followed  by  15  data  words.  The  lata  words  are  organized 

x,y,z,x,y,z, Since  the  sea  sys-em  has  only  two  sensor 

coils  at  this  writting,  ths  z  data  is  zsro.  The  land  system 
data  has  actual  values  for  )C,y,  and  z.  The  data  from  both  a 
land  and  sea  data  tape  is  altered  intD  the  IBM  3033  computer 
mass  storage  system  (MSS)  utilizing  tha  mass  storage  program 
shown  in  Appendix  3.  Curreatly  only  tia  x  and  y  coil  data  is 
read  into  MSS  where  it  is  available  for  access  by  the  Main 
Program. 
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3-      Main  E£22^31 

The  Main  Prograi  atilizes  th2  x  and  y  ceil  inputs 
from  simultaneous  land  and  sea  data  to  develop  various  plots 
for  comparision  and  analysis.  The  Main  Prograa  utilizes  a 
large  number  of  arrays  with  the  intention  of  naintaining 
each  bit  of  original  data  or  cilculation,  ready  for  recall 
at  any  time  for  futher  coaputat ion.  Equivalent  arrays  are 
employed  because  the  plotting  routine  is  not  able  to  handle 
complex  arrays.  Ocean  arra/s  are  iniicated  with  the  suffix 
'0'  and  land  arrays  with  the  suffix  'L*.  By  successive  calls 
of  "SOBBOUTINE  RD ",  an  array  of  the  x  data  and  y  data  called 
xxx(I)  and  yyy (I)  respectively  is  established.  One  such 
array  is  called  a  block  of  lata.  The  data  words  in  the  array 
are   integer    values  between   2 0U8    and    4D96. 

After  a  block  of  data  is  entered,  the  integer  values 
are  converted  to  complex  ^^alues  of  ^r^ltage  as  seen  at  -he 
front  of  the  PCM  board  between  -13. D  and  +10.3  volts  for 
systems  utilizing  PCM  boards  i/ith  3fnc  codes  Df  2620  and 
2653.  For  systems  utilizing  sync  coie  3155  the  voltages  vary 
between  -5.0  and  +5.0  volts.  The  xxx  and  yyy  arrays  repre- 
sent data  tak'en  from  the  x  and  the  y  coils  respectively 
without    application    of   the    system    traisfer    functions. 

The  voltage  values  (which  van/  over  time)  are  trans- 
formed into  the  frequency  lomain  by  utilizing  the  subroutine 
"FOUET"  which  performs  a  Fast  Fourier  Transfori  (FFT) .  It 
utilizes    the   Cooley-Tukey    FFT    algoritim. 

This  section  of  the  prograa  applies  the  system 
transfer  functions  to  the  data  that  has  been  transformed 
into  the  frequency  domain.  The  discussion  of  the  origin  of 
the  transfer  functions  is  given  in  Appendix  C.  The  actual 
transfer  functions  vary  as  a  function  of  frequency.  They  can 
however,  be  approximated  by  straight  line  segments  for 
different    frequency   ranges.    These    line    segments   are 


52 


qen^rdtsd  by  least  squares  a pprD ximation.  The  data  entering 
this  secxioii  enters  with  units  of  voltage  and  is  converted 
into    units   r)f   nano-Teslas    b^    the    traisfer    functions. 

The  previously  mentioned  sertLoQS  of  the  program  are 
embsdCGd  within  a  do-ioop  which  allows  the  analysis  of  a 
long  period  of  data  consisting  of  a  number  (NR)  of  blocks. 
The  Ir  h  el<2m6nt  of  each  block  is  tian  multiplied  by  its 
complex      conjugate,  averaged      and        converted      to        power 

(nap. O'-Tesla  squared)  and  stored.  Dnce  the  program  has 
passed  -hrough  the  averaging  loop  foe  the  last  time  it  has 
computed  the  arithmetic  average  for  aach  frequency  point  in 
the   arrays. 

The  next  section  of  the  program  calcalates  the 
Stcices  paraneters  0  through  3  of  taa  individual  site  ortho- 
ganal  ccmpoients  and  the  coherence  Df  the  planar  circular 
pclarizatio  :i      parameters      between        sLt=      componints.  The 

remai;ia-5r  o:  the  program  sets  ip  the  parameters  and  titles 
for    the    plo:s. 

A  more  detailed  description  of  these  programs  may  be 
found    in   Reference  [6]. 
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::V.       EXPERIMENrAL    RSSQLrS 

A.  DATA    QUaLITY 

Test  data  points  were  ievsloped  by  applying  a  5-0  volt 
peak -to- peak:  rinuHoidal  inpat  of  1,5,10,   and  15   Hz  to  the 

PCM  board.  Trie  data  was  saapled,  P-M  encoded,  transmitted, 
recaivedp  t'U-'^d,  ciecoded,  iigitalizei,  and  transcribed.  A 
computer  assisted  voltage  vs.  time  plot  of  tha  1.0  Kz  data 
resulted  ii'  ?.  1.0  Hz  sinusoidal  sigial.  A  plot  of  the  trans- 
fori  of  the  5.0  vclt,  5.0  Hz  signal  resulted  in  an  impulse 
function  at  5.0  Hz.  Thase  results  illustrate  that,  the 
system  funcT.i'jns  as  designed  from  th2  PCM  input  through  the 
data  proce.?Gir:g  w.:th  the  exception  of  the  application  of  the 
system  "crr-nsfer  'functions.  The  system  transfer  functions 
(which  take  ir.to  .iccount  tha  amplification  of  the  preampli- 
fier and  sior.El  <:ondition9r  s  as  well  as  the  transformation 
from  vol"J3  -':o   nano-Teslas)  are  described  in  Appendix  C, 

B.  DATA    QC<A.V-,'ITY 

Tha  follc'/ing  data  acqaisition  ?  =  ricd  resulted  in  the 
indicated  si  jiulta  ;ieous  laid/sea  data  from  La  :iesa  Village 
and  Monterey   Bay: 

1.  '0:3  2    AM    to    2:'i3    ?M,     25    Jily    1982 

2.  11:22    AM,     17    Aug.     1982    to    09:10    AM,    13    Aug.    1982 

3.  09:32  k'A ,  16  Sbd.  1932  to  09:58  AM,  17  Sep.  1982 
This  represents  50  hours  aid  30  minutes  of  accumulated  data 
during   three   collection  runs. 

Sections  A  and  B  abov2  illustrate  achievement  of  the 
stated  goal  to  develop  a  data  acquisition  system  with  a 
capability  to  gather  simultaneous  land/sea  data  with  high 
quality    in   great    quantity. 
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C.       PLOTS    OF    PRDCSSSED    AZTJAL    DMA 

Appendix  I  cor.:~i.sts  3f  a  number  of  actual  data  plots.  It 
is  lot  the  int^-nt  here  to  laalyze  taase  plots  but,  to  merely 
illustra-ce  the  progress  of  the  Naval  Posrgraduate  School's 
Geomagnetic   Recrearch      Teiin   sffort.  Dn-going    work      in    data 

analysis   includes    power    5p2:;tral    dansity,    coherance,    correl- 
lation    and   direct  i:5r -^litf    studiss. 
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V-       C^WCLUSIOtTS     AND    RECOM?1ENDATIONS 

A.  CONCLUSIONS 

The  Generation  II  sy^'hn  d<?sign93  for  this  rasearch  was 
successfully  emplcysd  in  the  contimiing  effort  of  geomag- 
netic data  collection  with  th^;  specific  task  of  recording 
quality  simultaneous  lan'j/sea  lata  with  increased  quantity 
over  its  predecessor  systeas.  The  data  collectsd  wich  this 
system  was  only  limited  by  power  supply  capacity  and  by  the 
schedule  cf  the  Research  Vssse!.  Acaaia. 

B.  RECOMMENDATIONS 

A  better  means  of  po  si:.i  oning  the  sensors  on  the 
seafloor  still  remains  a  pcDbl'^m.  rh=  ability  to  accurately 
position  and  detect  sensor  or-sntation  will  becDcne  critical 
in  future  data  acquisitioi  and  compDient  field  ::orrelation 
experiments.  The  exact  o::ien:ation  nathod  to  use  will  be 
diractly  related  to  the  ac'iiisition  system  chosea  for  future 
measurements.  A  flux  5^.-=  ma  gnetoa  star  may  possibly  be 
employed  to  indicate  crian-ation  of  the  antenna  coil  sensor 
assembly   to   the   earth's   [[la^neti^    fiali. 

Although,  the  rescar^ti  team  aid  the  ship's  crew  were 
well  versed  in  the  syste-j*  s  configuration  and  deployment 
techniques,  the  deployment  time  averagad  about  '-^5  minutes  as 
a  rasult  of  the  system's  oize  and  complaxity.  Furthermore, 
deployment  could  only  be  iccomplishad  in  calm  seas  at  a 
maximum  depth  of  200  feat.  To  alleviate  these  problems  the 
design  of  a  generation  III  systam  is  racommended.  The  design 
should  include  size  reduction  aad  systam  simplification.  The 
simplification  may  be  achiaved  by  zoncantration  of  several 
genaration   III    components    iito   i    siagla    package. 


56 


Additional  problGms  to  ov^rcoms  in  future  systs-ms  are 
the  lack  of  robustness  of  the  ;ipar  buoy  and  fiber:  optic 
systems.  The  antenna  groar^d  plane  elements  where  manufac- 
tured from  aluminum  and  Owissionally  fractured  under  stress 
during  transportation  or  daployment.  Future  sysuems  should 
utilize  stainless  steel  -la.-aenrs,  rhs  spar  buoy  sank  on  one 
occasion  due  to  compression  of  zha  Duoyancy  collar  during 
heavy  seas.  The  buoyancy  collar  shDald  be  replaced  by  a 
noncompressable  floatation  drvica.  During  recovery,  the 
fiber  optic  cable  was  crushed  oa  one  occasion  and  kinked  on 
another.  This  problem  may  be  corrects^  by  utilizing  a  fiber 
opric    cable   that    is    jackets!    with   a   haavy   rubber    sheath. 

A  final  recommendation  is  to  Incorporate  a  pressure 
sensor  to  measure  surfaca  ind  mternil  wave  activity.  This 
additional  information  will  allow  th=  analyst  to  determine 
the   contribution    by    waves    to    ths    leou i gn etic    dati. 
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APP3NDEX    A 


A.       EQUIPMENT    SCHEMATICS 

The  equipment  utilizai  in  this  project  i/as  locally 
designed  and/or  constru:;tai  ,  comnar^ially  available  or 
available  through  military  chaanels  and  locally  modified. 
The  system  equipment  schamatics  ira  presentad  in  the 
following    figures: 

1.  Preamplifier    (Figura    A.  1| 

2.  Magnetic    Activation    Circuit    (Figure   A. 2) 

3.  Pre-amp   and   PCM   Prero nditionar    (Figure   A. 3) 

4.  Pulse   Code   Modulation    System    (Figure    A. 4,    A. 5) 

5.  Optical    rransmittsr    (Figare    A..S) 

6.  Optical   Receiver    (Figure   A.7j 
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Figure  i,1        Praaaplifisr   Schematic 
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Figure    a.2        Magastic   Activition   Circuit 
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Figure  A. 3   Pre-aap-  and  PCM  PreconditisEer 


61 


4'{i\AA44'^4'iN4 
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Figure    A. 5        Pulse   Code   lodulatioa    System    (continued) 
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APPENDIX    B 


A.       SYSTEH    MID    EQUIPMENT    PREPARATION 

In  order  for  equipment  to  psrform  and  operata  as 
designed,  certain  pvsventi/a  maint3c=nc=  rasasurss  and  proce- 
dural steps  should  re  xaken.  Such  prudant  sxeps,  especially 
in  undersea  research,  will  most  C5rtainly  c.llow  maxiaium 
results   and   prclona    equ2pni=nt   life. 

The  fcllcwing  lis-  serves  is  a  ciliinatio.i  oi  lessons- 
learned  in    previous   and   on   going    projrct    work: 

1.  Spherical  gla£;s  iistrunent  spheres  are  composed  of 
two  Benthos  glasi;;  heaispheres  which  hav=  precision 
ground  mat:.nq  surracBs.  As  *:he5e  surfaces  are  susep- 
table  to  chipping,  cire  must  ba  taken  in  handling. 
During  assertibi  y  oi:  tae  housi.ig  prior  to  deployment, 
cer-ain  oioc*id':res  ars  necessary  to  maintain  watertight 
integrity : 

a.  The  g„ass  spa  ere  should  be  dry  and  warmed  to 
room  ts  ffipt^rature  using  a  tisat  gun  to  ensure  adhe- 
sion   of.   the   sealing   tape, 

b.  All  Kating  surfaces  lust  be  cleaned  with  a 
solvent  such  as  alcohol  to  remove  lay  foreign 
niatter,    grease,    dc    oil. 

c.  Prior  to  sealing  the  sphsrs,  careful  inspection 
of  all  equipmenc,  electrical  connections,  and 
packing  material  should  b=  nade.  The  use  of  check 
off  sheets  and  cDlor  coding  i/as  found  to  be  parti- 
cularly     helpful   in  avoiding    operator    srror. 
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2.  Slirctrcnic  'equiprtsnt ••  A  gansral  concern  is  -he 
proper  applicat.ior.  of  power  supply  voltages.  Assembly 
of  a  large  number  of.  alectrical  compoLents  within  the 
confines  of  a  small  glass  sphere,  necessitates  careful 
planning  and  utiliriation  of  United  power  resources. 
Improperly  applied  voltacjss  nay  result  in  damage  to 
electronic  circuitry  and/or  1d=3  of  data.  Certain 
pieces  of  electrc-nir  equipment  requirs  particular 
consi'leratio  n  as    follDws. 

a.  Sensor  sysr.eni-  Afier  the  sensor  coils,  preamp- 
lifiers, and  their  battery  power  supplies  are 
sealed  in  their  respacxiv^  spheres,  a  voltage 
measurement  shDUj.d  ae  talcei  to  ensure  operation 
and  then  the  dea:^  xiva  tion  nignet  positioned  over 
the  reed  switch  sue:!  that  the  voltage  drops  to 
zero. 

b.  PCM  systea-  Ex  :i:  a  cars  should  be  taken  to 
ensure      the    x-roil      is      connected      to    the      x-input 

pre ro nditiD asr/amplif ier    and      its 


x-catput   is 


.iBctBCL    to      ths    x-input      channel    of 


the    PCM      systen.       Th a   same    holds      for    the    y      and    z 
channels.         An      srror   in 


result 


tha   *rong    transfer   function   being 


applied  to  the  iita  iuring  tha  computsr  assisted 
data  analysis.  Tie  procedaces  identified  for  the 
sensor  system  aoove,  hold  as  well  for  the  PCM 
system . 

c.  Tape  recorder-  The  racDrder  heads  and  capstan 
drive  should  be  :^leaaed  aid  demagnetized  before 
each  acquisition  run.  This  action  will  anhance  the 
signal  to  noise  ratio  and  maximize  the  decoding 
speed.  High  quality  1300  foDt  tapes  are 

recommended. 
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d.  Fiber  optic  lata  liak-  Before  packaging  the 
fiber  optic  da^i  lirA  enc].Dsure,  a  system  check 
should  be  run  t:>  datsriiria  that  the  output  is 
200-2'40mv.  Before  connecting  the  power  cables  care 
should  be  taken  to  ensure;  the  transniirer  power 
lead  has  the  gcDund  on  th.B  pin  marked  with  a 
dimple  en  the  connector  cashing.  The  diaple  should 
be  aligned   with  the   indent  on   the  female   power 


connecter  mounted 


i:he   5nd  cap 


the 


trans- 


mitter enclosure.  Similar  steps  should  be  takei 
regarding  the  PC:i  data  in:)ut.  This  concern  isn't 
necessary  with  tri2  receiver  pDwer  input  as  it  is  a 
keyed  connector.  One-  th-.'  connections  are  made 
remove  the  activating  magast  from  the  PCM  sphere 
and  observe  the  idle  PCM  s:.gnal  on  an  Dscillosope. 
Adjust  the  transit  j.  tt  er  pD.'?r  adjust  potenriometer 
so  thaT  200-2^3  ip,v  i.>  observed.  Set  the 
O-degree/1 80-deg:-a  e  phrase  switch  to  the  0-degree 
position.  Package  thf;  '^ncl:)3ures  as  follows: 

1)  Clean  "D  "-rings,  "0"-ring  grooves  and 
making  surfa::es  v;ith  =ilcDhol. 

2)  Lightly  lubricate  "C'-rings  and  pipe 
mating  surface  between  "3"-ring  and  pips  end 
wi-h  vacuuji  grease. 

3)  Open  veit  nl'igs  ani  gently  slide  end  caps 
onto  the  enclosure  ensuring  th^  electrical 
wires  are  clear.  CAUIIDN:  care  should  be 
taken  not  to  kink  tha  fiber  optic  cable. 

4)  Both  sncLosures  should  be  secured  to  the 
dispensing  r^el  as  shown  in  Figure  B.2  and 
transported  as  a  separite  compon=at. 
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3.  Baxteries-  Gelycs  battsiiss  ire  ussd  foe  most  power 
requirements.  Thc-se  batterias  iihould  be  checked  prior 
to  system  usage  and  a  systematic  approacii  regarding 
service  use  and  charging  interval  should  be  taken.  The 
mercury  batteries  usei  in  rhe  sBnsor  spheres  provide  a 
constant  power  soarc=  througtioii-:  their  recommended 
life,  but  should  be  rsplacsd  at  lore  frequent  intervals 
to   ensure   successful   systes    op?rit.on. 

U.  Connector/Cables  -  Befors  a.T.i  af-er  each  use  in  a 
corrosive      environment,  all      'jo.inectors        should      be 

thoroughly  cleaned  to  axtend  soirvice  life.  All  rubber 
connectors  should  be  lightly  lubricated  with  silicon 
grease  to  ensure  2l5Ctrical  c^ri.inuity  and  prevent 
water  intrusion.  The  use  of  zdIo":  coding  3a  cable  and 
connectors  lessen  thi  possibility  of  error  and  is 
highly   recommended, 

3.       SYSTEM    DEPLOYMSNT/RECCVSRY 

On  the  scheduled  day  :>f  system  Qsployment,  the  previ- 
ously prepared  system  components  ars  taken  to  the  Monterey 
Bay  Coast  Guard  Pier  by  pick  up  -ru::k  where  they  are  on 
loaied  to  the  2/7  Acania.  The  on  load  is  accomplished  with 
assistance  of  the  ship's  crew  ana  its  on-b3ard  crane. 
Figure  B.I  shows  the  shipboard  lay^ac  of  system  components 
that  has  proven  to  be  ths  uDst  5ffici5.;t  starting  point.  Not 
shown  in  the  figure  is  the  fiber  optir  data  link  and  storage 
reel.  It  is  temporarily  stored  in  tie  aft  deck  house  with 
its   interconnecting    electrical   cables. 

The  following  steps  are  taken  to  facilitate  final  prepa- 
ration and  deployment.  The  numbers  in  parenthesis  refer  to 
the   numbered   items  in    Figure    3.1. 
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step  1:  The  crew  passes  tha  main  s^uppor':  lina  (1) 
through  the  A-frame  boon  block  to  chs  wi.ricn  on  top  of  the 
aft  deck  house.  The  line  is  reelai  en  -o  t.ha  winch  drum 
allowing  enough  slack  so  ths  line  ina/  easily  be  run  down  the 
port  side  of  the  ship  out  of  the  way  ^f  th^  follDi/  on  opera- 
tions. The  other  end  is  secured  to  the  sensior  system  (2) 
harness    by   means    of   a    bowline. 

Step  2:  The  spar  buoy  is  a3S2iDl=d  startinj  with  the 
power  supply  section  (6)  aid  the  frei  floodiiig  section  (7). 
The  interconnecting  cable  (11)  is  ooanscteC  to  ths  recep- 
tacle in  tha  power  supply  section's  water  "igh-r  plug  and 
passed  through  the  free  flooding  section.  RenovB  six  1/4-20 
scraws  from  the  freeflooiiag  section,  \lign  the  mating  index 
marks  on  each  section  ani  mate  the  sections.  Additional 
alignment  may  be  achieved  Dy  applying  a  long  shifted  screw 
driver  as  a  lever  in  the  free  flooding  ports.  Secure  the 
sections  with  the  screrfs  previously  reaovel.  Before 
connecting  the  power  cable  to  the  trar.smitt -i^-j  section  (8), 
ensure  the  dummy  load  is  in  place  on  the:  sectio'i's  antenna 
output  coaxial  cable.  Make  the  power  oable  connec:ion.  Mate 
and  fasten  .the  free  flooding  and  t  r  ansmi  t -ler  s  actions  as 
before.  Do  not  assemble  the  antenna  uast  (d)  jV.  the  spar 
buoy    at   this    point. 

Step  3:  Take  the  fiber  optic  traa saitter,  csceiver  and 
storage  reel,  packaged  as  shown  in  Figur-3  5.2  from  the  aft 
deck  house.  In  the  vicinity  of  the  oan  buoy  (5|  remove  the 
electrical  tape  securing  the  traasnitter  to  tna  storage 
reel.  Using  care  not  to  break  or  kink  the  optic  oable,  reel 
off  enough  cable  so  that  the  transmitter  may  be  positioned 
on  the  instrumentation  seotion  (4| .  Carefully  align  and 
thread  the  axle  into  the  hole  proviied  in  the  center  of  the 
electrical  connection  side  of  the  oan  buoy.  Tighten  the 
fitting   until   it      is    snug.       Secure    the    oan   buoy      to    the   deck 
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A -Frame 
Boom 


J 


Figure  B«1    Shipboard  system  component  layoat« 

winch  with  line  to  prevent  the  buoy  from  rolliag  when  t ne 
ship  is  under  way.  Do  not.  fasxan  tha  alectricai  connections 
at  this  point.  Position  the  optic  transmitter  agains-  tha 
instrumentation  subsection  (4)  fraaia  batween  the  two  spheres 
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(on  quick  release  clamp  siie)  with  ths  fiber  :;ptic  cable 
protrudir.g  from  the  power  supply  sph-re  end  of  th£  trans- 
mitter enclosure.  Ensure  there  is  sufficient  ro^m  fc^:  elec- 
trical cable  connections  it  the  P3M  sphere  i-a.i  -;f  the 
enclosure  and  securely  fasten  the  tc  ansmitter  to  r/-e  fraoie 
with    two    inch   Benthos   glass    spha re    s=aling    tape. 


alectrical    tape 


.r  ufabe  r 
wedge 


.stainless    steel 
banding  strap 


J 


Figure  B. 2 


Fiber  Optic  Enclosures  (1|  and  (2) 
reel  (3)  . 


i?.id  storage 


Step  U:  Connect  data  transmission  cable  (3)  to  the 
sensor  subsection  (2)  usiig  color  cdSs  for  easy  dad  error 
free  connection.  Prior  to  making  connection  apply  a  thin 
coat  of  vacuum  grease  zo  all  connector  shafts.  Make 
coniiecticns  and  secure  them  with  thi  locking  slaeves.  Do  not 
scrBw  locking  sleeves  too  tightly  or  damage  to  the  sleeves 
may  occur.  Fasten  the  othsc  end  of  the  cable  to  the  cables 
protruding   from  the   pulsa  cods   modulation  (PCM)    sphere. 
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Connect  the  PCM  data  cabls  to  ths  ^utic  transoQitter  data 
inpat  connector.  Use  car=  to  natch  the  dimple  on  the  cable 
connecter  with  the  indent  on  the  CDnnec Lor  housed  in  the 
transmitter  enclosure.  This  is  esseitic.i  with  ail  the  two 
conductor  Brantner  connections  in  ti-."=  system.  Failure  to 
follow  this  procedure  rfill  resai!:  in  darrage  to  the 
electrical  components  involved.  Coiinsct  the?  power  cable 
between  the  power  supply  sphere  and  the  optio  transmitter 
enclosure.  Now  connect  thi  power  cabls  oe-Tweer.  the  can  buoy 
and  the  receiver  enclosure. 

Step  5:  rest  for  idle  PCM  signal  by  ccnnecting  the  test 
cable  between  the  optio  raceivac  enclosurs  and  an 
oscilloscope.  Remove  tha  oower  actuating  magnet  on  the  top 
of  the  PCM  sphere.  Adjust  iihe  Dscillosrope  to  observe  a  200 
mv  PCM  signal.  The  signal  should  sync  rea.«;onably  well. 
Remove  the  actuating  magnet  from  tha  ssrisor  sphece  (one  at  a 
time)  and  observe  the  PCM  pulses  rua.  The  iiystem  is  now 
che::k.ed  out  from  sensors  through  tas  fiber  optic  link. 
Remove  xhe  optic  receiver  ia-.a  output  t-^z  cab'.a  and  install 
the  data  cable  between  thr  can  buoy  ind  the  ;)ptic  receiver 
enclosure.  Install  the  spar  buoy  cabLs  ('i'))  to  the  can  buoy 
and  test  for  the  PCM  signal  oy  pl>.  :;ii.-.a  the  oscilloscope 
probes  on  pins  2  and  3.  Seplaca  tha  r.aq.iets  oa  the  sensor 
sphere  and  the  PCM  sphere.  Disconnart  tr-=  spar  buoy  cable, 
the  optic  receiver  power  ::3.ble  and  tiie  optic  receiver  data 
cable.  The  system  is  read/  for  ieplDynant^ 

Step  1  through  step  u  n=y  ba  acoo nolished  along  side  the 
pier  or  en  route  to  the  iaploymant  Icoation. 

Step  5:  Fasten  the  lina  and  cable  to  the  spar  buoy.  When 
on  location  (130  foot  depti  for  this  axaaple)  and. the  ship's 
anchor  is  properly  set,  tha  ship's  oraw  will  fasten  the  spar 
buoy  to  the  on-board  crana  and  hoist  it  over  tha  starboard 
side.   When  the  fastener  is  at   the  nil,   a  tiabar  hitch  is 
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mada  about  the  spar  buoy  shaft  'dn-i  th9  spar  baoj  is  mada 
fast  to  the  ship's  rail.  Tha  crane  than  may  b*?.  ii.-icon?.scted. 
RemDve  the  four  1/4-20  has  head  bolts  froii!  the  top  of  the 
spar  buoy.  Remove  the  duramy  load  from  the  end  dZ  -he  roaxial 
cable  and  connect  the  cable  to  the  c^Dle  in  the  aiitenn.a  mast 
section  (9).  M ign  the  nating  inisx  marks  lad  mate  the 
sections.  Secure  the  sa^tion  in  place  v«ita  tha  belts 
previously  removed.  Tha  trench  utilized  for  thiij  prDcedura 
should  be  fastened  to  the  operator  wica  a  string  to  prevent 
loss  over  the  side.  When  tha  se^'ions  ace  secu:  i^,-  the  timber 
hitch  may  be  removed  and  the  spar  Duoy  lower?!  into  tha 
water  by  the  nylon  tending  line.  Whan  stable  in  the  water, 
remove  the  tape  on  the  spring  loadad  antenna  9i:::)und  plane. 
The  antenna  ground  plane  prasents  an  aye  hazz?.rd«  Use  care 
in  tending  the  spar  buoy.  The  spar  buoy  cable  should  be 
faired  along  the  starboard  siia  outboard  of  any 
obstructions. 

Step  6:  secure  the  sensor  cable  oDP.nectors  t:>  one  of  the 
sensor  sphere  hardhats  using  electrical  tape.  Remove  the 
actuating  magnets.  Allow  some  sialic  and  tape  :he  roaxial 
cable  (3)  to  tha  main  support  line.  I^a  sensor  sabsect-ion  is 
hoisted  up  and  the  A-frame  adjusted  a^t  until  r.h~  subsection 
is  clear  of  the  ship's  stern.  The  sensor  suL- sect  ion  is 
lowered  in  ten  foot  incraiants.  A  libaral  am:;j.-it  of  slack 
must  be  allowed  between  fixed  poiits  along  tee  coaxial 
cables  between  the  sensor  sabsection  a  .i  d  tne  instrumentation 
subsection.  Electrical  tape  is  applied  at  ten  f:>ot  intervals 
of  the  polypropylene  main  sapport  lina.  A  tab  is  placed  on 
the   tape   to   allow   for    easy   removal    daring   recovery. 

Step  7:  When  approximately  fiftsan  feet  of  the  cable 
remain,  the  A-frame  boom  is  adjustai  in  so  that  the  main 
support  line  is  very  close  to  the  stern  of  the  ship.  An 
ample    amount   of    fiber    opti::    cable      is    raeled    off    the    storage 
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reel  to  allow  the  inst  rums  nt2.tiD:i  SMbs^ctioi:  "o  bs 
positioned  on  the  stern  witi  PC)!  end  iDwn  and  giick  rdsase 
claaps  toward  the  line.  Ths  main  sapoort  line  is  ioweiid  to 
where  there  is  sufficient  slacic  in  ths  last  section  of  the 
coaxial  cable.  The  cable  connectors  sre  fastaaed  to  the 
Instrumentation  subsection  frame  with  electrical  tr.pe  and 
the  actuating  magnet  is  reiovei  fron  the  PCM  sphero,  Ths 
quick  release  clamps  are  male  fast  to  tr.e  aain  support  line. 
The  main  support  line  is  hoisted  off  the  deck  and  agai?  ths 
A-frame  boom  is  adjusted  aft  until  t'ls  load  is  clear  of  ths 
stern.  The  load  is  lowered  again  in  t=a  foot  inorements  and 
the  optic  cable  is  fastened  to  ths  niin  support  iir.-:-  with 
electrical  tape  as  before.  Care  uast  be  taken  to  prevent 
breaks  or  kinks  in  the  fiber  optir  cable.  Coord.-.n? tier, 
between  lowering  the  load  aid  raeling  off  the  fiber  optic  is 
a  mast.  When  the  170  foot  marker  (two  black  stripes  (15C 
ft. I  followed  by  two  succ=3sivs  singir  black  stripes  or  the 
main  support  line)  is  at  the  A-frame  boom  block,  the  vinch 
lowering  the  system  is  3t.opped.  The  Crmaining  fiber  optic 
cable  is  secured  to  the  storage  reel  with  a  rabber  wedgsi 
inserted  in  the  storage  groove  and  taped  in  place  with 
electrical  tape.  The  storage  C5  =  I  is  iffluobj.iizcd  by 
fastening  a  piece  of  parachute  oocd  from  tu3  ca:'i  buoy 
strongback  through  the  storage  reel  spokes.  The  pow^.r  and 
data  cable  connections  are  Bade  between  the  can  baoy  and  the 
optic  receiver  enclosure.  The  spar  baoy  tendiig  lir.v  is 
fastened  to  the  can  buoy  by  means  of  a  shackle  ~in.z  the  spar 
buoy  cable  is  connected  to  the  can  bjov.  The  nain  support 
line  is  lowered  until  the  170  foot  nark  is  at  deck  level. 
The  excess  fiber  optic  cable  nust  bs  carefully  tended  at 
this  point.  The  can  buoy  jiick  relsass  clamps  are  made  fast 
to  the  main  support  line.  Tie  Load  is  hoisted  and  the 
A-frame   boom   adjusted   aft    so    the    load    is    clear   of    the   stern. 
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The  can  buoy  is  lowered  intD  the  water  with  the  3xce3s  fibic-r 
optic  cable.  Simultaneoasl/  ths  spar  baoy  cabls  is  r-sl^assd 
over  the  side.  The  end  d£  the  polypropylene  r.ain  support 
line  is  unfastened  from  tha  winch  cible  and  fastened  to  a 
marker  bucy  then  released  Dver  the  strrn.  The  system  is  no*: 
deployed. 

Step  8:  Check  with  t^2  receiviaj  station  personnel  t;) 
ensure  they   are   receiving    a    sea   sida    PCM    signal.. 

Step  9:  The  recovery  Dperatioi  is  just  th=  reverse  of 
the  deployment  operation  omitting  the  systroa  checkoii". 
procedures. 
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APPENDIX    C 


A.       SYSTEM    TRANSFER    FUNCIIDNS 

The  triaxial  coil  n  agnetooiBt  5r  systsm  i=  dicgrarned 
in  Figure  C.I.  Presently  only  the  lani  based  system  employs 
three  coils,  the  ocean  system  uses  two  coils.  The  dati 
sensed  by  this  magnetometer  system  is  given  by  V!c(t),  Vy(tji 
and  Vz  (t) ,  the  output  vDltages  g  =  n?rated  by  the  input 
magnetic  fields.  The  magnetic  field  coraponeats  are  the 
desired  results  -herefore  the  system  transfer  function  xust 
be  determined  so  that  it  can  be  cenoved  froa  the  output 
voltage.  In  the  time  domain  the  relationship  between  the 
ou-put  voltage  and  the  magnetic  field  is  a  convolution 
integral: 


Voi(t)=  jhi 


hi(t-tM  bi(t')  it' 


(1) 


where    i=x,y,z.      It  is    more      convenient    to    express    equation    1 
in  the   frequency    domain   by   performing    a    Fourier   transform: 

Voi=Hi(u>)  Bi(cJ)  .  (21 

The   transfer      function   Hi (w  )       can      easily    be      determined    in 
teras    of   the   ratio  of    output    to   input. 

The  next  sec+ion  discjsses  the  components  of  the 
system  and  attempts  to  develop  a  system  model.  The  las- 
section  describes  the  calibration  aeasurements  and  the 
experimental   determinatiDn   of   the    transfer    function. 
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Figure   C.  1        Systsa    Comoonent    (Configuration 


2-       System   Description 

The  coil  magnetometers  are  composed  of  18  gauge 
copper  wire  with  ap proximi t ely  5463  turns.  rh=  individual 
transfer   function   for   the    coil   is: 

VcM=j    KB  M  (3) 

where    Vc  (u))       is   the    spectrin      of    the    /oltage    directly    out    of 
the   coil   and   K    is    a    gain    factor. 

The  next  component  in  -che  system  is  rhe  preamplifier 
low  pass  filter.  The  preamplifier  essentially  introduces  a 
gain  of  1000  so  the  vol-age  signal  can  be  senc  -o  remote 
instrumentation.         The      low      pass      filter      is      a      five      pole 


Chebyshev  with  a  3dB  cutDff  at  2D  Hz.  The  filter  is 
priaarily  used  to  remove  53  Hz  signals  from  the  system  and 
to  prefilter  the  data  before  it  is  iigitalized.  An  addi- 
tional pole  is  introduced  becaase  of  a  capacitive  coupling 
between  the  second  and  third  stages  3f  the  filter.  Figure 
C,2  shows  a  pole-zero  iiagram  of  the  preamplifier.  The 
transfer    function  takes  the    fcra; 

H(jw>)  =K/(Ja)+a)  (ju>+b)  ij^j*^'!^)  (JoJ+ojii  (jui+e-jf)  (juj^e+jf)     (U) 

where  K  is  the  gain  of  the  preamplifier  and  the  ooeficients 
a  -  f  are  taken  from  the  pole-zero  iiagram. 

The  preco nditionsr  has  additional  gain  of  approxi- 
mately 1000  and  has  zener  iiodes  to  Limit  zhe  input  voltage 
to  the  PCM  board.  Typically  the  cutoff  is  around  7  to  8 
volts. 

The  magnitude  and  phase  respDnse  of  the  preamplifi- 
er-filter and  precondition er  are  shown  in  Figures  C.3 
through  C.8.  The  input  test  signal  generating  ttiese  curves 
had  flat  magnitude  and  phase  character  as  shown  in  Figure 
C.9.  The  figures  then  represent  the  actual  transfer  func- 
tions of  the  preamplifier- filter  prer onditioner  system  for 
the  x-y-z  channels  of  lani  and  oceaa  sensor  systems.  The 
figures  show  that  the  channels  are  well  matched  in  magnitude 
and  phase.  The  phase  response  is  linear  and  virtually  iden- 
tical in  every  case.  This  is  an  important  point  since  the 
phase  part  of  the  transfer  function  is  not  removed  from  the 
output  voltage  signal.  However,  since  the  data  analysis  is 
in  terms  of  power  spectral  densities  and  coherences  this 
fact  has  no  bearing  on  the  results. 

The  total  transfer  function  of  the  system  is  the 
product  of  equations  3  ani  4.  Based  on  these  formulas  and 
the  figures,  the  resulting  transfer  function  should  be 
fairly  linear  at  the  low  frequency  eni. 
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Figure  C.2        Pole-Zero  Diagram   for   i!Ddsl    13-10A   Amplifier 
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Figure  C  .5    Response  of  the  Land  Z  System 
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3.   DeterminatijOn  of  zh^    System  Iransfsr  Function 

The  magnitude  of  tae  systsm  transfer  fanction  was 
detarmined  experimentally  through  the  use  of  a  Helmholtz 
coil,  a  source  of  sinusoidal  current  and  a  spectrum 
analyzer.  The  experimental  set  up  is  illustrated  in  Figure 
C. 10.  A  one  turn  Helmholtz  coil  with  a  0.61  meter  radius  and 
a  C.61  meter  separation  between  coils  is  used  to  provide  a 
calibrated  magnetic  field  to  the  sensing  ceil.  The  sensing 
coil  is  placed  concentric  to  the  3elmholtz  coil  and  is 
cer.tered  vertically  to  ensire  a  uaifDrm  magnetic  field  is 
created  about,  the  sensing  ::oil.  The  magnetic  field  in  this 
region  of  the  Helmholtz  coil  is  directed  along  its  axis  and 
is  given  by  [7]: 

Bj^=  ^SNI/sVs'r  (5) 

to  terms  on  the  order  of  (L/R)  where  L  is  the  linear  dimen- 
sion of  the  region,  N  is  the  number  of  turns,  R  is  the 
Helmholtz  coil  radius  and  I  is  the  driving  current  of  the 
coil.  The  driving  current  is  prcducai  by  a  Wavetek  signal 
generator  model  number  142.  The  cirrent  is  measured  by 
reading  the  voltage  across  %  995  ohm  r=sistor  id  series  with 
the  Helmholtz  coil.  This  voltage  is  measured  by  a  spectrum 
analyzer  (HP-3582A)  to  avoid  any  spurious  signals  generated 
by  harmonics.  Once  the  driving  currant  is  known,  then  the 
calibrating  magnetic  field  can  be  determined.  The  voltage 
output  of  the  sensor  system  is  then  measured  by  a  spectrum 
analyzer.  3y  repeating  the  measuren an ts  at  a  number  of 
different  frequencies  and  applied  fields,  the  system 
transfer   function   can    be    measurad. 

The  actual  calibration  was  conducted  at  the  La  Mesa 
village  site  using  the  satup  illustrated  in  Figure  CIO. 
The  Helmholtz  ceil,  sensing  coil  and  preamplifier  are 
located    approximately    40    meters   from      the    shack    which    housed 
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Figure   C.  10        Systsin    Calibration    ExperiisQt 

the  remainder  of  the  experiment.  with  the  Heimholtz  coil 
disconnec-ed ,  the  backgroan.i  magnetic  noise  of  this  location 
was  repeatedly  examined  to  snsure  a::oirate  readings.  Figure 
C.11  shows  a  typical  bacJcgcound  noise  plot  from  the  spectrum 
analyzer  directly  (therefore  the  zraasfer  function  has  not 
been  removed).  The  disriict  featuris  are  Schumann  reso- 
nances, otherwise  the  bao< ground  was  low  enough  and  stable 
eno'j  gh  -o  adequately  calibrate  the  system.  Only  at  the 
Schumann  frequencies  was  the  background  subtracted  from  the 
measured  output  voltage.  Zalibratioci  runs  were  conducted 
fron  0.05  Hz  to  20.0  Hz  rfith  applied  calibrated  magnetic: 
field  of  0.02,  0.2  and  1.0  nanoteslas  (gamma).  Figures  C.12 
through   C. 16      show  the   resulting      calibration   curves      of   the 


39 


land  and  ocean  magnetometer  systems.  rii=  true  transfer  func- 
tion is  represented  by  tha  0.02  nl  field  of  the  Helahol-cz 
coil.  The   larger      appiisi    field      siow      saturation    by      the 

system  caused  by  the  protective  zenec  diodes  in  the  precon- 
ditioner.  At  the  low  frequency  end,  all  the  systems  are 
closely   matched.       At    1    Hz    an    output    vDltage    of    1    v    indicates 


_j 


Figure   C.11        Background   Gaomagaetic   iJoise   f  (Hz)    05/20/82 
1513  hrs. 


a   measured   magnetic    field    Df    0.4    nT.    \    1    nl    fieli    is 
saturated   above    2    Hz. 

The  computation  of  the  actual  equations  representing 
the  magnitude  of  the  traisfsr  functions  presented  a  distinct 
problem   because    fitting   th=    data    to   ti2    theoreti::al    transfer 
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Figure   C.  12        Land   System   X-CdII   Calibration 

fun::tion  cf  the  previous  S53tic!i  is  difficult.  Ideally  this 
is  desirable  since  the  phase  character  cf  the  sysrem 
transfer  function  can  then  be  inferred.  However,  since  the 
priaary  area  of  interest  1=  in  the  iDwer  frequency  region, 
the  transfer  function  was  modeled  in  terms  of  linear 
segments.  The    equations      derived    ii      this      manner    have      a 

maximum  error  of  ^%  in  the  range  of  D  to  10  Hz  and  5%  in  the 
10   to    20    Hz   range. 

In  the  actual  computation  each  curve  was  broken  into 
convenient  segments  that  could  be  represented  by  a  straight 
line.  The  range  of  the  sagnents  were  0  to  3  Hz,  3  to  5  Hz,  5 
to  7.5  Hz,  7.5  to  10  Hz,  10  to  15  Hz  and  15  to  20  Hz.  For 
eaci    segment   the    slope   was    firsc    computed    by: 
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Figure   C-  13        Land   System   Y-Coil   Calibration 


M=  (Y2-Y1)/(X2-  XI) 


(6) 


Once  the  slope   of  the  line  was  found   the  intersection  with 
the  Y  axis  (b)  was  found  from: 


Y=MX+b. 


(7) 


These    computations  were   performed    for  each    segment    to    yield 

a  series      of    five   equations    ro      model  sach    sensing      coil    and 

amplification    system.      The    Lndividaal  transfer    fincticns   can 

be   found    in    the   computer   program,    see  Table    I. 
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Figure   C.  14        Land   Systsm   Z-CdII   Calibration 
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Figure  C.  1  5   Ocean  System  X-CdH  Calibration 
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Figure  C. 1 6        Ocean   System   Y-CdII   Calibration 
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TkBLE    I 
Transfsr    FuQcticn   3ompat2r    Program 


SEA* 

1  iFfFRQ.LE.15.)  GO    TO    2 

XXO  (L)  =XXO(L)  /  (-3  1.36+17.  06  6*FE?2"3.  5  3  89*  FRQ*FS  3) 
fYO  (L)  =YYO(L)  /  (-63.435  +  11  .7106*FSQ-3.  23  27  9*FRQ*FRQ) 
3  0    TO    8 

2  IF  (FRQ.  LE.10.  )  GO    TO    3 

XXO  (L)  =XXO(L)  /  (-18.486  +  i.  9029*FR2J 

YYOJL)  =YYO(L)  /  (-6  3.435  +  1  1  .  7  10  6*FRQ  -  D  .  23?79*FR2*FRQ) 

'^  n    T  n    fl 

3  rF(FRQ.LE.7.5)GO    TO    4 

XXO  (L)  =XXO(L)  /(3.9a7*FR2-9.35  8) 
YYOJL)  =YYO(Li  /  (4.  2  95*FR2-  11.837) 
GO    TO    8 

4  IFfFRQ.LE.S.)  GO   TO    5 

XXO  (L)  =XXO(L)  /  (2.  8  105*Fa  3 -0  .  5  657i 
YYOJL)  =YYO(L)  /  (2.  8  00*FR2-  .5  0) 

5  XXO  (L)  =XXO(L)  /  (2.723*FRQ| 
YYO  (L)  =YYO(L)  /  (2.7  17*FR2) 

GO    TO    8 

LAND: 

1  IF  (FEQ.  LE.15.  )  GO    TO    2 

XXX  (L)  =XXX(L)  /  (10  5.5-3.  n  =*FR3) 
YYY]L)=YYY(L)  /  (181.32-7.  5  83*? RQ) 

2  IF (FRQ. LE.10. )  GO    TO    3 

XXX  (L)  =XXX(L)  /  (5.  9  58*FR3-30.9  7) 
YYY|LJ  =YYY(Li  /  (7.  1  66* FRQ- 39  .  9  9) 
"^  n    td    8 

3  IF (FEQ. LE  .7.5)  GO    TO    4 

XXX  (Li=XXX(L)  /  (3.  49  2*FR2-6.31) 
YYY  (L)  =YYY(L)  /(4.252*FR3-10.35) 
GO    TO    8 

4  IF  (FRQ.LE.5.)  GO   TO    5 

XXX  (L)=XXX(L)  /  (2.631 1*FRQ +0.1  466  7 •> 
YYY  (L)  =YYY(LJ  /  (3.  0  12*FR2-  1  .  55  ) 
30    TO    8 

5  IF  (FRQ.LE.3.)  GO   TO    6 
XXX  (L)  =XXX(L    /  (2.6  311*FR2+0.14667| 


iil 


7      YYYJL)  =YYY(L)  /  (2.7G2*FR2 
"  0    TO     8 

6      XXX  (L)  =XXX(L)  /  (2.72*FRQf 


96 


APPENDIX    D 


A.       PULSE    CODE    MODULATION    SISTEM 

The    PCM   circuit    componsits    =rs    ii-jiited   on    a   3  inch    by    15 

inch    printed    circuit      boiri    wi-h    a    15    pin      dcubla  sidi^d   edge 

connector    at   each      end.      The    pii    assignnents      i-::L  the    signal 

connector   and  the  control   connector   ice    listed   iq  tabl.es    D.I 

and    D.  2,      respectively.      The    circuit    is    diagramed  in    ;:igures 
A. 4    and   A. 5. 


TABLE    II 
Signal    Connector  Pin    Assignment 


PIN  ASSIGNMENT 

— 5  cEanner""   T  moit  (anal^gl 

B  channel      2    inbut  (inalDgi 

C  channel      3    input  (inalcgi 

D  channel      U    inojt  (inalDai 

E  channel      8   input  (inalD'gt 

F  channel      7   irioir.  (analDg) 

H  channel      6   inpat  (inalcgi 

J  channel      5    inDit  (analDgl 

K  channel    15    inbut  (inalDgi 

L  channel    14    inoit  (analDg) 

M  channel    13    input,  (inalcgi 

N  channel    12   inoit  (inal^g) 

P  channel    11    input  (i".alDgi 

R  channel    10   inoit  (analDg) 

5  channel      9    input  (analogi 

1  -15V    (max.) 

3  common     (sigaali 

15  +15v    (max.) 


The  PCM  system  incorporates  a  crystal  oscillator  and 
associated  CMOS  integrated  circuitry  to  develop  the  clocking 
pulses,  a  16  channel  CMOS  iaalog  raulti  plezer ,  a  12  bit  CMOS 
analcg-to-digital  converts::  and  associated  circuitry  to 
provide  the    pulse  coding. 
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TkBLE    III 
Control    Connector   Pin    Assigninen' 


PIN 

ASSIGNMENT 

"T" 

EirrTTr~ITJ~(CHNL    14) 

2 

Wo 

3 

EXTDW2    IN 

H 

EXTDW2    SEL 

5 

SA!1P.    RATE    SEL    IN 

6 

DCL 

7 

POWER   GRD. 

8 

+  157 

9 

NRZ-    L 

10 

Bi0-L 

1  1 

EXT     (T 
CI        ^ 

12 

13 

C2 

ia 

C4 

15 

C8 

A 

TX    FILTER    INPar 

B 

SIS     GND 

C 

TX    FILTER    OaTPQT 

D 

SIG     G'^d 

K 

-15  7 

L 

+12  7    INPUT 

M 

SAiiPLE    RATE    123 

N 

SAMPLE    RATE    64 

P 

SAMPLE    RATE    32 

R 

SAMPLE    RATS    16 

S 

SAMPLE    RATE    8 

ThG  crystal  clock  oscillator  oo£ra*:23  at 
24.576  KHz  producing  a  sqiare  wave  Datpat  wi-h  3.  1c?l-  bit 
rata  of  one  bit  in  one  million.  This  palsa  trai^  with  ;ulse 
duration  of  20.35  niicrosscDi  ds  triggers  the  clock  inpat  of  a 
binary  up  counter  217(1/2)  and  is  also  rourei  directly  to 
pin  M  of  the  control  connector.  Pin  2  of  Z17  is  -he  counter 
enable  input  which  is  connected  to  7+  thereny  constantly 
enabling  the  counter.  Pin  7  is  tne  reset  input  ^tfhich  is 
grounded.  The  counter  is  composed  of  four  -oggle  flipflops 
thdc  are  incremented  on  the  positive  going  edge  of  the  clock 
pulses.  The  flip-flop  outputs  Q1  through  Q4  ace  routed  to 
control  connector  pins  N,P,R,  and  S,  respectively.  The  pulse 
repetition  rate  at  pins  M,N,P, R,  iii  S  represent  sample 
rates   of    128,    64,      32,      16,       and    8.       Selection    of    the    sample 
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rat=  is  dependent  upon  the  Jlyquist  sanpling  criteuiori.  "If  f 
is  the  highest  frequency  of  intarsst  in  ~hG  analog  dat--^ 
being  sampled,  than  the  simple  rat=  should  ne  equal  to  or 
greater  xhan  2f  for  unbiased  resoluti^a.  Curr-jQ-Iy  a  sample 
rats  of  64  per  second  is  employed  to  ensure  the  Nyquisc 
sampling  criterion  is  mat  over  ths  frequency  range  of 
intsrest       (0-20    Hz),  A   connection      between      the    pin      that 

represents  the  selected  simpling  rits  and  pii  5  cf  th^ 
control  pin  connector  will  establish  that  pu.\~s  repeti'cion 
rate    as   the   system  clock. 

The  clock  triggers  Z8,  a  fcur  stage  binary  up  counr.er. 
Z8  inputs  consist  of  a  single  clDck,  carry-in  (clock 
enable),  binary/decade,  up/dorfn,  pcaset  enabla,  and  four 
individual  jam  signals.  Foir  seperate  buffered  2  signals  -:\nd 
a  carry  out  signal  are  provide!  as  outputs.  y>  input  on 
pins  9  and  10  configure  ths  circuit  as  an  up  olnary  counter. 
The  jam  signal  is  selected  to  ievslop  the  data  ^Mord  length. 
The   counter   starts  at    a   coint      of   '4    as    confiaured    and   counts 


to  1 5  (12  counts).  CCO  goas  high  sitting  Z12.  The  prc-.s-st 
enable  signal  on  Z8  receives  the  hija  set  output  from  Z12 
(BTCR)  and  the  Z 12  counter  is  asynca ronously  preset  to  the 
jam  signal  count.  The  presat  count  is  gated  thcDugh  Z9  and 
Z10  and  the  convert  start  signal  tD  Z3,  the  15  char^nal 
analog-to-digital  converter.  Ihs    carry      out    signal       (CCO) 

froa  the  Z8  counter  acts  is  the  clock  input  to  Z17(1/2) 
another  four  stage  binary  up  countsr.  Its  outputs,  0^- 
through  QU  represent  the  multiplex  chmnel  snabir  signals  to 
Z2,  the  16  channel  analog  mult ipl  =  ic  3r .  Each  channel  is 
enabled  for  12  clock  pulses.  rha  analog  multiplexer, 
depending  on  the  state  of  ta =  multiplex  enable  signal, 
switches  a  common,  output  to  one  of  sixteen  iiputs.  Its 
output  is  amplified,  conditioned  and  passed  to  the  10  Volt 
input    of    Z3,      the   analog- to- digital    converter.      The    parallel 
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digital  output  of  Z3  is  passed  to  Z%  and  Z5  for  parailcl-to- 
serial  conversion.  The  analog  data  word  (ADW)  is  ga^ed  to 
NAND  gate  Z1 1  along  with  w3,  DCL,  aid  EXTWD2.  Wo  and  Wo  ara 
developed  in  NOR  gate  Z9(1/2).  Mo  is  aigh  when  Z^  through  C8 
are  zero  and  low  otherwisa.  Then  usiag  Z17's  30un"er  zero 
couQt  as  intial  sequencing  point,  tha  frame  s/nc  word  is 
passed  through  Z 13  pin  13  followed  by  channel  1  through 
chaanel      15    analog      data   wDrds.  Z12    and      Z14    develop      the 

signals  to  be  biphase  logic  or  noa-rsturn  to  zero  (IIRZ) 
logic.  This  system  feeds  the  biphass  logic  BI^-L  through  a 
low  pass  filter  and  a  signal  conditioning  circuit  where  the 
modilation  index  is  set.  En  ths  ssa.  system  tha  offset  and 
gai.1  potentiometers  are  aij  ustsa  sd  the  signal  conforms  to 
TTL  levels.  This  final  output  is  the  lata  input  to  the  optic 
transmitter. 
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APPENDIX    E 


A.       OPTICAL    FIBERS 

I*      Fiber   types 

Althoagh  many  types  of  fibairs  iiave  been  cried  in  the 
recanr  past,  presenx.  ta^tinoiogy  favors  three  classes  of 
fibers  for  optical  communications:  single-mode  step-index, 
multi-mode  step-index,  aad  milti-tnode  graded- index  (see 
Figure   E. 1) . 

a.  Single-Mode   Step-Index 

Single-mode  fibers  are  typlified  by  a  very  thin 
cor5  (1.5  to  6  micrometersi  with  unifom  index  of  refraction 
surrounded  by  a  cladding  of  considerably  greater  thickness 
whose  index  of  refraction  differs  by  a  very  snail  amount 
from  that  of  the  core.  Tie  "single  mode"  feature  can  be 
a-ctributed  to  the  abrupt  interface  between  core  and  cladding 
as  well  as  the  snail  changa  ia  the  refrac-ive  index.  In- 
cladding  material  is  lossy  so  that  any  modes  Lz  accepts  will 
guickly  attenuate  and  not  coupla  back  into  the  core.  Since 
the  difference  in  refractiva  index  from  core  to  cladding  :.s 
small,  the  numerical  aperture  is  snail  also  on  the  order  of 
0.14.  The  small  aperture  makes  coupling  to  souc ce/det ect or 
diff icul-  [3,  9]. 

b.  Jlulti-aode   Step- index 

The  core  diamat-er  of  multi-moda  step-index 
fibers  are  larger  (25-153  micrcmet ars) .  The  refractive 
indax    is    uniform    across   rha    cora    diamater   [9].      The    cladding 


101 


refractive  index  is  slightly  iDwer  than  that  of  the  cere. 
The  numerical  aperture  is  siall  but  ths  core  sizs  is  larger 
thaa  that  of  the  single  aois  fiber,  alleviating  somewhat  the 
coupling  problem.  The  loss  characteristics  of  these  fibers 
can  be  extremely  low  and  in.  general  3an  be  linked  directly 
to  the  lower  numerical  apertures  aad  manufacturing  tech- 
niques employed. 

c.   Multi-Mode  3raied-Index 

Graded  index  fibers  appear  to  be  the  optimum 
fiber  (of  those  currently  availiblei  for  long  line  ipplica- 
tioas  [8],  The  fibers  have  core  diameter  ranging  from  20  to 
150  micrometers  with  an  indsx  of  refraction  that  varies  as  a 
function  of  the  distance  from  the  core  center  as  shown  in 
Figure  E.2.  The  grade  in  tii=  index  irenoves  the  abrupt  change 
that  is  typical  of  the  step- index  fiosr.  The  net  effect  is 
that  all  modes  cross  the  ^pti^al  axis  at  aboit  the  same 
places  and  at  about  the  sam=  tin?  [10]. 

2«   ^i^snuaticn  in  Fiber  Wa/eg^uidas 

a.   Losses  in  General 

The  principal   :::auses  f^r  attenuation  losses 
are: 

1)  Material  absDrption:  Tiis  loss  can  be  mostly 
attributed  to  Hydroxy!  and  Transition  metal 
ions  in  glass. 

2)  Scattering:  This  loss  is  primarily  caused  by 
geometric  irregularities  at  the  core-cladding 
interface. 

3)  Radiation  losses:  This  Loss  is  a  result  of 
bending  the  fiber  cablB  at  small  cadii  of 
curvature  [11]- 
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Figure  E. 1        Tfpes  3f   Optical   Fibers 
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b.  Scattering  Lossss 

Ml  transpar  =  nt  materials  scatter 
due  to  frozen  in  thermal  fluctaatioas  of  constituent  atoms. 
This  intrinsic  scattering  represents  tiie  fundamental  contri- 
bution to  attenuation  in  optical  wav9gaides.  The  glasses 
utilized  in  fiber  waveguiies  generally  have  more  than  one 
oxiie  present  that  cause  a:iDther  fori  of  scattering.  Phis  is 
due  to  concentration  fluctuations  in  ths  constituent  oxides. 
For  high  silica  glasses,  3 oncentra ti on  fluctuations  typi- 
cally account  for  only  approximately  25^  of  the  observed 
scattering  loss. 

In  addition  to  xhese  tMO  intrinsic  scattering 
losses,  induced  scattering  through  ion-linear  affects  can 
exist  such  as  stimulated  Raman  aaa  Bcillouin  scattering. 
Because  of  the  snail  core  size,  the  confined  guidance  and 
the  long  interaction  length,  relatively  low  absolute  power 
levels  are  required  to  observe  such  affects.  Actually .  2000 
watts  has  been  injected  into  a  2  Km  liigth  of  guide  having  a 
75  nicroraeter  core  diametar  with  no  non-linear  attenuation 
obsa  rved  [11]. 

c.  Radiation   Lossas 

These  losses  are  associated  with  the  waveguide 
structure.  The  cladding  has  in  practice  been,  a  thickness  of 
a  few  tens  of  micrometers.  If  the  jacket  is  lossy  (to 
minimize  crosstalk),  soma  fraction  of  the  mode  field  can 
reach  this  region  and  be  attenuated.  Additional  radiation 
losses  are  experienced  when  the  optical  wave  guide  is 
curved.  The  presence  of  outsiie  parturbations  on  a  fiber 
cause  light  loss  through  tna  cladding.  In  step-index  fibers 
this  mechanism  causes  a  change  in  tha  angle  of  incidence  of 
the  light  rays  until  a  critical  loss  angle  is  reached.  In 
graded-index    fibers    loss    occurs    from    node    mixing    and   loss    of 
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the  higher  order  iiiodes  in  a  simular  fashion.  These  losses 
can  become  appreciable  in  fiber  cables  for  underwater  appli- 
cation. Losses  as  great  as  1400  IB/K.n  have  bsen  reported 
f  12]. 

d.  Loss    MinimizatiDn 

Fortunately  it  is  known  that  fibars  can  be 
protected  from  micr obeniing  and  ancased  in  CDastructions 
capable  of  high  pressure  [13].  The  size  of  thase  losses  is 
typically  about  1  dB/km.  MinimizatlDn  of  microbending  losses 
IS  achieved  by  encasing  tha  fiber  ii  alternataly  soft  and 
hard  jackets.  A  transversa  force  applied  to  the  stiff  outer 
jacket  will  be  resisted  ani  spread  o'/ar  a  larga  area.  The 
soft  inner  jackat  will  cDmoly  with  tha  force,  attanuating  it 
afid  spreading  it  over  a  still  largar  area.  The  total  force 
will  then  be  able  to  transmit  to  tha  fiber  which  is  itself 
dciff.  For  axampla  ITT's  T-1211  fibar  is  buffered  with 
STVISii,  a  silicon  rubber,  luring  tha  drawing  pr:)cess.  This 
has  the  advantage  of  protecting  tha  fiber  from  cDntaminants 
early  on,  particularily  OH  ions.  k  second  buffer,  Hytral 
72U6,  a  hard  polyester  elastomar.  Is  than  applied.  This 
combination  makas  the  fibar  a  gcoi  choice  foe  anderwatar 
cabling.  Hytrel  also  gives  the  fibeirs  hundreds  Df  hours  of 
inimanity  from  moisture  [13]. 

e.  Losses  due  to  Hydrostatic  Pressure 

Hydrostatic  pras suras  can  also  affact  tha  atten- 
uation in  optical  fibers.  Silicon  ribber  clad  silica  (PCS) 
is  compressed  in  sea  watar.  At  2000  psi,  the  numarical  aper- 
tura  (NA)  for  this  type  fiber  approaches  zero  and  no  light 
is  transmitted  due  to  compression  of  the  silicone  and 
increase  of  the  silicona  index  of  refraotion  [1^]. 
Therefore,  glass  clad  fibars  must  oa  used  in  underwater 
cabl es. 
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3-      Static   Fatigue  and   Strenaill 

a.  Since  the  fibsr  optic  utilized  in  a  data  link 
for  m.?.cnetic  data  collection  need  Qot  be  a  supporting 
elsaent  of  the  system  ths  cansile  strength  is  not  a  signifi- 
cant   fac-:or. 

While  the  theoretical  strength  of  optical  fiber  is 
mors  thea  600  !<:g/sq  mm  ]15]  the  actual  strength  is  related 
to  lainiuiim  strength  points  due  to  fl^ws  and  crarks  on  tha 
fib2r  s'\rfac=.  For  a  brsaking  strain  of  1  percent,  th  = 
ainimuQ  flaw  size  must  ba  lass  than  1  nicroraetar,  which  is 
difficult  to  assure  in  production  runs  '16].  Furthar  more, 
the  growth  if  flaws  is  ilrectly  ralated  to  tamperature, 
applied  strain,  relative  himidity  ani  aven  pH.  A  pH  of  7-8 
causes  the  mosT.  severe  flaw  growth  i17].  This  value  approxi- 
mates   the   pH   of   the    ocean. 

^'      ■'^la.ki.i  Configuration 

Cabling  of  optical  fibers  las  the  objective  of 
providirg  useable  strength  at  minimun  weight  while  ensuring 
environmental      protection.  Extensive    tes-s      upon      several 

different  configurations  of  miai-cablas  for  underwater  use 
witQ  and  without  conductors  has  been  conducted  by  NOSC, 
Hawaii  [18].  The  following  presents  some  of  these  cables 
whose  characteristics  may  nake  tham  suitable  for  underwater 
magnetic   data   collection   operations. 

figure  E. 2  repraseits  three  possible  cable  arrange- 
ments, dll  similar,  except  for  the  composition  of  the  load 
bearing  structure.  This  may  be  all  Kevlar,  a  hybrid  mix  of 
Kevlar    and   S-Glass,    or  all    3-Glass    in    an    epoxy    matrix. 

a.      Lightweight   Siigle   Fiber    Table 

1)  Cable      5358      (all      Kevlar) :  This      cable 

survived    500,000    flexure    cycles   at      23    percent    loading.       The 
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1. 

2. 
3. 
4. 


Low  loss,  graded  index."  optical  fiber  to  0,005"  O.D 

Annulus  of  Dow  Corning  silicone  (RTV  134)  to 

noraina 1  0.0  14"  0.0. 

Annulus  of  polyester  elsstomer  (Hytrel  7246)  to 

nomin  a  I  0.020"  O.D. 

Loadbearing  annulus  to  C.D.  of  0 

aay  be: 

a.  5058,  twelve  3gO-Oenier 

b.  Ten  38C.0cnier  yarns  of 
5  ENDS  HT5-901  S-Glass 

c.  32  ENDS  S-Glass  (nominally  6528  filaments) 


,04  0" 

yarns    of    )CEVLAR-49 
KEVLAR-49    and 


Figure  E. 2   Lightwsight  Siagle  Fiber  Cable 
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breaking    strength   of   this    ^able   is   225    pounds   with    a    tensile 
Strang th    cf    18  0,2  00    psi. 

2)  Cable  5383  (Hybrid) :  This  cable  survived 
2,370,000  fl€»xure  cycles  without  failure  at  20  percent 
loading.  Its  breaking  strangth  is  227  pounds  and  tensil=» 
streng-ch    is    16^0,0  00    psi. 

3)  Cable  5093  (all  3-5Lass)  :  This  is  the 
strongefe'-c  combination  with  a  breaking  strength  of  346  pounds 
and  a  tersile  strength  of  275,000  psi.  This  caale  survived 
750,000  flexi.r?  cycles  at  20  psrcsnt  loading  without 
breaking   befoi-e   the   test    was    stopped. 

Th.ose  cables  tasted  by  :J33C  Hawaii  provide  the 
ex^3nsion  of  fiber  optic  ts^hnology  to  the  underwater  envi- 
ronment and  tLose  raentionai  above  allow  for  utilization  in 
und9rwaL.er    magnetic    data    collection. 

B.       LIGHT    S0DLC2S    AND    DETSCrORS 

Light  sources  for  fibar  optic  systems  require  certain 
chacriCter:  sties  including  Id  ng-lif  e- t  i  is  in  use,  high  effi- 
ciency, reasonably  low  oost,  sufficient  optical  power 
output,  cap=jbiiity  for  various  typas  of  moduiaticn  and 
physical  compatibility  witii  fiber  s:ids.  There  are  three 
basic  corf igurations:  ligtit  emitting  diode,  solid  state 
lasar,    anci    sen.iconductor    injection    laser. 

"f-      Liakl  Sources 

?...      Light    Emitting    Dioda     (LSD) 

LED' s  generate  light  whan  carriers  injected 
across  a  p-n  junction  radiatively  racDDbine,  emaitting  radi- 
atiDn  over  a  solid  angle  of  27T  staradians.  The  actual 
output  power  amitted  into  a  unit  solid  angle  per  unit  emis- 
sion     area    is      low  for      LEO* s      aven    tiough      the    total      power 
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output  is  rsasor.r^bly  h:.gi.  Thers  ire  two  types  of  LED 
sources  available,  edge  Eimltters  ani  loned  surface  emitters. 
Doffisd  emitters  a  re  rather  complex  to  manufacture  resulting 
in  edge  eniittsr  LPJD  •  s  b^iiig  more  comnoaly  utilized.  Noise 
introduced  by  LED's  into  the  optical  fibers  is  relatively 
small  in  most  cases  an:l  is  proportional  to  the  spectral 
width  and  the  number  of  i:ransmission  aoies  propagated  in  the 
fib=r.  Typical  pcver  veirses  drive  response  curves  for  LED's 
are  usually  linear  (see  Figure  E.3).  Most  available  LED 
sources  have  peak  »aveli?ngths  of  33)  to  850  am  and  power 
outputs  of  appro  xij.sTely  1  nw.  LED*3  are  constructed  to  last 
in   excess    of    1    niiijion   h:)urs    [19/    8]. 
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50         KX>       150      20O      250      300      TSO 
DIRECT    CURRENT    (mA) 
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Dome  emitter, 
.4mm  dla 


Figure   E.3         Typical   Power   ts    Drive    Response   Curves 
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bo      Lasers 

Doufale-hetero janction  DH  laser  structures  can  be 
mads  with  relatively  long  lifeuimes  when  the  lattice  parame- 
ters of  the  electrical  composition  in  the  various  layers  and 
bouadaric-s   are      carefully    Batched.  In    the      DH    laser      with 

stripe  germetry,  the  lasing  region  occurs  only  in  the  trans- 
verse   direczion   under    the    tfidth   of   th=    contact   stripe. 

The  basic  semiconductor  mterial  for  both  LED*3 
and  injection  lasers  is  the  GaAs  combination.  Doping 
eleaenti?  like  Al ,  As,  In,  ^e,  and  Sb,  are  used  to  shift  the 
emission  wavelength.  The  threshold  current  for  typical  DH 
lasBrs  range  from  0.7  to  1^  kA/sq  en  and  the  quantum  effi- 
ciency from  1C  to  50%  compared  with  3%  for  an  LED  resulting 
in  a  larce  ir.cremental  increase  in  output  power  for  small 
increas?:=  in  current  beyoii  the  cirrent  threshold  (Figure 
E.4I  . 

c.       Ccmp  arisen 

Virtually  all  applications  of  fiber  optic  tech- 
nology for  local  data  coin unications  to  date  have  chosen 
LED' s  for  use  as  the  source.  Bandwidth  distance  rates  have 
been  within  the  capabilities  of  this  device  and  its  cheaper 
cost,  feai-e  of  use  and  op?Eational  axpsdiency  have  made  it 
competitive    with    lasers. 

2-      ?.S;i§ctors 

In  fiber  optic  systems,  the  most  conaonly  used 
receivers  utilize  photodiodes  (either  PIN  or  avalanche 
types)    to   convert   incident   light   into    electrical    energy. 
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a,  PIN   Photodiode    (PIN)  = 

Th«  PIN  dioie  consists  Df  a  relatively  large 
(very  ligh^.ly  doped)  region  sandwictiei  between  p  and  n  doped 
semiconductina  regions.  Ptiotons  absorbed  in  this  region 
create  elect ron-hcle  pairs  that  are  then  separated  by  an 
electric  field,  thus  generating  an  ilsctric  current  in  the 
loai  circuit,  I  ^i  e  quantum  sfficiency  of  the  photodiode  is  a 
fun;:tion   of    wave    length  ani    temperaturs    [3]. 

b,  Avalanchs    PhotDdiodas    (APO) 

The  APD  is  designed  for  applicatiois  requiring 
greater   sensitivity.  Beoause   of      a    strong      electric    field 

arising  within  it  as  a  result  of  external  biasing,  the  APD 
exhibits  an  internal  gain  n e chaaisn.  They  require  a  consid- 
erably higher  bias  voltage  than  PIN  diodes  (bias  voltages  on 
the    order    of    300    V  are  not    jncoiiiiion)     13]. 

c,  Ccraparision  of    PIN    and    APD 

As  r-.  result  of  the  high  reverse  bias  require- 
ments of  the  APD  and  the  power  constraints  of  the  inagnetic 
data  collection  system  the  PIN  diode  is  the  obvious  choice 
for   the   system '3    detector. 


C.       TBANSfllTTSri/RECSIVrlR 

There  is  a  v^ide  variety  of  o ommer^ially  available  trans- 
mitter/receiver sets  with  tie  full  range  of  source/detectors 
previously  mentioned.  These  preolide  the  necessity  of 
designing  the  associated  circuitry.  The  Burr-Brown  3713T/R 
is   such   a    set. 
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D.       SPLICING    AKD    CONi^  ECTIN5 

A  prcblsm  in  Eany  appl  L::at  iDns  is  in  joining  one  length 
of  fiber  or  cable  to  anot-hsr,  either  by  splicing  (relatively 
perianent)  or  rising  ^oiaectDrs  (intended  for  repeated 
disconnection   and   reconneciion) . 

1.      5£licin3 

Improvements  in  n  anuf  actar ing  techr.igae'=;  that 
continue  to  lower  fiber  atienuation  lave  increased  the  need 
for  low-loss  field  splicing  methods  tiat  will  pernit  greater 
spacing  between  repeaters..  Ths  coipLing  loss  of  a  splice 
depends  on  two  types  of  fa::tors:  thoss  that  are  intrinsic  to 
the  fibers  being  spliced  a:id  extrinsi::  ones  which  depend  on 
the  splicing  tecaniquG.  I:itrinsic  factors  inclai?  core  dia- 
meter, numerical  aperture,  and  rsfractive  index  profile. 
These  can  be  expected  to  vary  ran.iDJily  because  they  are 
influenced  by  the  iiianu  :  a  ctur  ing  and  quality  control 
processes.  Their  effect  on  splice  loss  can  be  treated 
statistically.  The  factors  introduced  oy  the  splicing  method 
include  lateral  in  isalignaeii  t  ,  end  separation,  and  end  prepa- 
ration quali-y.  Of  these,  lateral  misalignment  is  the  most 
critical. 

The  ideal  splicing  method  is  a  field  usable  tech- 
nique which  miniiLizes  b  hese  fiber  extrinsic  factors. 
Splicing  methods  now  -.ivailible  depeni  on  either  fusion  cr 
adhesives    [20']- 

a.      Fusion  Method 

The  fusion  or  "hot  splice"  method  uses  an  elec- 
tric arc  or  other  heat  source  to  fuse  two  fiber  ends 
together   after   they   have   bsen   accurately    aligned. 
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b.       Adhesive    Methoi 

In  t  hs  adhc5i(*€  Dr  mechanical  approach,  a 
mechanical  deyics  is  usei  to  aliga  th=  fibers  thcDugh  refer- 
ence tc  their  outer  diaaefesrs.  The  fibers  then  are  bonded 
into   this    position  permanen uly,    usually    with   an   epoxy 

adhesive. 

2«      Connecting 

The  same  intrinsic  and  extrinsic  factors  apply  to 
the  design  of  fiber  optic  ::3nnec-Drs  with  the  added  require- 
ment that  the  connector  system  peralt  numerous  disconnec- 
tions and  rematings  withojh  an  a ppi:e:;iable  increase  in  the 
coupling  loss.  There  are  -wo  basic  types  of  ::onnectors: 
those  which  take  the  aiigaaent  reference  from  the  outside 
diameters  of  the  fibers  thansel/es,  aid  those  whi::h  encapsu- 
late the  fibers  then  refer=nce  the  Dater  diameters  of  the 
encapsulating   mem-bers. 

a.  Capillary   Type 

The  mcsT.  cbvio.i  s  design  based  on  alignment  of 
the  fibers  thiSiselves  li.  tiie  mating  Df  -wo  fibers  inside  a 
precision  capillary  with  a;.air.eter  only  slightly  greater  than 
that  of  the  fibers,  a  vari-:t:icn  of  tha  capillary  approach  is 
the   quasi-capillary    formed   oy   three    preoision   rois. 

b.  Other    Types 

The  other  commercial  terminations  fall  into  the 
second  category.  Some  use  a  molded  taper  or  cylindrical 
ferrule  to  align  one  fiber  rfith  another.  Most  of  the  avail- 
able connectors  which  use  mechanical  alignment  Df  machined 
members  are  "dry"  (they  do  not  ase  an  index  matching  fluid), 
and  many  are  based  on  the  3 MA-type  connector  ieveloped  for 
radio    frequency   transmissiDi    over   cDpper    coaxial    cable   [20]. 
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E.       OPTICAL    DATA     LINK    ENSINESRIilG 

The  fcllowing  paragraphs  oitiitir  procedurss  for  engi- 
neering a  fiber  optic  iata  linli  utilizing  th=  Burr-Brown 
3713T    transmitter    and    371 3R    rec3iver    [21]- 

1«       Determining;   the   3 able   T/pe 

Initially,  the  primary  .ia^t^r  which  determines 
subsequent  design  of  a  fibsr  optic  system  is  the  amouni  of 
attsnuaticn  expected.  Fibers  have  bsen  developed  which 
exhibit      attenuation    of      lass     tha-:      1     iB/kra.  Comasrcial, 

fibers  are  available  with  attenuatioi  of  less  than  6  dB/kra 
a-  3.85  micrcaieter  [22].  Tia  rugq^dizid  graded-iadex  optical 
cables  that  were  tested  by  NOS: ,  Hawaii  had  attenuation  of 
approximately  7  d3/k.m  [23].  In  aost  applications  the  cost  of 
the  fiber  is  inversely  proportional  tD  .ts  att  eaiiat  ion .  For 
this  reason,  selection  of  i  fiber  is  nore  directly  related 
to  dollars  than  d5  [22],  Thrre  ara  also  envircnmental 
concerns.  The  mechanical  stresses^  ^nyironmenral  temperature 
extremes,  and  exposure  to  airsh  cneiioals  musr  be  considered 
when  choosing  a  particular  -able.  Sach  things  as  jacketing 
material  and  the  type  of  added  strangtii  menbers  (see  Figure 
E.  2)     determine   a    cable's   me^h?.  nical    aai    envircnnantal 

capabilities. 

There  are  three  basic  ryp js  of  optical  fibers  to 
choose  frca:  step  index,  single- mod  2,  and  graded-index.  Each 
type  specifies  the  profile  Dr  variation  of  the  fiber's  index 
of  refraction  as  a  function  of  radial  distan::e  from  the 
fiber's  center  as  previously  diagrammed  in  Figure  E.I. 
While  step-index  fiber  with  a  irore  diameter  of  200  (iiicrome- 
ters  or  larger  is  recoimeaded  for  ase  with  tha  3713T  and 
3713R  as  it  is  the  best  compromise  between  cere  diameter, 
bandwidth,      coupling    loss,      and    cost,       the    graded    index    type 
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5083  (Hybrid)  cable  prrviDusiy  issrribed  is  iacorporatad 
intD  the  subject  system  design.  Ihh  outside  diameter 
disparity  is  reduced  by  dDplicition  of  several  jackets  of 
heat  shrink  tubing.  For  short  ieigths  of  fiber,  the 
increased  coupling  loss  assDciated  '^fith  the  strallsr  diameter 
cors    is  tolerable, 

2.      Determining   L^unchsi    £2Si^i 

Once  the  mschanical  cha  ract  sc  if:  t  ics  and  rhs  cable 
xypa  have  been  chosen,  it  becomes  la-s-sary  to  analyze  link 
performance  with  various  :;abl=s.  Vhb  first  analysis  to 
undertake  is  the  power  lauacjhed  into  a  particular  cable  from 
the   optic  transmitter. 

Beginning  with  tis  LSD  soiree,  four  important 
optical  parameters  must  be  known  in  or.i~r  to  caloulat*^  accu- 
ratsly   the   device's    intersotion    with    th2    system.    These    are: 

a.  Total    optical    po'^er 

b.  Peak    wavelength   of   emission 

c.  Physical    size    of    output    poet 

d.  Spatial  profili    of    tna    aralssior    pattern 

Total  output  powar  ( Po)  is  asaally  gi\ien  in  units  of 
microwatts.  For  convenianoa  tha  optio?.!  power,  as  well  as 
the  losses  throughout  the  sysca-ji  i.il  the  detaotor  sensi- 
tivity, can  be  expressed  ii  dS.  The  total  system  loss  is  tha 
sum  of  these  individual  components.  To  oonvert  taa  L2D  total 
powar    output   to    dBm,    refai    it    to    a    la/al    of    one    uilliwatt. 

LED    output    power  (dSra)     =    10    log  (Po  ( uw) /10G0)  (1) 

The  peak  wavelength  of  aaissioa  ^  is  necessary  for 
determining  the  loss  of  tia  system  fiber  (which  is  a  func- 
tion of  wavelength).  It  is  not  unusual  for  fioar  attenua- 
tion to  double  or  even  triple  ovar  only  a  13D  nanometer 
range    in    wavelength. 
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A  powsr  loss  can  occur  if  tbsre  is  a  mismatch  of 
areas  of  the  emitting  port  anii  thi  fiber  sria.  This  loss 
expressed   in   dB,     is: 

LED-to- fiber   area   los3(dB)    =    231d  g  (D2/D1;  (2) 

where  D1  and  D2  are  the  ilamsters  of  the  systsa  cores  and 
the  LED  output  port,  respeot  ivel  y .  If  D1  is  grc-fi-er  than  D2, 
the   area    loss   is    zero. 

The  spitial  emission  pattsri  of  ^-r\e  '..SD's  output 
must  be  known,  because  a  loss  can  orcur  if  tr.B  high-crd==r 
moda  rays  of  the  LED*s  oacput  cannot  be  captnvad  into  the 
numarical  aperture  (NA)  of  the  systsm  fiber.  This  is  called 
NA  loss.  Most  data  shaets  give  information  oi  numerical 
aperture.  If  the  NA  of  the  LED  is  lass  than  that  of  tha 
fibar  then  no  NA  loss  ocoirs  at  that  interfi::e.  If  tha 
system    cable  has    the    smailar    NA,     than    the    loss   is: 

NA   loss  (dB)     =  201og(NA    sourca/NA    system   fib=jr)  (3) 

If  there  are  no  fibar  splicss  or  taps  in  the  system 
the  remaining  losses  will  be  =t  tha  fibc-r  aa':r::tcr  inter- 
face.   These    losses  may   ba    oalculated    is    described    above. 

Typical  gap  and  misalignmant  displac^^nant  in  tha 
various  connector  systems  available  for  the  fiber  optic 
active  component  package  are  D.15  am  and  0.05  mm  respec- 
tively. Angular  misalignraant  in  th=s~  connecters  is  very 
small    and    can    be    ignored. 

Beginning  at  the  da  tec  tor  and  working  na::k  wards,  the 
reguired    LED   drive   current    oan    be    calculated   as    follows. 

Generally  receiver  sensitivity  is  a  function  of 
banlwidth  and  desired  bi t-a rror-rata  (BER) .  The  detector 
data  sheet  specifies  tha  minimum  powar  to  be  launched  into 
the  detector  input  port  foe  a  given  BER.  Therefore,  it  is 
necessary   to      calculate   tha      coupling    loss      fro:n    the      system 
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fiber    to      the   detector.      Ttiis  will    yisld   the      opticr.l    Dow  =  r 

required   at      the    detector    sad  of   ths      fiber.      It    is      "hsn   a 

simple    matter   to    subtract    the  loss   of    the    cable    to    cirrivs    at 

the   required   launched    power.  Detector       ininimum   power    can   be 

converted   to   dBm    by    equation  (1).      ::aDie    loss   is    calcalated 
by      multiplying      cable      Leagth    (km)          by      cable      it T.^nuaticn 
(dB/km). 

PL    =    PD.M    -    NALD    -    ALD    -    CL    -    FL    -    3lB  (U) 

where, 

PL   =    launched    power 

PDM   =    minimum   detector    power 

NALD   =   detector   NA    Loss 

ALD   =   detector   area,    loss 

CL   =    cable  loss 

FL  =   Fresnel    loss 

For  design  ease  the  3713r/R  speirif ici t ion  shee-  in:;ladis  ti:3 
graph  illustrated  in  Figire  E.  U.  Power  l3iinohed  may  be 
determined   simply    by    enteririg   rhe    gripti    with    th2    oor--^ 

diaaeter. 

3-      Determining   Haximui    Cable   Lenjth, 

Once  the  actual  laiachei  power  in-o  a  sppoific  cable 
(PL)  and  the  attenuation  of  a  fiber  a.t  a  specific  wavelength 
(L  I  are  known,  the  maxiaiim  cable  length  say  be  obtained. 
First    the    loss    margin     (LM)    is    found    b/: 

LM(dB)    =    10   log(PL/Pin,mini  (5) 

Pin, min  is  the  minimum  rited  ourpit  power  (with  power 
adjusted      set    for      maximum      power)  isually    in      nicrowatts. 

Pin, min  is  specified  in  ths  transmitter  specif icition  sheet. 
Then  Xmax,  the  maximum  cable  length  that  will  jusx  present 
Pin,  min  to   t,he   transmitter   is   found    by: 
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Figure  3.4    Launchsi  Power  vs.  Core  Diameter 


Xmax  =  LM  (d3) /I^(i3m)  (6) 

If  this  length  is  too  short,  a  cable  with  a  largsr  diameter 
(larger  PI)  or  less  atteaaation  (sialler  Lv)  dt  a  trans- 
mitter   with   a    higher    Pin,asit\    must    be    isad. 

The  type  5088  cabla,  foe  example,  has  L^  =  7dB/kip  at 
660  nm  and  a  core  diameter  0.005  iioh  or  127  liorometars 
whizh    corresponds   to    0,9    mi::  rcwa  tts.       Pin,min   =    30    nanowatts 

-9 

worst    case    for    10         BER.    This: 


LM    =    10lDg(0.9    rni^rowitts/SD    nanowatts' 
LM   =    101og30   =    r4.  77dB 
Xmax    =     14.  77d3/7x1  O'^dBm 
Xmax   =    2.  1  1    km 


(7) 

(8) 

(9) 
(10) 
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This  length  is  loncrsr  '-.har.  requira:!  and  is  therefors 
sufficient. 

U,      Determininj    Miniaiuis    Cable    Lsnith 

Short  C2,bl€  lengths  wi^h  large  diameter  fibers  may 
cause  some  receiver  slew  rate  liirit  vrhich  will  appear  as 
duty  cycle  distortion,  Li.iiit.iDg  "hs^  power  couplsd  into  the 
recsiver    from   the   cable   will    pr=v3nt    th:-S    from   occuring. 

The  minimum  cable  length  nay  be  deterained  for  a 
specific  cable  by  calculating  'he  powe::  launched  into  it  by 
usiag  the  typical  curve  fr^n  Figur?^  S.U.  The  typical  value 
of  launched  power  should  d=  us2d  sii::.-;  this  will  give  the 
most  power  into  the  recei/er,  cssjlting  in  the  desired 
minimum  cable  length.  Continuing  tha  example  of  the  type 
5083  cable,  che  typical  launcned  oowsi:  is  0.9  microwatts. 
The  minimum  loss  margin  is  found  using  the  typical  power 
launched  into  the  cable  and  th^  iia.cimum  input  to  the 
receiver  which  is  1  niicrowat't  f.-?r  the  3'M3R  as  stated  in  the. 
electrical   specifications    fDr   the    r.iC2iTer. 

LM   fflin   =    lOlog    ( ?l/?in, uaxt  (11) 

LM    min   =    lOlog    (0.9    microwatt/I    mi-ro wii.z*)     =   -0.46dB       (12) 

The  minimum  cable  length  for  :.  he  sp3;:ified  cable  in  this 
example   is: 

Xmin    =   LH,rain/L^  (13) 

Xmin    =    -O.a6dB/7x10"^    dBs    =    -65.7m  (lU) 

If  this  minimum  cable  langth  is  too  long  for  a  particular 
application,  use  a  smaller  3ore  diametsr  cable,  a  cable  with 
a  higher  loss,  or  the  pow=r  launched  into  the  cable  may  be 
reduced  as  described  in  ssotion  II.^.I.c  by  adjusting  the 
transmitter  power  adjust.  The  negative  sign  in  -he  minimum 
length    indicates      there  is    lo      niniajii    length      limitation   on 
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this    fiber      optic   cable      -^hsv.    used      ii    '":on junctiDn      with   the 
3713T/R    set. 

F.       TYPICAL    APPLICATIONS 

The  3713T/R  fiber  optir  link  s.-JlViS  such  data  transmis- 
sion problems  as  crosstalk,  cirxgiag,  aad  echoes. 
Ele::tromagnetic  radiation  in terf arsncs  i?  avoided  when  using 
a  fiber  optic  data  link  in  higi  noise  sivironments. 
Ligtitning  damags  to  cablas  and  coanaotid  eqjipmsnt  can  be 
eliminated  where  fiber  opti;  cables  r?pl?c9  mstaliic  conduc- 
tors. Fiber  optic  cablas  ire  not  sibie.ct  to  shorting  out 
when    emersed    in    water. 
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APPENDIX    F 
TSSr     PRiGRA)! 

//J0EM1MB    JOB    (1592,0129)  ,«  FISHER    SMC,  1  399  •,  CL4SS=C 
//*:1AIN   OBG=NPGVM1.15  92P,LISES  =  (5  0) 
//^FORMAT    PR,DDNAJ1E=PL0T.  SZ3  VECP  R,  DESI  =  LOCAL 

//   EXEC    FRTXCLGP,  FARM.LKED='  LIS!  ,  M  A?,  !C  REF'  ,  REGI^U  .  GO  =2700K 
//SySUTI  DD      aNIT=SYSDA,3PACS= (CYL, (3,3)) 

//SrSUT2         DD      0Nir=SYSDA,3? ACE= (CIL, (8, 8) ) 
//SfSLIN         DD      SPACE=  (6080,  (80,30) )  ,aSir  =  SYSDA 
//FDRT.SYSIN    DD    D  SN=MSS.  5  YS3  .  NONIMSL.  5  DQRCE  (FCUH?)  ,D:CSP  =  SHR 
//  DD       * 

C  A    LARGE   NDMBER    3F     ARRAYS    ARE    JTILISSO 

C  WITH    THE    INTENTION    OF    dAINTAINING    EACH 

C  BIT    OF    ORIGINAL    DATA    OR    CALCaLATION, 

C  READY     FOR    RECALL     AT    ANY    TI^E    FOR    FaRPtlER 

C  COMPUTAITON.    EQOrVALENT    ARRAYS    ARE    EMPLOYED 

C  BECAOSETHE    PLOTTING    ROUTI^IE    •DPAWF»     IS 

C  NOT    ABLE    TO    HANDLE    COMPLEX    ARRAYS. 

C  ARRAYS    'ITS', 'RTB  •  ,  '  ALAB*     ^  >I  D    '  TITLE  ^'     ARE 

C  USED    TO   GENERATE    THE    VSRSATEIC    PLOTTER    DUTPDT. 

C  OCEAN     ARRAYS    ARE     INDICATED    VIA    THE    SbTFIX 

C  'O',     AND    LAND    ARRAYS    ARE    IllDICATED    WITH     'L'. 

COMPLEXES     XXO  (8192)  ,y  YO(3192)  , 

C    CSO(8192)  ,  COO  (8192)  ,ZXO  (3  192)  , 

C    CRO  (8192)  ,CTO  (8192)  ,3U0  (8192»  , 

C    C70  (3192)  ,CWO  (8192)  ,:P0(8192>  , 

C    ZYO(8192)  ,COMOL(8192|  ,COL(8192i  , 

C    CUL(8192)  ,CyL(8192)  ,:WL  (3192)  , 

C    COFOL  (8192), XXL  (81921   ,YYL(8192l  ,ZXL(8192), 

C    ZYL  (819  2)  ,  CSL  (8  192)  ,  3  PO  (3  1  92)  , 

C    BMO  (8192)  ,  BPL  (8192)  ,3ML  (8192)  , 
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C    SPOL(8192)  ,SMOL(8192)  ,SPD(8192j  , 
C    SPL(8192)  ,  SMO  (8192)  ,3TL(3192|  , 
C    CPL(8192)  ,CRL  (8192)  ,SML  (3192) 

DIMENSION    TIME  (8192)  ,  FREQ  (8192i  , 
C    WORK(1638i*)  ,FRQ2(3192) 

INTEGER*!*     ITB(12)  /12*0/ 

REAL*4  RTB  (28) /23*0.0/, 

C    CROO(16384),CPCO(1633a)  ,3UOO{1S334)  , 
C    COOO(1638a),CTOO(1633a)  ,C  VOD  ( 1  5  38  !4)  , 
C    CWOO( 16384) ,CTLL(16334) , 
C    CPLL(16  38  4),CRLL(163  3a) ,:OLL(15  384) , 
C    CULL( 16334) ,CVLL(16334) , 

C    CWLL(16  38  4),COMOLL (15  38  4)  , C3?DL L (  1  6384)  , 
C    SPOO  (16384), SPLL(16334)  , 
C    SMOO(16384) ,SMLL(16334) 

REAL    aLA3 (22) /'PDXO' ,  • PD?0 • , ' ST  1 3 '  , 
C    'ST20', 'ST30'  , 'srOO' ,  'PHA3 », 
C    'GAMO', «?DXL' , 'PDYL' ,  »STlL»,«Sr2L', 
C    'ST3L» , '3T0L' , '?HAL' »  'GAML* , 
C    »5PL0', '3?  IL'  , 'SMDO'  ,  »SMLL» ,':?3L' , 'CMOL'/ 

REAL*8    TirLS(12) 

EQOIVALSNCE  (TITLE  (1)  ,  RTB  (  5)  )  ,  (3  ROO  ( 1 )  , 
C    CRO  (1)  )  ,  (CPOO  (1)  ,CP3(  1)  )  , 

c    (cooo  (1)  ,coo(i) ) ,  (crDO(i)  ,crD  (D  i  ,  (cuoo  (ii  ,cio  (^) )  , 

C  (CWOO  (1)  ,CWO  (1)  )  ,  (C7D0(1)  ,073  (1)  )  ,  (CTLL(1)  ,  CTL  (' )  )  , 

C  (CPLL  (1),CPL(1)),(CSLL(1)  ,CRL(1))  ,  (COLL  (1)  ,COL  ( '0  )  . 

C  (COLL  (1)  ,CUL(1)  )  ,  (C7LL(1)  ,C7L(1j  )  ,  (CWLL(1)  ^C'^Li^))  , 

C  (COPOLL  (1)   ,CDP0L(1) )  ,   (CCM0LL(1|  , 

C  COMOL  (1)  )  ,   (SP00(1)  ,3P0  (1|  )  , 

C  (SPLL  (1)  ,  SPL(1)  ),  (3.13  0(1)  ,SM3(1)  )  ,  {SMLL(1»  ,SML  (1)  ) 

DATA    XXO,YYO, CS0,C03,  ZXO, Z YO, CR 3 , CTO , CUD, C70 , CWO , 

C  XXL,YYL,ZXI,ZYL,C3L,  3PO,BM0,BPL, 

C  BML,SP0L,SM0L,S?O,SPL,SM0, 

C  SML,CTL,CRI,C0L,C7L,C0PCL,CWL, 
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C    COPOL,COMOL/278523*(D.0,3. 0) / 

DATA    TIME, FREQ,FR22/2a575*0-0/ 

IFEAME=8192 

NR=30 

FNE=FLOAT (NR) 

P=0.0 

DO    70    L1=1  ,NR 
C  THE    DO    LOOP    SNDI.'JG    WITH    SriPEMENT    70    ENABLES 

C  THE    PROSRAI'!    TO    PROCESS    A    LiRGE    AMOUNT    3?    DATA 

C  3Y    REPEATING    THE    PROCESS    I^    BLOCKS. THE    DATA 

C  POINTS    FROM    EACH     RON    THROQ^H    THE    DO    LOOP    ARE 

C  ADDED    TOGETHER    AS D    EVENTJALL?     AVERAGED    BY    THE 

C  NUMBER    OF    RUNS    THROUGH    THE    DO    LOO?.     »NR' 

C  REPRESENTS    THE    SJMBES    OF    DITA    SEQUENCES    TO 

C  BE    AVERAGED.    1     32QUE:^CE    CURRENTLY    EQUALS 

C  8192    DATA    P0IST3     FOR    EACH    CHANNEL    OH    256 

C  SECONDS   OF    DATA,     FOR    THIS    TEST    SINE    AND 

C  COSINE    WAVES    DE    4  HZ    WIL    3E    USED. 

DO    60    JJ=1  ,IFRAME 

XXO  (J J)  =SIN(P*.785  39  3  1631 

YYO(JJ)  =COS(P*.735393  163| 

XXL  (J J)  =SIN(P*.  785  393  16  3) 

YYL  (J J)  =COS(P*.785  39  3  163l 

P  =  P+1 
5  0       CONTINUE 
C  THE    FOLLOWING    SEITION    GENERATES    THE    TIME 

C  AND    FREQUENCY    ARRAYS    AND    NORMALIZES 

C  THE    INPUT    PCM    DATA    TD    V0LT\3E    FORM    IN 

C  PREPARATION    FOR    THE    FAST    FDURISR    TRANSFORM 

C  TO    THE    FREQUENCY     DOMAIN. 

N=8192 

FN=FLOAT(N) 

DELTAT=1,/32. 

T=FN*DELTAT 
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DELTAF=1./T 

DO    20    J  =  1  ,  N 

TIME  (J)  =DELTAT*FL3Ar(  J) 

FEEQ(  J)  =DELTAF*FLDAr  (J) 

FRQ2(J)  =AL0310  (FREQ(J)  ) 
20      CONTINUE 
C  THE    NEXT    FOUH    STATEMENTS    PSSFORM 

C  AN    FFT    ON    THE    li^IPUT    TIME    SERIES 

C  DATA. 

CALL    FOORT  (XX0,N,1  ,-1  ,0,;?ORK) 

CALL    FOURT  (yY0,N,1  ,-1  ,0,;*ORK| 

CALL    FOURT  (XXL,N,1  ,-1  ,0,MORK> 

CALL    FOURT  (Y'/L,N,1  ,-1  ,0,-^OHKl 

DO    40    Ka=1 ,N 

XXO(KU)  =XXO(Ka)/FN 

YYO  (K4)  =yYO(KU)  /FN 

XXL(KU)  =XXL(KU)  /FN 

YYL(Ka)  =YYL(KU)  /FN 
'4  0       CONTINUE 
C  THE    NEXT    LOOP    CALCULATES    THE 

C  UNNORMALIZED    SPECTRAL     DE^^SITIZS 

C  FOR    SINGLE    SITS    DRTHD30NAL 

C  COMPONENTS    OF    TiE    GEOMAGNETIC    FIELD, 

C  THE    INDIVIDUAL    SITE    :R0S    SPECTRA 

C  BETWEEN   COMPONENTS,    THE    INDIVIDUAL 

C  ORTHOGONAL    COMPONENTS    OF    RIGHT    AND 

C  LEFT    CIRCULARLY    POLARIZED    SPECTRA, 

C  AND    TWO   SITE    CR333    SPECTRA    FOR    RIGHT 

C  AND    LEFT    CIRCULAR    POL  ARIZ  AT  IDN. 

DO    30    11=1 ,N 

ZXO  (II)  =ZXC(II)  +  (XX3  (II)  4'C0NJG(XX0  (II)  )  ) 

ZYO  (II)  =ZYO(II)  +  (YYD  (II)  *CONJG  (YYO  (II)  )  ) 

CSO  (II)  =CSO(II)  +  (XXD  (II)  *CONJG  (YYO  (II)  )  ) 

ZXL(II)  =ZXL(II)  +  (XXL(  II)  *CONJG(XXL  (II)  )  ) 
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ZYL(II)  =ZYL(II)  -••(YYL(II)  *C0NJ3(Y7L  (II)  )  ) 
CSI  (11)  =CSL(II)  +  (XXL(II)  *C0NJ3(rrL  (II))  ) 
3P0  (II)  =XXO(II)+((0.,  1.)*YY0(ir)  ) 
BMO(II)  =XXO(II)  -((0.  ,  1.)  *YYO(II|  | 
BPL  (II)  =XXL(II)-f(  (0.,  1,)  *YYL(ir)  | 
BML(II)  =XXL(II)-((0.  ,  1.)*YYL(II)  ) 
SPOL(II)=SFOL  (II)  ♦  (3P0(II)  *CDNJ3(BPL(II)  )  I 
SMOL(II)  =SMOL(II)  *■  (BJIO(II)  *C3NJ3  (BML  (II)  |  | 
SPC{II)  =SPO(II)  +  (BP3(II)*CONJ3(BPO  (II)  )  ) 
S?L(II)  =SPL(II)■^(BPL(II)  *C0NJ3(BPL  (II)  )  ) 
SaO  (II)  =SaO(II)+  (BMD  (II)  *C0NJ3  (BMO  (II)  )  ) 
S;iL(II)  =S  ML  (I  I)  +  (BML  (II)  *C0NJ3(BML  (II)  )  ) 
3  0       CONTINUE 
70      CONTINUE 
C  THE    NEXT    LOOP    NDRMALIZES    THE       ABOVE 

C  SPECTRA   AND    CALZULATES    POWSR    SPECTRAL 

C  DENSITIES. 

DO    33    13=1  ,N 
ZXO  (13)  =ZX0(I3)  *r/FNR 
ZYO  (13)  =ZY0(I3)  *T/FN3 
CS0(I3)  =CS0(I3)  *T/F:JR 
ZXL(I3)  =ZXL(I3)  *T/FNR 
ZYL  (13)  =ZYL(I3)  *T/FNS 
CSL(I3)  =CSL(I3)  *r/FNR 
SP0L(I3)=SP0L  (13)  *r/FNR 
SPO  (13)  =SP0(I3)  *T/FNa 
SPL(I3)  =SPL(I3)  *r/FNR 
SMO  (13)  =SM0(I3)  *T/FNa 
SML(I3)  =SML(I3)  *T/F;JR 
SM0L(I3)=S  MOL  (13)  *T/FNR 
3  3      CONTINUE 
C  THE    NEXT    LOOP    CALCULATES    STDKES 

C  PARAMETERS    0    THRDUGH    3    OF    THE 

C  INDIVIDUAL    SITE    OHTHDGONAL 


125 


C  COMPONENTS    AND    THE   CDHEREN3E    OF 

C  THE    PLANdR    CIRCaLAR    POLARIZATION 

C  PARAMETERS    BETWEEN    SITES, 

DO    44    14=1  ,N 

CTO  (IU)  =  (ZX0(I4)+ZYD(IU)  I  *2./T 
CP0(I4)  =  (ZX0(I4)-ZYD(I4)I  *2./(r*CTO(I4)  ) 
CR0(I4)  =(4*CS0  (14)  )  /(CT0(I4)  *T) 

COO  (14)  =CSC(I4)/(CS2RT  (ZXO  (14)  I  *CSQRT(ZYD  (14)  )  ) 
CW0(I4)  =ATAN2  (AIMA3(:00  (14)  )  ,  Ri  AL  (COO  (14)  )  ) 
COO  (14)  =CSQRT  (COO  (14)  *C0SJG(C03  (14)  )  ) 
CTO  (14)  =4.  3429448*CL3G  (CTO  (14)  t 
cue  (14)  =4,  342  9  448*CL3G(ZX0  (14)  I 
CVO  (14)  =4.  34  29  448*CL3G(ZY0  (14)  i 
CTL(I4)=(ZXL(I4)>ZYL(I4)|  *2./T 
CPL(I4)  =(ZXL(I4)-ZYL  (14))  *2./(r*CTL(I4)  ) 
CRL(I4)  =  (4*CSL  (14)  )  /(CTL(  14)  *rt 

COL  (14)  =3SL(I4)  /(C33aT(ZXL(I4)  I  *CSQRT  (ZYL(I4)  )  ) 
CTL(I4)  =4.  3429448*CLDG(CrL  (14)  | 
CUL  (14)  =4.  34  2  9  448*CL3G  (ZXL  (14)  ) 
CVL  (14)  =4.  342  9448*CLDG(ZyL  (14)  ) 
CWL(I4)  =AT  AN2  (AIMA3(C0L  (14)  )  ,  SEAL  (COL  (14)  )  ) 
COL  (14)  =CSQRT  (C0L(I4)   *C0NJG  (ZOL  (14)  )  ) 
C0POL(I4)  =SP3L  (14)  /(:  SQRT  (SPD  (14)  )  *CSQRT  (3?L  (14)  )  ) 
COPOL  (14)  =CSQRT  (CDPDL  (I4|  *CD^12  (COPOL  (14)  )  ) 
C0K0L(I4)  =saOL  (14)  / (Z  SQRT  (SMD  (14)  )  *CSQRT  (3ML  (14)  )  ) 
COMOL  (14)  =CSQRT  (COMDL  (14)  *C0NJ3  (COMOL  (14)  )  ) 
SPO  (14)  =4.  342  9  448*CLDG  (SPO  (14)  i 
SPL  (14)  =4.  3429448*CL3G(SPL  (14)  I 
SMO  (14)  =4.  342  9  448*CLDG(S!10  (14)  )  ^ 
SHL  (14)  =4.  342  9448*CLDG(S;iL(I4)  i 
44       CONTINUE 

C  VERSATEC    PLOT    OF     CALCULATED    QUANTITIES 

NPTS=10,/DELTA?-H. 

NPTS '     DETERMINES     NUMBER    DF    POINTS 
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I    >T 


NECESSARY    TO   PLDT    THE    0    Id    13    HEETZ    RANGE, 
ITB  (2)=0 
ITE(3)=20 
ITE  (a)  =  10 
ITB(6)=1 
ITB  (7)=1 
ITE(12)  =1 
ETB  (1)=0 
RTB  (2)  =0 
RTB  (3)=ALAB(1) 
EEAD(5,3000)  TITLE 

CALL    DRAW?  (NPTS,FRQ2,  CUOD  ,118,  RIB) 
ITB(6)  =1 
RTB  (3)=ALAB(2) 
READ  (5,  30  0  0)  TITLE 

CALL    DRAW?  (NPTS,FRQ2,  CVOD  ,ITB,  STB) 
ITB  (6)  =1 
RTB  (3)=ALAB(3) 
READ  (5,  30  0  0)  TITLE 

CALL    DRAHP  (N?TS,FRE3,  CPOD  ,  ITB,  STB) 
ITE(6)  =1 
ETB  (3)=ALA3(ti) 
READ  (5,  30  0  0)  TITLE 

CALL    DRAWP  (N?TS,FHE3,  CR03  ,IT3,  ^TB) 
ITB(6)  =1 
RTB(3)=ALAB{5) 
READ  (5,  30  0  0)  TITLE 

CALL    DRAHP  (N?TS,FREQ,  CROO  (2)  ,113,  RTB) 
ITB  (6)  =  1 
RTB(3)  =ALAB(6) 
READ  (5,3000)  TITLE 

CALL    DRAWP  (NPTS,FRQ2,  CTOD , ITB, STB) 
ITB  (6)  =  1 
RT3(3)  =ALAB(7) 
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READ  (5,  30  0  0)  TITLE 

CALL    DRAWP  (NPTS,FRE2,  CWOD  ,IT3,  SIB) 

ITB(6)=1 

RTB  (3)=ALAB(8) 

READ  (5,  30  0  0)  TITLE 

CALL    DRAWP  (NPTS,FRE3,  COOD  ,ITB,  3  IB) 

ITB(6)  =1 

RTB  (3)=ALAB(9) 

READ  (5,  30  0  0)  TITLE 

CALL    DRAWP  (NPTS,FRQ2,  CULL^ITB,  XTB) 

ITB(6)  =1 

RTB  (3)=ALAB(10) 

READ  (5,  30  0  0)  TITLE 

CALL    DRAWP  (NPTS,?RQ2,  CVLL,ITB,  STB) 

ITB(6)  =1 

RTB  (3)=ALAB(1  1) 

READ  (5,  30  0  0)  TITLE 

CALL    DRAWP  (NFTS^FREQ,  C?LL,ITB,  RIB) 

ITB  (6)  =1 

RTE  (3)=ALAB(12) 

READ  (5,  30  0  0)  TITLE 

CALL    DRAWP  (NPTS,FRE2,  CRLL,IIB,  RTB) 

ITB(6)  =1 

RTB  (3)=ALAB(13) 

READ  (5,  30  0  0)  TITLE 

CALL    DRAWP  (NPTS^FREQ,  CRLL(2)  ,118,  RTB) 

ITB(6)  =1 

RTB  (3)=ALAB(14) 

READ  (5,  30  0  0)  TITLE 

CALL    DRAWP  (NPTS,FRQ2,  CTLL,ITB,  RIB) 

ITB(6)=1 

RTB  (3)=ALAB(15) 

READ  (5,  30  0  0)  TITLE 

CALL    DRAWP  (NPTS,FRE3,  CWLL,ITB,  RIB) 
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ITB  (6)  =  1 

RTB  :(3)=ALA3(16) 

READ  (5,  30  0  0)  TITLE 

CALL    DRAWP  (NPTS,FRE2,  COLL  , ITB,  RIB) 

ITB  (6)  =1 

RTB  {3)=ALAB(17) 

READ  (5,  30  0  0)  TITLE 

CALJ.    DRAWP  (N?TS,FRQ2,  SPOD,ITB,  RTB) 

ITB  :6)  =  1 

RTB  (3)=ALAB(18) 

READ  (5,  30  0  0)  TITLE 

CALL    DRAWP  (NPTS,FRQ2,  SPLL,ITB,  STB) 

ITB  [6)  =1 

RTE  (3)=ALAB(19) 

rea:)(5,30  0  0)  title 

CAL."-    DRAWP  (NPTS,FRQ2,  SMOD,  ITB,  STB) 

ITE  :6)  =1 

RTB  (3)=ALAB(20) 

READ  (5,  30  0  0)  TITLE 

CAL;.    DRAWP  (NPTS,FRQ2,  SMLL,  ITB,  RIB) 

ITE  (6)  =1 

RTE  (3)=ALAB(21) 

READ  (5,  30  0  0)  TITLE 

CALL    DRAWP  (NPTS,FREQ,  C0P0LL,ir3,  RTB) 

ITB  (6)  =  1 

RTB(3)  =ALAB(22) 

READ(5,3000)  TITLE 

CALL    DRAWP  (NPTS,FRS3,  COiDLL,  I  TB  ,  RTB) 
3000       rOEMAT(6A8) 

STOP 

END 
/* 

//GD.3YSIN    DD    * 
TEST    POWER    SPECTRAL    DENSITY    OF    SIN    '4  HZ 
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IN    XXO,30    AVGS,(REF    DB    VS    LOG    FREQ) 

TEST    POWr;H    SFECTHAL    DENSIT7    OF    COS    iEZ 

IN    "£10, 30    .rvGS,(REF    DB    75    LOG    FREQ) 

TEST    STOKES    VSTOKES     0    OF    SIN    4HZ 

IN    XXO    AND    CCS    ^UZ   IN    YyO,3D     AV3S 

TEST    STOKZS    2/SrOKES     0    OF    SEN    43  Z 

IN    XXO    AND    CCS    4 1:  Z   IN    YY0,3D     AVGS 

TEST    STOKES    P/STOKES     0    OF    SIN    UHZ 

IN    XXO    A^d    COS    a  HZ   IN    YYO,3  3     AVGS 

TEST    STOKES    0    OF     SIN    4HZ    i:^     XXO    AND 

COS    4HZ   IN    YY0,3C    A7GS,(REF     DB    7S    LOGi 

TESr    PHASE    OF    SIF    4HZ    IN    XXD 

AND    COS    UflZ    TN    Y^O,30    AV3S 

TEST    COHEP.    Or     SIN    UHZ    IN    XXO 

AND    COS    4HZ    IN    Y^O,  30    AV3S 

TEST    PGWIIR    SPECTx:AL    DENSITf    OF    5  IN    '4H2 

IN    XXL,  30    AVG3,(r:EF    DB    VS    LOG    FREQ) 

TEST    PCWIR    SPSCliiAL    DENSIT'f    OF    COS    4H3 

IN    YYL,30    «.VGS.(E:EF    DB    VS    LOG    FREQ) 

TEST    STOKES    VSTOKES     0    OF    31 N    4ciZ 

IN    XXL    AND    COS    UHZ   IN    YYL,3  3     AVOS 

TEST    STUK2S    2/STOKSS     0    OF    SIN    UHZ 

IN    XXL    AND    COS    UHZ   IN    YYL,3  3     AV3S 

TEST    STOKES    3/SrOKZS    0    OF    31 N    UHZ 

IN    XXL    AND    COS    UHZ   IN    YYL,3  3     AV3S 

TEST    STOKES    C     OF     SIN    UHZ    I.^     XXL    AND 

COS    UHZ    ir<    •n'L,30    A7GS,(REF     DB    VS    LOG)  ■ 

TEST    PHASE    OF    SIN    UHZ    IN    XXL 

AND    COS    UHZ    IN    YYL,  30    AV3S 

TEST    COHER    OF    SIN    UHZ    IN    XXL 

AND    COS    UHZ    IN    YYL,30     AV3S 

TEST    RT    CIRC    POLARIZATION    PSD    OF    SIN    4  HZ    IN 

XXO    AND    COS    UHZ    IN   YYO,     (DB    VS    LOG),  33     AVGS 

TESr    RT    CIRC    POLARIZATION    PSD    OF    SIN    4  HZ    IN 


130 


XXL    AND   COS    nilZ    IN    lYL,     (03     VS    LOG), 33     AVGS 
TSSr    LEFT    CIRC    POIAalZAriON    PSD    OF    SI^'    4  HZ 
IN    XXC    AND    COS    U  H  Z-   IN    YJO    (DB    VS    LOG),  3  0    AVGS 
TEST    LEFT    CIEC    POURIZAIION    PSD    OF    SIM    '4  HZ 
IN    XXL    AND    COS    4 H Z   IN    YfL     (DB    V3    LOS), 3D    AVGS 
TSSr    COHEP    RT    CIPC    .?OLARIZAr  ION    OF    3I!J    UHZ 
IN    XXO,XXI    AND    COS   4HZ    IN    YYO,YYL,    30    A7GS 
TEST    COHEE    LEFT    CIRC    P0LARI3 ATIDN    OF    SIN    UHZ 
IN    XXO,XXL    AND    COS   4HZ    IN    YYO,YYL,    33    AV3S 
/* 
// 
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aP?E:5DrX    G 

//JDEB1MB    JOB     (15  92,  012'3'  ,' FI:;KER    S!1C.  1  399  '  ,  CLA33  =  F 

//*«AIN    OF.G^flPGVJI  1.  IS?:? 

//   EXEC    FORTXCLG 

//FDET.SYoIIJ    DD    D  SN^MSS.  S  5f  S3  , '^ON  IMSL.  3  DQRCS  (FOURD  ,DISP=SHR 

//  DD       * 

C  THE    ARRAY     «  IN*     WILL    BE   USED    TO 

C  RECEIVE   THE    dk7A    PASSED    FRD3 

C  THE    saBROUTINE    '  3  SAD'     AND    THEN 

C  TRANSFERRED    TO    THE    APPROPRIATE 

C  XXX    OR    YYY    Ai'aAY. 

INTEGER*2    IN(16i 

COMPLEXES    XX:C(8  192)  ,  YYY  (31  92) 

DIMENSION    FREQ {310k)  ,TIHE(8  192)  , SORK  ( 163  8U) 

DATA    XXX,YYY/1538u*  (3.  0„  0.  0)  / 

DATA    TIMS,?RE^/16354*3.()/ 
C  THE    EGLLOi^II'iG    JSCTTION    READS    THE 

C  FIRST     SIX    SECOiOS     OF    COMPaiER    TAPS 

c  AND  discar:.?  ?his   data. 

ITL=192 

DO    ^"5    JJ  =  1  ,  ITL 

CALL    RD{2Q, IN,2  00,IH23 ,IRR) 
5  5    CONTINUE 

IFRAME=8192 

NE=20 

FNR  =  FLOAT(NR) 

DO    70    L1=1 , NR 
C  THE    NEXT    LOOP    READS    THE    ZOllpaTER 

C  TAPE    aSING    THE    PROGRAM    PSD7IDED 

C  BY    MR.     TIM    STANTDN    OF    NAVAL    POST- 
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C  GRADUATS    SCK03L. 

DO    6  0    JJ=1 , IFRAMB 

CALL    RD(20,  IN,  1000,I3EC,ii:R) 

XXX  (JJ)=IN  (2) 

YYY  (JJ)=IN  (3) 
5  0    CONTINUE 

N=8192 

FN=FLOAr  (M) 

DELTAT=1 ./3  2. 

T=fN*DELTAT 

DELTAF=1./T 

DO    20    J=1,N 

TIME  (J)=DELTAr*FLOAr(J) 

FPEQ  (J)=DELTAF*FLOAT(J) 
C  THE    NEXT    U     STEPS    C0N7ERT    I'-iS 

C  DATA    TO    VOLTASE    AND    £NSUaE3    THAT 

C  NO    ERRONEOUS    DATA    HAS    BEEN    INTRODUCED 

C  INTO    THE    ARRAYS. 

XXX  (J)  =  (XXX  (J)  -204  5.5)  =^iO./204  5.  5 

XXX  (J)  =REAL  (XXX  (J)  ) 

YYY  (J)  =  (YYY  (J)  -20U5,  5)  *10-  /2045.  5 

YYY  (J)  =REAL  (YYY  (J)  ) 
20    CONTINUE 

CALL    FOUR!  ( XXX,N,1,-1,  0,W3RK) 

CALL    FOUR!  (YYY,N,1,  -1,  0,>:3RK) 

DO    40    K4=1 , N 

XXX  (KU)=XXX  (K4)  /FN 

YYY  (K4)=YYY  (KU)  /FN 
'4  0    CONTINUE 
C  THE    NEXT    LOOP    APPLIES    THE    SYSTEM 

C  TRANSFER    FUNCnON    TO    THE    TRANSFORMED 

C  FREQUENCY    DOMAIN    DATA.     THE    TRANSFER 

C  FUNCTION    CONVERTS    VOLTS    TO    NANDTESLAS. 

DO    9    L  =  1  ,  N 
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c 
c 
c 
c 
c 


FSQ  =  FBEQ(L) 

1  IF  (FEQ.LE.  1  5.)  GO    TD    2 

XXX  (L)  =XXX  (L)  /(  J05.5-3  .  m*FRQj 
YYY(L)  =YYY  (L)  /(131.  32-7.53  8*FR2l 
GO    TO    8 

2  IF(FRQ.LE.  10.)  GO   TO    3 

XXX  (L)  =XXX  (L)  /(  5.95  3*FRQ-3  0.97) 
YYY  (L)  =YYY  (L)  /  (7.  166*F  RQ-3  9  .  99) 
GO    TO    8 

3  IF  (FEQ.LE. 7. 5)  GO    TD    '4 

XXX  (L)  =XXX  (I)  /(3,'49  2*?RQ-5.31) 
YYY(L)  =YYY  (L)  /  (U.  25  2*  FRQ-1  0  .  35) 
GO    TO    8 

4  IF(FRQ.LE.  5.)  GO    TO    5 

XXX  (L)  =XXX  (L)  /(2.6311«'FEQ+0.  UG57) 
YYY  (L)  =YYY  (L)  /  (3  .01  2*P  EQ-1  .  55j 
GO    TO    8 

5  IF     (FRQ.LE,  3.)  GO    13    5 

XXX  (L)  =X  XX  ( L)  /  (  2 .  6  3  11  *  FF  Q<-  0 ,  1  '4  i  S  7  J 

7  YYY  (L)  =YYY  (L) /(2.702*FRQ) 
GO    TO    8 

6  XXX  (L)  =XXX  (L) /(2.72*FaQ) 
GO    TO    7 

8  CONTINUE 

9  CONTINUE 

THE    NEXT    WRITE    STATE-IEIITS    SEND 
THE    CONVERTED    DATA    TO    y.SS 
FOR    FUTURE    MANIPJLATION    AIID     RECALL. 
THIS    NEXT    SET    OF    STAIEMENId    ALSO 
PROVIDES    THE    US2S     A    DIAGNOSnC. 

HRITE(21) XXX 

WRITE  (6,*)  XXX  (1)  ,XXX(3  192) 

WRITE(21) YYY 

WRITE  (6,*)  YYY  (1)  ,yyY(3  192) 
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7  0  CONTINUE 

ENDFILS  21 

STOP 

END 
SUBROUTINE    RD  (lUN^ID,  IRS,  IREC,r.RQ} 
C 

C  THIS    PROCEDURE    FJRNISflSD    3Y    ^B.     TIM    STANTON, 

C  DEPARTMENT    OF    3C3A NOGR APHY . 

C 

C  READ    DATA    FROM    UN,    ALLI3N    ,    CHECK    5    RETURN 

C 

C  IUN=TA?E     NUW3ER,    E3    20 

C  IC=INTSGER*2    ARRAY,      16    LONG, 

C  (VALUES    0-4095,    SUBPRACr    20431*5/2028.    GIVES    VOLTAG] 

C  IES=    NUMBER    OF    RSSIN'CS    ALLOWED     (ERRORS) 

C  IREC=    COUNTER    OF    R3 CORDS     (FRA.1ES    OF    DATA} 

C  BLOCK    512    BITS,  32     BUS    =     REC05D 

C  300    3PI    TAPE    U^ILABLSD 

C  IRQ=    NUMBER    0?    ACTJAL    RESINZS     (ERRORS) 

C 

c 

INTEGER    *     2   10  (16)  ,1?  (16| 

DATA    IRR    /O/ 

IF     (IREC.EQ.O)     IS  =  0 

IER=0 
20  FORMAT     (16A2) 

IF     (IS.  NE.  0)     GO    TO    53 

READ     (lUN,  20,END=900|     IP 

IEEC=IREC<-1 
40  IS=IS+1 

IF     (IS.  LT.  17)     GO    ro    5  0 

READ     (lUN,  20,END=900)     IP 
IS=1 
IPEC=IRSC*1 
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50 
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ICH=IMASK(IP(IS)  ,3,0|  <•  1 
WRITE     (6,55)     ICH,IS,IUN,IEEC 
FORMAT     (•        RESYNCIN3    ICH,  IS,  lU:^ ,  IREC    ' 


UI3) 


IF     (ICH.  NE.  1)     GO    T3    '4  3 
DO    100    1  =  1 ,  16 

10(1)  =ISHIFT(TF(IS)  ,4) 
ICK  =  i:iASK(IP(IS)  ,3,0)  f  1 
I?     (ICH,  SQ.  I)     GO    ID    3D 
IER=IER+1 

WRITE     (6,70)        IUN,IRE:     ,I,  ICH,IER 
70  FORMAT     (' U  NIT'  ,13,  •  R3  CORD  '  ,15  ,  *  39  AN    5    DATA    C;i    ',2IU, 

$     'ERRORS    ',17) 
80  IS=IS+1 

IF     (IS.LT.  17)     GO    ro    1  00 
READ     (lUN,  20,  SND=903|      I? 
IS=1 

IREC=IREC+  1 
100  CONTINUE 

IF     (lER.ZQ  .0)      GO    rO    1  50 

IHR=IRR+1 

IF     (IRR.LT.IRS)     GD    TD     123 

WRITS     (6,  1  10) 
113  FOR^.AT     ('1    STOPPED    IN     SUB    RD    BECAUSE 

C    OF    IRR.GT.  ',16,'     AT    L110') 

IRQ=IRR 

STOP 
120    CONTINUE 

WRITE  (6,130)  IREC,IRR 
133    FORMAT  ('RESYNC  AT  ?  S  AME'  ,16  ,  '  "-JIIH  TOTAL  ERRORS', 17) 

IER=0 

IRQ=IRR 

GO  TO  50 
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153 

90D 
910 


C 
C 
C 


30 


C 

C 


50 


CONTINUE 

RETURN 

WRITE     (6,910)      lUN^IREC 

FORMAT     ('1  END    OF    UNIT    ',13,'     AT    EEC    ',17) 

STOP 

END 

FUNCTION    ISHIFT     (IN,  NPLQ  = 

RETURNS    SHIFTED    VALUE    OF    1*2    WORD    IN 

-VE    LEFT, +72    RI3HT    SHIFT 

INTEGER    *     2   IN 

IP=IN 

I?     (IP.  LT.O)     I?=IP-«-55  536 

IF     (NPLC.LT.O)     GO    TD     30 

ISHIF?=IP/  (2**IA3S  (N?  LC)  ! 

RETURN 

ISHIFT=I?*  (2**IABS  (N?LC)) 

IF     (ISHIFT  .ST. 65535)      ISHI  FT=!1D0  (ISHIFT  ,.65  5  36  ) 

RETURN 

END 

FUNCTION    I  MASK     (IN,I3L,I3R) 

aASX    1*2    WDRD    IN    OUTSIDE     31x3    I3L    o    I3R 

INTEGER    *     2   IN, 10 

IO=IN 

IF     (I3R.EQ  .0)     GO    TO    5  0 

IT  =  ISHIFT  (IN,I3R) 

IO  =  IT 

IP=ISHIFr  (I0,I3L-15-I3R) 

I0  =  IP 

IMASK  =  ISHIFT(I0,15-IBL) 

RETURN 

END 


/* 
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//GD.FT2  1F00  1  DD     UNIT=333 07, MSV3 P=PU34 A , DISP= (NE« , CAILG) , 

//  DSN=MSS.S1592,  GMDri A, 

//  DCB  =  (HECFM  =  VB3 ,BLXSIZE  =  4096,LRECL=4092)  , 

//  SPACE=  (CYL,  (8,  8)  ) 

//GD.FT20F001  DD    UNIT=343  0-'4  ,  VOL  =SER  =  3aDr  1  A,  DISP=  (OLD,  ?  ASS) 

//  LABELS  (1,NL,,r  N)  , 

//  DCB  =  (R2CFM  =  FB,  LR2CL  =  32 , BLKSIZE  =  5 12 , DEN=2) 

// 
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APPENDIX    H 


MAIN     PR3GRAM 


//J3EK1M3   JOB    (1592,  0129)  ,' F  ISHER    SMC.  1  399  •  ,  CLASS  =  F 
//*MAIN    ORG=N?GVai.  1  5 92P,  LI.^  ES=  ( 6 0) 
//♦FORMAT    PR,DDNAME  =  PLOT.  SYS  VECr P, DESr  =  LOCAL 

//    EXEC    FETXCLGP,  FARM. LKED  =  '  LISr, MAP, XREF*  , REGION.  GO=2'^00K 
//SYSUT1  DD      UNIT  =  SYSDA,SPACE=  (CYL,  (8,8)  ) 

//SYS0T2         DD      aNIT=SYSDA,3?ACE=  (CYL,  (3,8) ) 
//SYSLIN         DD      SPACE=  (6030,  (80,80)  )  ,a^ir  =  SYSDA 
//FDRT.SYSIN    DD    DSN  =  MSS.SYS3  .  NONI  MSL.  3  DIJRCE  (FOURI)  ,  DISP=SHR 
//  DD       * 

C  DETAILED    DESCRI  ?  TIO^J  S    OF    THE 

C  VARIOUS    STEPS    AND    LOOPS    ARE 

C  CONTAINED    IN    APPENDICES    B    AND    C. 

INTEGER*2     IN(16) 
COMPLEXES     XXO  (8192)  ,YY0(3  192)  , 

C    CSO(8192)  ,  COO  (8192)  ,ZXO  (8192)  , 

C    CRC  (8192)  ,  CrO  (3192)  ,  COO  (3  1921  , 

C    CVO  (819  2)  ,CWO  (8192)  ,CPO  (3  192)  , 

C    ZYO(8192)  ,C0MOL(8192)  ,C0L(8192|  , 

C    CUL  (8192)  ,CVL  (3192)  ,:WL(3192)  , 

C    COPOL  (3192)  , XXX  (8192)   ,YYY(8192|  , 

C    XXL  (819  2)  ,  YYL(8192)  ,Z  XL  (3  19  2)  , 

C    ZYL(8192)  ,  CSL  (8192)  ,  BPO  (3  1921  , 

C    BMO  (819  2)  ,  EPL  (8192)  ,3ML  (3  192)  , 

C    SPOL(8192)  ,SMOL(3192)   ,SFD  (8192)  , 

C    SPL(819  2)  ,  SMO  (8192)  ,CTL(3  19  2|  , 

C    C?L(8192)  ,CRL  (8192)  ,SML(3192) 
DIMENSION    TIME  (8192)  ,  FREQ  (81921  , 

C    WOEK(16384) ,?RQ2(8192) 
INTEGER*U     ITB  (12)  /12*0/ 
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REAL* 4  BTB(28)  /23*  0.  0/,  CR33  ( 1  6  38  4)  , 

C    C?00(16384)  ,CUOO(16334)  , 
C    COOO(16  38  4),CTOO(163  3  4)  ,3VOO(15  384)  , 
C    CWOO(16334)  ,CTLL(16334)  , 
C    CPLL(16384),3RLL(16334)  ,:OLL(15384)  , 
C    CDLL(16384),CVLL{16334)  , 

C    CWLL(16384) ,COMOLL (15 384)  , COPOLL ( 16384)  , 
C    SFOO( 16384), SPLL(16334)  , 
C    SHOO  (16384), SMLL(1633  4) 

REAL    ALA3 (22) /«PDXO» , »PDyO« ,'SriO' , 
C    •ST20', •ST30« , 'SrOO* , 'PHAO' , 
C    'GAiO' , «PDXL» , •?DIL» , •ST1L','Sr2L«, 
C    •ST3L', 'STOL*  ,  «?HAL' ,  'GAaL'  , 
C    'SPLO',  'S?IL'  ,  'SMDO*  ,  •  S  MLL  •  ,  '  3?  DL  '  ,  'CMOL'  / 

REAL*8    TITLE  (12) 

EQDIVALSNCE(TITLE(1)  ,  RTB  (5  )  )  ,  (3RD0  ( 1 )  , 
C    CR0(1))  ,  (CE00(1)  ,CP3(1)  ), 

c    (ccoo  (1)  ,coD(i) ) ,  (crDO(i)  ,CT3  (1 ) )  ,  (cuoo  (1)  ,cao(i)  )  , 

C     (CWOO  (1)  ,CW0(1)  )  ,  (C7D0(1)  ,C7D  (1  )  )  ,  (CT  LL  ( 1 )  ,  CTL  (  1 )  )  , 

C     (CPLL  (1)  ,CPL(1)  )  ,  (CRLL(1)  ,C?.L(1)  )  ,  (COLL  (1)  ,COL  (1)  )  , 

C     (CULL  (1)  ,CUL(1)  )  ,  (CVLL(l)  ,CVL(1  |  )  ,  (CWLL(1)  ,CWL(1)  )  , 

C     (COPOLL  (1)  ,C0P0L(1)  I  ,   (C0J10LL(1|  , 

C    COJ^OL  (1)  )  ,   (S?00(1)  ,3?0  (1»  )  , 

C     (SPLL  (1)  ,SPL(1)  )  ,  (S}1D0(1)  ,SMD  (1)  )  ,  (SMLL(1)  ,SML(1)  ) 

DATA    XX0,YY0,CS0,C03, ZXO,ZYO,CRD, 
C    CTO,CUO,cyC,CWO,XXX,f YY, 
C    XXL,YYL,ZXL,3YL,CSL, BPO,BMO,BPL, 
C    BML,S?OL,S  K0L,SP0,S?L,SM3, 
C    SML,CTL,CRI,COL,C7L,COPOL,CWL,:3rOL, 
C    COMOL/2949  12*  (0.0,0.  0)  / 

DATA    TIME, FREQ,FR22/24575*0.D/ 

ITL=192 

DO  55  JJ=1  ,ITL 

CALL  RD  (20,IN,200,IREC,IRR) 
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55       CONTINUE 

IFRAM2=819  2 

NR  =  20 

FNE=FL3Ar (NH) 
C  THE    FIRST    SPAZEMENrS    OF    LD3P    70 

C  RECALLS    DATA    PREVIOUSLY    STORED    IN 

C  THE     IBM    3033    iASS    ST0RA3E    SYSTEM. 

DO    70    L1=1  ,NR 

READ(21)XXX 

READ(21)  YYY 

DO    H3    11=1  ,IFRAME 

XXL(II)  =XXX(II) 

YYL(II)  =YYY(II) 
14  3      CONTINUE 

DO    60    JJ=1  ,IFRAME 

CALL    RD  (20  ,IN,  1000,r3EC,IRR) 

XXO(JJ)  =IN  (2) 

YYO(J.J)  =IN  (3) 
5  0       CONTINUE 

N=8192 

FN=FLOAT(N) 

DELTA  T=1./ 32. 

T=FN*DSLrAT 

DELTAF=1./T 

DO    20    J=1 , N 

TIMS  (J)  =DELrAT*FL3Ar  (  J) 

FREQCJ)  =DELrAF*FLOAr  (J) 

FRQ2  (J)  =AL03lO  (FREQ(J)  ) 

XXO  (J)  =  (XXO(J)  -204  5.5)*5./20^5.  5 

XXO(J)  =REAL(XXO(J|  ) 

YYO  (J)  =  (YYO(J)  -20'45.5)  *5.  /2045.  5 

YYO(J)  =REAL(YYO(J)  ) 
20       CONTINUE 

CALL    FOURT  (XX0,N,1  ,-1  ,0,wORK) 


1U1 


CALL    FOUET  (YY0,N,1  ,-1  ,0,»ORK) 
DO    aO    KU=1 ,N 
XXO  (KU)  =XXC(K4)  /FN 
YYO  (K4)  =YY0(K4)  /FN 
4  0       CONTINUE 

DO    9    L=1,N 
FRQ=FREQ(L) 

1  IF  (FRQ.LE.  15.)  GO    TO    2 

XXO(L)=XXO  (L)  /  (-81.53  +17.065*FR2  -  0  .  5389  *FRQ*?RQ) 
YYO  (L)  =YYO  (L)  /(-6  3.'4  3  5  +1  1  .  71  3  3  *FRQ-0- 23279*FRQ*FF.Q) 
GO  TO  8 

2  IF  (FRQ.LE.  10.)  GO  TO  3 

XXO  (L)  =XXO  (L)  /(-13.436  ■••  U.9029*FRQ) 

YYC(L)  =YYO  (L)  /  (-63.435  +11.7105*FRQ  -0.  23279  *FFQ*FRQ) 

GO  TO  8 

3  IF(FRQ.LE.  7.5)  GO  TO  4 

XXO  (L)=XXO  (L)  /  (3.947*  FRQ-9.358) 

YYO  (L)=YYO  (L)  /  (4.  2  95*  FRQ- 1  1 .  33  7  | 

GO  TO  8 
U   IF  (FRQ.LE.  5.)  GO  ID  5 

XXC  (L)  =XXO  (L)  /(2.8105  *FRQ-0.5557) 

YYO(L)  =YYO  (L)  /  (2.3  03*  FRQ-.50) 

GO  TO  3 
5   XXO(L)  =XXO  (L) /(2.723*FRQ| 

YYO  (L)  =YYO  (L)  /  (2.  7  17*  FRQ) 

GO  TO  8 

8  CONTINUE 

9  CONTINUE 

DO  30  11=1  ,N 

ZXC  (II)  =ZXO(II)  +  (XX3(II)  *CONJG  (XXO  (II)  ) 
ZYO(II)  =ZYO(II)  +  (JYD  (II)  *CONJG(£YO  (II)  ) 
CSO  (II)  =CSO(II)+  (XXD  (II)  *C0NJ3  (YYO  (II)  ) 
ZXL(II)  =ZXL(II)  +  (XXL  (II)  *CONJG  (XXL  (II)  ) 
ZYL(II)  =ZYL(II)  -^(lfYL(II)  *CONJG(ryL  (II)  ) 
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CSL  (II)  =CSL(II)  +  (XXL  ( 
EPO  (II)  =XXO(II)  +  ((0.  , 
3M0(II)  =XXO(II)  -(  (0.  , 
BPL(II)  =XXL(II)  +((0., 
BML  (II)  =XXL(II)-(  (0.  , 
SPOL(II)  =SFOL  (II)  +  (BP 
Si^OL(II)=SMOL  (II)  ♦  (33 
SPC  (II)  =SPO(II)  +  (BPD  ( 
S?L(II)  =SPL(II)  ■f(BPL( 
SaO  (II)=SMO(II)  +  (BMD( 
SML(II)  =SML(II)-^  (3ML( 

30      CONTINOE 

7  0       CONTINUE 

DO  33  13=1  ,N 
ZXO  (13)  =ZX0(I3)  *T/FNR 
ZY0(I3)  =ZYG(I3)  *T/FMR 
CSC  (13)  =CS0(I3)  *r/F^is 
ZXL(I3)  =ZXL(I3)  *T/FMS 
ZYL(I3)  =ZYL(I3)  *I/?liR 

CSL (13) =csL(i3) *r/F^a 

SP0L(I3)  =SFOL  (13)  *T/F 
S?0  (13)  =S?0(I3)  *r/F!^El 
SPL  (13)  =SPL(I3)  *T/F:^R 
SM0(I3)  =SMC(I3)  *T/FNa 
SML  (I3)=S:1L(I3)  *r/FMS 
SM0L(I3)  =S-^OL  (13)  *T/F 
3  3       CONTINUE 

DO  44  14=1  ,N 
CTO  (14)  =  (ZX0(I4)  +ZYD( 
CPO  (14)  =  (ZXO  (I4)-ZY3( 
CR0(I4)  =(4*CS0  (14)  )  /( 
COO  (14)  =CSO(I4)/(CS2a 
CW0(I4)  =ATAN2  (AIMAS(: 
COO  (14)  =CSQRr  (COO  (14) 


II)  *CONJS(YYL  (II)  )  ) 

1.) *YYO(II) I 

1.) " YYO(ir) ) 

1.)  *YYL(II)  I 

1.) *YYL(ir)  > 

0(11)  *C3NJ3 (BPL  (II) )  I 

0(11)  *C:NJ3 (BML  (II) )  I 

II)  *C0NJ3 (BPO  (II) )  ) 

II) *CONJG (3FL  (II) )  ) 

II)  *CDNJ3  (BMD  (II)  )  ) 

II) *C0NJ3 (BML  (II) ) ) 


NR 


NP. 


14)  )  *2./r 

14)  I  *2./(r*CT0  (14)  ) 

CT0(I4)«'T) 

T  (ZXO  (14 1  I  *CSQRT(ZYD  (14)  )  ) 

00  (14)  )  ,R5AL(C00  (14)  )  ) 

*COS JG (COD  (14) )  ) 
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CT0(I4)  =4.  342  9aa8*CL3G(CT0  (14)  I 
cue  (14)  =4.  3429448*CL3G  (ZXO  (14)  ) 
CVO  (14)  =4.  3429448*CLDG(ZY0  (14)  ) 
CTL  (14)  =  (ZXL(I4)+3YL  (14)  I  *2./r 
CPL  (I4)  =  (ZXL(I4)-ZYL(I4)|  *2./(r*CTL(I4)  ) 
CEL(I4)  =(4*CSL  (14)  )  /(CTL(I4)  *T| 

COI  (14)  =CSL(I4)/(CS3RT  (ZXL(I4)  i  *CSQRT  (ZYL  (I  4)  )  ) 
CTV  (14)  =4.  3429448*CLDG(CTL  (14)  ) 
CUI(I4)  =4.  3429448*CL3G(ZXL  (14)  ) 
CVL  (14)  =4  .  3429448*CL3G  (ZZL  (14)  ) 
Ca'L(I4)  =ATAN2  (AIMAG(:0L  (14)  )  ,?5AL  (COL  (14)  )  ) 
COL  (14)  =CSQRr  (C0L(I4)  *C0.^  JG(COL.  (14)  )  ) 
CCPOL  (14)  =SPOL  (14)  /(:  SQET  (SPD  (14)  )  *CSQRT(SPL  (14)  )  ) 
COFOL  (14)  =CSQRT(CDPDL  (14)  *  CON  J  5  (COPOL  (14)  )  ) 
C0i!':0L(I4)  =SflOL  (14)  /{ZSQRT  (S.^D  (14)  )  *CSQET(SaL  (14)  )  ) 
CO^JOL  (14)  =CSQRT  (CDMDL  (I4|  *C0SJ3  (COMOL  (14)  )  ) 
SFC(I4)  =4.  342  9448*CL0G(S?3  (14)  I 
SPL  (14)  =4.  3429448*CLDG(SPL  (I'4)  i 
SMC  (14)  =4.  3429448*CLDG  {S'AO  (14)  ) 
5:hL(I4)  =4.  3429448*CLDG(S)!L  (14)  I 
4  4       COl'TINUS 

NPTS=10./DELrAF  +  1. 

IT?  (2)  =0 

ITS  (3)=23 

ITE(4)  =  10 

ITS  (6)  =  1 

IT3(7)=1 

ITS  (12)  =1 

RT5  (1)  =0 

RTB  (2)=0 

RTB(3)  =ALAB(1) 

READ  (5,3000)  TITLE 

CALL    DRAWP  (NPT  S  ,FRQ2  ,  CU03  ,IT3,Rr3) 

ITB  (6)  =  1 
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ETB  {3)-iLAB(2) 

READ  (5,  30  0  0)  T"  TLB 

CALL    DRA't^P  (N?::S,FRQ2,  CVOD  ,IT3,  RIB) 

IT5(6)  =1 

RT5  {3)=ALAB(3) 

EEAD(5,3C00)  T!!  TIE 

CALL    DRAW?  (NP::S,FRE2,  CPOD  ,ITB,  SPB) 

IT3  (6)  ^'! 

aTB  (3)=ALAB{U) 

HEAD  (5,  30  0  0)  T"  TLE 

CALL    DRAW?  (NP::S,FRSQ,  CROD  ,113,  STB) 

ITB  (6)  =  1 

RTE(3)  ^ALAB(5] 

HEAD  {5,30  0  0)  TITLE 

CALL    DRA^?P  (NPTS,FRE3,  CPOD  (2)  ,113,518) 

ITB  (6)  =  1 

RTB  (3)  ==ALAB(6' 

READ  (5,  30  0  0)  TETLE 

CALL    DRAWP  {NPTS,FRQ2,  CTC3  ,ITB,RrB) 

ITB  (6)=1 

RT5(3)  =ALAB(7) 

READ  (5,  3'^  0  0)  TITLE 

CALL    D:%AWP  (NPrS,?RE3,  CWOD  , ITB,  RIB) 

ITB  (5)  =  1 

RTB(3)  =AL;^B(8) 

HEAD  (5,  3J  DO)  TITLE 

CALL    DRAW?  (N?TS,FRE2, C003  ,IT3,RrB) 

ITB  (6)  =  1 

RTE(3)  =ALAB(9) 

READ(5,3000)  TITLE 

CALL  DRAWP  (NPTS,FRQ2, CULL  , ITB, R IB) 

ITB  (6)  =  1 

RTB  (3)  =ALAB(10) 

READ(5,3000)  TITLE 
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CALL    DRAvJP  (N?TS,FRQ2,  CVLL  ,IT3,  S  IB) 

ITB  (6)^1 

RTB(3)=ALAB(1  n 

READ  (5,  30  GO)  TITLE 

CALL    DRAWP  (NPTS^FRE^, CPLL ,ITB, RTB) 

ITB  (6)  =  1 

HTB  (3)  =ALAB(12) 

READ  (5,  300  0)  TITLE 

CALL  DRAW?  (NFTS^FRE;),  CRLL  , ITB,  STB) 

IT3(6)=1 

RTB(3)  =ALA3{13) 

READ  (5,  30  0  0)  TITLE 

CALL  DRAWP  (NPTS  ,  FRE3,  CRLL  (2)  ^IIB^RTB) 

ITB(6)  =1 

RTB  (3)=ALAB(14} 

READ  (5,  30  0  0)  TITLE 

CALL    DRAHP  (N?T:J,FRQ2,  CTLL,ir3,  RIB) 

ITB  (6) =1 

RTB(3)=ALA3(1:5) 

READ  (5,  30  0  0)  TITLE 

CALL    DRAW?  (N?T3,FREi,  CWLL,ITB,  RIB) 

ITB  (6)  =1 

RTB  (3)=ALAE('I6) 

READ  (5,  30  0  0)  TITLE 

CALL    DRAW?  (NPTS,  FRE2,  COLL, 113,  STB) 

ITB  (6) =1 

HTB(2)=ALAB(1'?) 

READ  (5,  30  0  0)  TITLE 

CALL    DRAW?  (N?TS,FRQ2,  SP0D,IT3,  RIB) 

ITB(6)  =1 

RTB  (3)=ALAB(18) 

READ(5,3000)  TITLE 

CALL    DRAWP  (N?TS,FRQ2,  SPLL,ITB,  RTS) 

ITB  (6)  =1 
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c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 


RTB(3)=ALAB(19) 
READ(5,30aO)  TITLS 

CALL    DRAW?  (NPTS..FBQ2,  s:MOD,irB,arB) 
ITB(6)  =1 
RTB(3)=ALAB(20) 
READ  (5,  30  0  0)  TITLE 

CALL    DRAW?  (NPTS,FHQ2,  SMLL,  ITB,  5  T  3) 
ITB(6)  =1 
RTB(3)=ALAB(21) 
READ  (5,  30  0  0)  TITLS 

CALL    DRAWP  (NPTS,?REQ,  COPOLL,irB, RTB) 
ITB  (6)  =1 
RTB  (3)  =ALAB(22) 
READ(5,3000)  TITLE 

CALL    DRAW?  (NPTS,FRE3, COMOI L, ITB , RTB) 
3030       FORMAT (6A8) 
STCP 
END 
SUBROOTINE    RD (lUN, 13 , IRS, IREC, E RQ) 

THIS    PROCEDURE    FUS UISHZD    3Y    MR.     TIM    STANTON, 
DEPARTMENT    0?    GCEft  HOGRAPHlf . 

READ    DATA    123  A    UK,     ALLISN    ,     CHECK    &    RETURN 

IDN=^TAPE     NUMBER,     E3     20 

IO  =  INTE3SE-*2    ARRA.Y,      16    L0N3, 

(VALUES    0-4095,    SUBIHACr    20'4  3 1  *  5 /2028.     3IVES    VOLTAGE 

IES=    NUMBER    OF    RESr^ICS    ALLOWED     (ERRORS) 

IREC=    COUNTER    OF    RECORDS     (FRAIIES    OF    DATA) 

BLOCK    512    BITS,  32     BITS    =    RECORD 

800    3PI    TAPE    UtTLABLED 
IRQ=    NUMBER    OF    ACT3AL    RESINC3     (ERRORS) 
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INTEGER    *     2   10  ( 16i  ,  I?  C  6)  ■ 

DATA    IP.R    /O/ 

IF     (IREC.  EQ-0)     13  =  0 

IEB=0 
20  FORMAT     (16A2) 

IF     (IS.NE.  C)     GO    TD    53 

READ     (lUN,  20,  SND=903|     I? 

IREC=IREC+ 1 
40  IS=IS+1 

IF     (IS.LT.  17)     GO    TO    5  0 

READ     (lUN,  20,  SJID=903|      I? 
IS=1 

IREC=IREC+1 
50  ICH=IMASK(I?(I3)  ,3,0»  «■  1 

C  WRITE     (6,5  5)     ICK,  IS,  I  UN,I  EEC 

55  FORMAT     (•       RESYNCIN7j     IC  H,  IS  ,  IQ^I ,  IREC    '^UIS) 


C 


IF     (ICH.NE.  1)     GO    TD     n 
DO    100    1  =  1  ,  16 

10  (I)  =1  SHIFT  (IP  (IS)  ,4  ) 
ICH=IM ASK  (IP (13)  ,3,0»  ^  1 
IF     (ICH.E2.  I)     GO    TD    3D 
IEP=IEH+1 

WRITS     (6,70)       inN,IRE:     ,I,IC:i,IEEi 
70  FORMAT     (»  3  I^ilT  '  ,  13,  *  RE  CORD  '  ,  J  5  , '  CHAN    &    DATA    CH    ',214, 

$    'ERRORS    '  ,17) 
80  IS=IS+1 

IF     (13.  LT.  17)     GO    10    1  00 
READ     (ION,  20,END=90D|     IP 
IS=1 

IREC=IRSC+  1 
103  CONTINUE 


C 
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c 
c 
c 


IF     (lER.EQ.O)     GO    TO    1 50 

IRIl=IRR-^1 

IF     (IRR.LT.IRS)     3j    T:^     123 

WRITE     (6,  1  10) 
113  FORMAT     (M    STOPPED    I^    SUB    RD    5F,CAUSS 

C    OF    IRR.GT.M6,'    AT    L110M 

IRQ=IRR 

STOP 
123         CONTINUE 

WRITE    (6,130)     IRECrSR 
133  FORMAT     (»aESYNC    AT    FP.  AME«  ,  16 ,  »  »  ITH    TOTAL    ERRORS', 17) 

IEE=0 

I3Q=I3R 

GO  TO  50 
153    CONTINUE 

RETURN 
903    WRITE  (6,9  10)  IUN,IR2C 
913    FORMAT  (»  1    END  OF  J  NIT  ♦,I3,'  AT  ?.SC  ',17) 

STOP 

END 

FUNCTION    I  SHI  FT    (IN,N?LC) 

RETURNS    SHIFTED    VALQE    D?    1*2    WORD    IN 
-VE   lsft,+ve:    RISHT   SHTF? 

INTEGER    *     2   IN 
IP=IN 

IF     (IP.LT.O)     IP=IP  +  65536 
IF     (NPLC.LT.3)     GO    TD     30 
ISHIFT=IP/  (2**IABS  ('JPLC)  ) 
RETURN 
30  ISHIFT=IP*  (2**IABS  (NPLC)) 

IF     (ISHIFT.GT. 655351      ISHIFT=M3D (ISHIFT, 65 5 36 ) 

RETURN 

END 
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FUNCTION    IMASK     (IN,I3L,I3H) 
C  MASK    1*2    WDRD    IN'    OUTSIDE    BITS    I3L    &    IBR 

C 

INTEGER    *     2   IN, 10 

IO  =  IN 

IF     (IBR.EQ.O)     GO    TO    50 

IT  =  ISHIFT  (IN, IBR) 

IO=IT 
50  IP=ISHIFT (I3,IBL-15-r BR) 

IO=IP 

IMASK=ISHIFT(I0,15-rBL) 

RETURN 

END 

/* 

//GD.FT20F001    DD     UNIT  =  343  0-4  ,  VOL  =SSR  =  3  MDT  1  3  ,  DISP'-^  (OLD,  PASS)  , 

//  LABEL=  (1,NL,,rN)  , 

//  DCB  =  (RECFM  =  FB,  LRS:  L  =  32  ,  3  LKSIZS=5  1;. .  DEN^2) 

//GD.FT21F001    DD    UNIT=3330V,  MS V3P  =  ?aa4  A,  DIS?=  (OLO,  KZ  iiP)  , 

//  DSN=MSS.S1592.  GMDZIA, 

//  DCB  =  (RECFM=VB3  ,3LK  SIZE  =  -i  0  96  ,  LRECL^4  3  92) 

//GD.SYSDUMP    DD    SYSOUT=A 

//GD.SYSIN    DD    * 

X   COIL    PSD,     MIRY     BAY  ,  1 122-1  2  54    LOCAL  ,  1  7  A  iJGS2  , 

20    AVGS,    32S/SEC, 5V0LT,     (DB    NT**2    V3    LOG    FREQ) 

Y    COIL    PSD,     HTRY     EA  Y ,  1  122-1  2  54    LOCAL  ,  1  7/\  UG82  , 

20    AVGS,    32S/SSC,  5V0LT,     ( DB    NT*?'2    V3    LOS    FREQ) 

STOKES    NO.    ONE,    MTRY    BAY,     1122-1254 

LOCAL,     17AUG82,     20   AVGS,    32S/SEC,    SVOLT 

STOKES    NO.    TWO,    MTRY    BAY,     1122-1254 

LOCiL,     17AUG82,    20   AVGS,    32S/SE:,    S^OLT 

STOKES    NO.    THREE,     MIRY    BAY,      1122-1254 

LOCAL,     17AUG82,    20   AVGS,    323/SEC,    SSTOLT 

STOKES    NO.     ZERO,     MTRY    5AZ,     1122-1254    LOCAL, 

17A0G82,     20    AVGS,    32S/SEC,    5  VOLT ,     (L03     VS    LOG) 
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PHASE    OF    X    AND    Y,    MTRY    BAY    1122-1254 
LOCAL,     17AUG82,    2 0   AVGS,    323/SEC,     S/DuT 
COHER    OF    X    AND    Y    COILS,    MTRY    BAY,     1122-1254 
LOCAL,     17AUG82,    20   AVGS,    32S/SE3,    5yOLr 
X    COIL    PSD,    LA    MESA, 1122-12545L3CAL,17AaG82, 
20    AVGS,    32S/SEC,  107OLT,     (DB    NT**2    VS    LOG    FREQ) 
Y   COIL    PSD,    LA    MESA, 1 122-1254    LOCAL, 1 7 AUG82 , 
20    AVGS,     32S/SEC,  10VOLT,     (08    NT»*2    VS    LOG    FREQ) 
STOKES    NO.    ONE,    LA   MESA,     112201254 
LOOM,     17AUG82,    20  AVGS,    323/SE:,     13V0LT 
STOKES    NO.    TWO,    LA   MESA,     1122-1254 
LOCAL,     17AUG82,    20   AVGS,    32S/SE:,     13VDLT 
STOKES    NO.    THREE,    LA    MESA,     1122-1254 
LOCAL,     17AUGS2,     20   AVGS,    32S/SE3,     13V3LT 
STOKES    NO.     ZERO,     LA    MESA,     1122-1254    LDZf^L, 
17AaG82,     20    AVGS,    32S/SEC,    10VOLT,     (L:)3    VS    LOG) 
PHASE    OF    X    AND    Y,    LA    MESA,     1122-1254 
LOCAL,     17AUG82,    20   AVGS,     323/SEC,     10VDLT 
COHER    OF    X    AND    Y    COILS,    LA    MESA,     1122-1254 
LOCAL,     17AUG82,    2  0   AVGS,     323/SEC,     1DVDLT 
RT    :iRC    POLARIZATION    PSD, MIRY    B\ Y , 1 1 22- 1 254 
LOCAL, 17ADG82, 20     AVGS ,32S/3 H C, 57 OLT, (OB    VS    LOG) 
ET    :iRC    POLARIZATION    PSD,    LA    MESA,     1122-125U 
LOCAL,  17AUG82, 20     AVGS , 32S/3EC, 10  VOLT,  ( DB    VS    LOG) 
LEFT    CIRC    POLARIZATION    PSD,     MTRY    BAY  ,  1  1  2  2- 1  2  5u 
LOCAL, 17AnG82,    20    AVGS ,32S/S EC ,5  VOLT , ( 05    VS    LOG) 
LEFT    CIRC    POLARIZATION    P3D,     LA    MESA,     1122-1254 
LOCAL, 17ADG82, 29     AVGS , 32S/32 C, 10  VOLT,  ( DB    VS    LOG) 
COHER    RT    CIRC    POLARIZATION    MTRY    BAY/L^     MESA, 
1122-1254    LOCAL,     17AaG82,    23     AV3S,    323/SEC 
COHER    LEFT    CIRC    POLARIZATION    MTRY    BAY/LA    MESA, 
1122-1254    LOCAL,      17AaG82,     20    AV3S,    325/SEC 
/* 
// 
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APPENDIX    I 


TYPICAL    COaPCirSR    SENERATED    PLOTS 
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Log  Freqency  [ log  Hz  ] 


X-SCflLL=5.00E-01    UNITS    INCH 
T-SCflLE=2,00E+01     UNITS     INCH 
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